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Abstract. The top quark was discoverd at the CDF and DO experiments in 1995. As the partner of the 
bottom quark its properties within the Standard Model are fully defined. Only the mass is a free parameter. 
The measurement of the top quark mass and the verification of the expected properties have been an 
important topic of experimental top quark physics since. In this review the recent results on top quark 
properties obtained by the Tevatron experiments CDF and DO are summarised. At the advent of the LHC 
special emphasis is given to the basic measurement methods and the dominating systematic uncertainties. 
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1 Introduction 

The aim of particle physics is the understanding of ele- 
mentary particles and their interactions. The current the- 
ory of elementary particle physics, the Standard Model, 
contains twelve different types of fermions which (neglect- 
ing gravity) interact through the gauge bosons of three 
forces |ij2iL3it4j- In addition a scalar particle, the Higgs 
boson, is needed for theoretical consistency |5j . These few 
building blocks explain all experimental results found in 
the context of particle physics, so far. 

Nevertheless, it is believed that the Standard Model is 
only an approximation to a more complete theory. First of 
all the fourth known force, gravity, has withstood all at- 
tempts to be included until now. Furthermore, the Stan- 
dard Model describes several features of the elementary 
particles like the existence of three families of fermions 
or the quantisation of charges, but does not explain these 
properties from underlying principles. Finally, the light- 
ness of the Higgs boson needed to explain the symme- 
try breaking is difficult to maintain in the presence of ex- 
pected corrections from gravity at high scales. This is the 
so called hierarchy problem. 

In addition astrophysical results indicate that the uni- 
verse consists only to a very small fraction of matter de- 
scribed by the Standard Model. Large fractions of dark 
energy and dark matter are needed to describe the obser- 
vations |6jX 8j 9j 1_0_ . Both do not have any correspondence 
in the Standard Model. Also the very small asymmetry be- 
tween matter and anti-matter that results in the observed 
universe built of matter (and not of anti-matter) cannot 
be explained until now. 

It is thus an important task of experimental particle 
physics to test the predictions of the Standard Model to 
the best possible accuracy and to search for deviations 
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pointing to necessary extensions or modifications of our 
current tfieoretical understanding. 

The top quark was predicted to exist by tlie Standard 
Model as tfie partner of the bottom quark. It was first 
observed in 1995 by the Tevatron experiments CDF and 
D0 pilfT^ and was the last of the quarks to be discov- 
ered. As the partner of the bottom quark the top quark 
is expected to have quantum numbers identical to that of 
the other known up-type quarks. Only the mass is a free 
parameter. We now know that it is more than 30 times 
heavier than the next heaviest quark, the bottom quark. 

Thus, within the Standard Model all production and 
decay properties are fully defined. Having the complete 
set of quarks further allows to verify constraints that the 
Standard Model puts on the sum of all quarks or particles. 
This alone is reason enough to experimentally study the 
top quark properties. The high value of the top quark 
mass and its closeness to the electroweak scale has inspired 
people to speculate that the top quark could have a special 
role in the electroweak symmetry breaking, see e.g. [131 
[M] . Confirming the expected properties of the top quark 
experimentally establishes the top quark as we expect it to 
be. Any deviation from the expectations gives hints to new 
physics that may help to solve the outstanding questions. 

Since the last review on top quark physics in this jour- 
nal [TS] the luminosity at the Tevatron has been increased 
by more than an order of magnitude. Now measurements 
of the top quark mass are no longer limited by statistical 
but by systematic uncertainties. Even discussions about 
the knowledge of the implicitly used theoretical mass def- 
inition become relevant. Moreover measurements of addi- 
tional top quark properties now reach a viable precision. 
In this review the recent results on top quark properties 
obtained by the Tevatron experiments CDF and D0 are 
summarised. At the advent of the LHC special empha- 
sis is given to the basic measurement methods and the 
dominating systematic uncertainties. Other reviews with 
different emphases are available in the literature [l6lfT7l 

After a short introduction to the Standard Model and 
the experimental environment in the remainder of this sec- 
tion. Chapter [5] describes the current status of top quark 
mass measurements. Then measurements of interaction 
properties are described in Chapter [3] Finally, Chapter 2] 
deals with analyses that consider hypothetical particles 
beyond the Standard Model in the observed events. 



1.1 Theory 



In the following a short description of the Standard 
Model shall be given to set the stage for later descriptions. 
For this the natural units where h — c = 1 will be used. 



1.1.1 The Standard Model 

Lagrangian The Standard Model is a quantum field the- 
ory with local gauge symmetry under the group S'C/(3) x 
SU{2) X C/(l). Its Lagrangian contains fields correspond- 
ing to three types of particles: Gauge or vector bosons, 
fermions and scalars. Gauge boson are described by vector 
fields, A^, scalars by complex fields, 4>. Fermions can be 
described by Weyl spinors, "01 with left and right handed 
helicity. Using Einstein's summation convention the La- 
grangian can be written as 



C 



1 



pA j^Ap. ^ ^^^0^^ + {D^cp'^)[D'^r) 



-^Vapi^ai^P^*^ 4- V{(t)) + (ghost- and gauge terms) 



(1) 



Here capital Latin letters run over the gauge bosons, lower 
case Latin letters over the scalar fields and Greek letters 
a, /3 index the fermions of the standard model, fi and u 
are Dirac indices. The field tensor is defined as 






d,A^ - gCABcA'^Ai 



(2) 



where Cabc are the structure constants of the gauge 
group. The covariant derivative is defined from the gauge 
symmetry to be 



D, = d,- ^'-^ t^Al 



(3) 



Here gj^j is the coupling strength for the gauge boson 
A^ , i.e. 5, g' or gs- t^ are the generators of the gauge 
symmetry in the representation that corresponds to the 
particle field on which the derivative acts. 



p^ at'Df, with a^' = (l,±cr) 



(4) 



where cr is the vector of Pauli matrices. The positive sign 
applies to fermions of right handed helicity, the negative 
for left handed ones. 

The Yukawa couplings, y^^, are free parameters that 
may be non-zero only when the combination of the 
fermions a and /? with the scalars is gauge invariant. V{(l)) 
is a quartic form. 



The physics of elementary particles is described by the 
Standard Model of particle physics. It describes three of 
the four known forces that act on elementary particles, 
the electromagnetic and the weak force are unified in the 
GSW-theory [TJ[31[3] , the strong force is described by quan- 
tum chromodynamics (QCD) [3]. So far the influence of 
gravitation could not be unified with the other three forces 
in a consistent quantum theory. Due to its weakness it is 
usually safe to neglect its influence in the context of par- 
ticle physics. 



Particle Content The particle content of the Standard 
Model is specified by defining the representations of the 
gauge symmetry for each of the fields contained. 

The Weyl spinors describing the fermions transform 
according to fundamental representations for each of the 
subgroups or may be invariant under a subgroup. Quarks 
transform under the three dimensional representation of 
the SU{3) group, 3. All left handed spinors transform un- 
der SU{2) according to the two dimensional representa- 
tion, 2, while the right handed spinors are singlets, i.e. 
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Table 1. Fermions of the Standard Model and their represen- 
tations. 



be chosen to exist in the lower component of the doublet: 



Gene- 
ration 


(3,2). 


Quarks 

(3,1). 


(3,l)_i 


Leptons 

(l,2)_i (1,1)_, 


1st 


(A 


UR 


dR 


BR 


2nd 


(A 


Cr 


sr 


W. - 


3rd 


r 1 


tR 


bR 


["A 

TR 



Table 2. Gauge bosons of the Standard Model and their rep- 
resentations. 



Symbol 


Representation 


Coupling strength 


g 

{w\w^,w'-') 

B 


(8,l)o 
(l,3)o 

(l.l)o 


9s 
9^ 92 

9 = V^Pi 



invariant under SU{2) rotations. This reflects the left 
handed nature of weak interactions. The transformation 
properties under U{1) are specified by the hypercharge. In 
the Standard Model all representations are repeated three 
times and build the three generation of fermions. 

In Tab. [T] the representations of the fermions of the 
Standard Model are summarised. The representations of 
the non-Abehan gauge groups SU{3) and SU{2) are spec- 
ified by their dimension (in bold- face), the hypercharge 
corresponding to the U{1) group is given as index. For 
a few years it is known from the measurement of neu- 
trino oscillations [25 ] EB l l^7 ] f5] that also the neutrinos have 
mass and thus right handed neutrinos should be added to 
Tab. [TJ The necessary extensions of the Standard Model 
are not unique and thus they are usually not considered 
part of the Standard Model. In the context of top quark 
physics neutrino oscillations do not play any role and can 
thus be ignored for the purpose of this review. 

The gauge bosons of a quantum field theory need to 
transform according to the adjoint representation of their 
sub-group. Thus we get eight gauge bosons for the SU{2>) 
symmetry, the gluons, three for SU{2) and one for the 
C/(l) symmetry, c.f. Tab. [21 In addition to the bosons and 
fermions described the Standard Model contains a com- 
plex scalar doublet, the Higgs doublet ^, which transforms 
according to (1, 2) i . It is needed for symmetry breaking. 



Symmetry breaking: Higgs mechanism In nature the 
symmetry of the Standard Model is broken. The symme- 
try breaking is implemented by the Higgs mechanism [5], 
which assumes that the scalar isospin doublet (p has a vac- 
uum expectation value. This is achieved by proper choice 
of parameters in the most general potential, V{(j)), for the 
scalar field in Eq. ^. According to the symmetry this can 



<P 



with 



m 







(5) 



By expanding the complex scalar field around the vacuum 
expectation value first four real scalar fields are specified. 
Three of these can behave like longitudinal components 
of the SU{2) gauge bosons, W^. Usually the theory is 
thus written in terms of three massive vector bosons W^ , 
Z, W~ , a massless vector boson, the photon A, and the 
fourth real scalar field, the Higgs boson H: 



W^ 



1 



V2 



(W^zpiW^) 



q'B - qW'^ o 

Z = ^ , ^ ^ = sinOw B- cos OwW^ 



A = 



gB + g'W^ 
v/5" + ff' 



= cos 9w B + sin 9w W^ 



(6) 



Here B is the gauge boson of the (hypercharge) U{1) sym- 
metry and the Weinberg angle 9w is defined by the ratio 
of coupling constants 



tan {9w) '■= — 



(7) 



After this rewriting the theory remains SU{3) x C/(l) in- 
variant. The C/(l) symmetry now corresponds to the elec- 
trical charge. So the Standard Model parts are quantum 
chromodynamics (QCD) and quantum electrodynamics 
with their SU{3) and f/(l) symmetries, respectively. 

The Higgs mechanism not only yields massive vector 
bosons, it is also responsible for the masses of the fermions. 
For the specified particle content, the Yukawa terms in 
Eq. (m may be non-zero for left-handed quark or lep- 
ton doublets paired with the corresponding right-handed 
quark and lepton S't/(2)-singlets. For the first generation 
these terms are 



(8) 
Corresponding terms can be written not only for the other 
generations, but in general also for fermion pairs between 
different generations. Unitary rotations in the three di- 
mensional space of generations are commonly used to re- 
define the lepton and quark fields such that Yukawa cou- 
plings occur only between particles of the same generation. 
This provides the mass eigenstates of the quark and lep- 
ton fields. These rotations cancel in most terms of the La- 
grangian. The only observable remainder of this rotation 
is the Cabibbo-Kobayashi-Maskawa (CKM) matrix p8l 
[^IB] which occurs in the coupling of the W^ bosons to 
quarks: 

(Vud Vus Vub \ 
Vcd Vcs Vcb (9) 

Vtd Vts VtbJ 

At the same time this is the only process in the Stan- 
dard Model that connects the different generations. Nu- 
merically this unitary matrix has diagonal entries close to 
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unity and off-diagonal entries that are around 0.2 between 
the first and second generation, around 0.04 between the 
second and third generation and even smaller for the tran- 
sition of the first to the third generation [5] . 



1.1.2 Perturbation Theory 



Predictions of the Standard Model for high energy reac- 
tions are so far generally performed in perturbation the- 
ory. The reactions are described as one or more point- 
like interactions between otherwise free particles. The al- 
lowed reactions can be read off the Lagrangian and are 
usually represented by Feynman diagrams. For example 
the Yukawa term yip(f>'ip yields an interaction vertex of the 
strength y with two fermions ip and the scalar field <j>. 

Each Feynman diagram serves simultaneously as a dia- 
grammatic description of the reaction and as a short hand 
notation for the corresponding computation of the transi- 
tion amplitude. The quantum mechanical amplitude of a 
given process that transforms a set of initial state particles 
to a set of final state particles is given by the sum of all 
possible diagrams with the corresponding initial and final 
state particles as external lines. 

Diagrams with few interactions usually give the largest 
contributions and higher order corrections are suppressed 
by factors of the additional coupling strengths. Thus cal- 
culations are usually performed in a fixed order of the 
coupling constant(s). In some cases, however, kinematic 
enhancements of logarithmic type may compensate the 
suppression by additional powers of the coupling constant. 
Notably these occur in cases of collinear or soft gluon ra- 
diation and for top quark pair production near threshold. 
In these cases resummation of the leading (or next to lead- 
ing) logarithms to all order of the coupling are performed. 

Diagrams of higher orders generally involve loops 
which require to integrate over all possible momenta of 
the internal lines. Naively, such integrals diverge. It is nec- 
essary to renormalise the theory in order to obtain finite 
predictions. There are several possible schemes to perform 
this renormalisation. Most commonly the so called MS 
scheme is used, which itself depends on a continuous pa- 
rameter the renormalisation scale /z. The dependence of 
results on the choice of this parameter is often used as a 
measure for theoretical uncertainties of a prediction. 

Perturbation theory described so far deals with the 
particles of the Standard Model named above, i.e. with 
quarks, leptons, gauge bosons and the Higgs boson. In na- 
ture, however, quarks have not been observed as free par- 
ticles, rather they are confined in bound states of colour 
neutral hadrons. The dynamics of quarks and gluons in- 
side hadrons cannot be described by perturbation theory. 



To describe the collisions of hadrons with perturba- 
tion theory the (soft) physics that governs the behaviour 
of quarks and gluons in the hadron needs to be fac- 
torised from the hard process in the collisions. The par- 
tons inside an incoming hadron are considered as a num- 
ber of "free" partons that may enter the hard interaction. 
The distribution of partons inside the incoming hadron 
is taken from parton distribution functions (PDFs) that 
are derived from experiments. With this the cross-section, 
cr{pp — > X]s), to produce final state particles X can be 
described in terms of cross-sections a of incoming quarks 
and gluons to produce X: 



a{pp -)■ X; s) 



i,j—g,u,u,d,d.. 



da;idx2 fi{xi)Jj{x2) cr{ij -> X; s) (10) 



Here s is the centre-of-mass energy squared of the incom- 
ing hadrons. At the Tevatron Run II y/s = 1.96 TeV. The 
colliding partons carry momentum fractions xi and X2, 
and s = X1X2S is their centre-of-mass energy squared. 
fi and fj are the parton distribution functions for the par- 
tons i and j, i.e. the probabilities to find this parton in the 
proton and anti-proton, respectively. Formally, the par- 
ton density functions /, / and the partonic cross-section, 
(J, depend on the factorisation scale, fxp, which specifies 
which effects are included in the PDFs and which remain 
to be described by the hard matrix element. They also de- 
pend on the renormalisation scale, /i. These dependencies 
are not written out in Eq. ([TU)) . 



1.1.3 Top Quark Production in Hadron Collisions 

In hadron collisions, like the proton-antiproton collisions 
at the Tevatron, top quarks can be produced singly or 
in pairs. The pair production occurs dominantly via the 
strong interaction. In leading order the quark anti-quark 
annihilation and gluon fusion processes shown in Fig. [1] 
contribute. In higher orders also quark-gluon processes ex- 
ist. The relative contribution of these diagrams depend 
on the parton distribution functions. At the Tevatron 
with a centre-of-mass energy of 1.96 TeV next-to-leading 
order predictions lead to an expectation of 85% contri- 
bution from qq annihilation and 15% from gluon fusion. 
The total cross-section of top quark pair production has 
been computed in perturbation theory using various ap- 
proximations 30,31,32.33; 34[[33] . Its value has a signifi- 
cant dependence on the top quark mass, which near the 
world average value is about — 0.2pb/GeV. For a top 
quark (pole) mass of 175 GeV Moch and Uwer [33] find 
ati = 6.90to:64pb, based on the CTEQ6.6 [3B] PDF. The 
uncertainty includes uncertainties of the PDF and the 
scale uncertainty. 

Single top quark production can only take place via the 
weak interaction. The leading processes are quark annihi- 
lation through a W boson, also called s-channel, the quark 
gluon fusion with a W boson in a t-channel, c.f. Fig. [21 and 
production of single top quarks in association with a (close 
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Fig. 1. Born level Feynman diagrams contributing to top quark pair production. The quark annihilation (leftmost diagram) is 
dominating top pair production at the Tevatron. The gluon fusion processes (three right diagrams) contribute about 15%, only. 



b q 




Fig. 2. Leading order Feynman diagrams contributing to sin- 
gle top quark production at the Tevatron. According to the 
structure of the diagrams the left process is called s-channel 
and the right t-channel production. 



to) on-shcU M^-boson. Charge conjugate diagrams apply 
for anti-top quark production. At the Tevatron Run II the 
cross-section to produce a single top or anti-top quark 
is 3.4 ± 0.22 pb. The s- and t-channel contribute a little 
less than 1/3 and 2/3, respectively. The associated pro- 
duction contributes less than 10% ^37u38y39 ]. These num- 
bers were derived assuming rrit — 175 GeV and using the 
MRST2004 gD] PDFs. 

After production top quarks decay very rapidly 
through the weak interaction into a W boson and a h 
quark. In the Standard Model contributions from decays 
to light quarks are suppressed due to the smallness of the 
corresponding entries of the CKM matrix. The expected 
decay width of about 1.34 GeV corresponds to a lifetime 
of the order of 5 • 10~^^s ^B . Thus the top quark decays 
before it can couple to light quarks and form hadrons. 
The lifetime of tt bound states, toponium, is too small, 
Fii ^ 2 Ft, to allow for a proper definition of a bound 
state as already pointed out in the early 1980s [H]. 

The decay modes are defined purely by the W boson 
decays. W bosons may decay to quarks, i.e. hadronically, 
or leptonically. For top quark pairs this yields three basic 
decay modes: the all hadronic channel, the semileptonic 
channel and the dilepton channel. The all hadronic chan- 
nel has a branching fraction of 46%. Each of the three 
charged Standard Model leptons contributes 15% in the 
semileptonic channel. The dileptonic decays have a total 
branching fraction of 9%. See Fig. [3] for a graphical rep- 
resentation. Decays that involve r leptons are usually not 
considered in analyses of the semileptonic and dileptonic 
decay modes, because r leptons are difficult to identify. 
The analyses, however, include the events in which the 
tau decayed to an electron or muon. 



Following this experimental nomenclature the semilep- 
tonic and dileptonic channels are considered to include 
only electrons and muons. If an analysis considers iden- 
tified tau leptons this fact shall be explicitly stated. For 
single top quark analyses so far only the leptonic W boson 
decays (to electrons and muons) are considered. 



1.2 Experiments 

Up to recently only one collider provides centre-of-mass 
energies sufficiently high to produce top quarks: the 
Tevatron at the Fermi National Accelerator Laboratory 
(FNAL) near Chicago, IL, USA. At the two collision 
points typical general purpose experiments of present col- 
lider physics are positioned, CDF and D0. Each consists 
of a cylindrical part that covers particles produced at large 
angles to the beam and two end-caps that detect particles 
at smaller angles. Close to the beam tracking and ver- 
tex reconstruction components are placed, followed to the 
outside by calorimetry and muon detection systems. Some 
more details of the Tevatron and the two detectors shall 
be described below. 



1.2.1 The Tevatron 

The Tevatron collides beams of protons and antipro- 
tons. The protons and antiprotons are produced and pre- 
accelerated in a series of smaller machines and then filled 
into the Tevatron to circulate in opposite directions. In 
a ffi'st phase of operation between 1992 and 1996 the 



alljets" 46% 



T+jets 15% 



H+jets15% 




-c+t 



+jets 15% 



■dileptons" "lepton+jets" 

Fig. 3. Top quark pair branching fractions [42j , 
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beams were accelerated to 900 GeV. This phase is com- 
monly called Run I. After an upgrade the Run II started 
in 2001. The upgrade enables the Tevatron to accelerate 
the beams to a final energy of 980 GeV to yield a centre- 
of-mass energy of 1.96 TeV. Also the peak luminosity was 
gradually improved and now reaches a factor of approxi- 
mately 10 over the Run I performance [43ll44) . 

In Run I an integrated luminosity of about 160 pb^ 
was delivered to both experiments. This was sufficient 
to discover the top quark in proton-antiproton collisions 
with a centre-of-mass energy of 1.8 TeV [TT1|12| . With the 
increased centre-of-mass energy a 40% increased cross- 
section was achieved for top quark pair production. At 
the time of writing an integrated luminosity of more than 
6 fb~ was delivered to both experiments. Preliminary re- 
sults using up to 3.6 fb~ have been made public by the 
experiments. 

It is currently foreseen to continue running the Teva- 
tron at least until the end of 2010 and an extension into 
2011 is being discussed. The total integrated luminosity 
is expected to increase by about 2.5 fb~^ for each year of 
running. Thus until the end of the Tevatron program the 
total luminosity will more than double compared to what 
has been analysed so far [45i . 



1.2.2 The CDF Detector 

The Colhder Detector at Fermilab (CDF) |461I47[ H5] is one 
of two detectors recording collisions at the Tevatron. The 
vertexing and tracking components, the calorimetry and 
muon detection systems as used in the Run II measure- 
ments shall be described shortly in turn. 

The tracking system of CDF is placed in a 1.4 T 
solenoidal magnetic field. The heart of the CDF track- 
ing for Run II consists of three separate silicon detectors. 
The innermost (LOO) consists of a single layer of single 
sided sensors at a radius of about 1.5 cm. The upgraded 
silicon vertex detector (SVX II) consists of five cylindri- 
cal double-sided layers that along the beam axis reach to 
±45 cm from the centre of the detector at radii between 
2.5 cm and 10.6 cm [l^l50j . Finally the intermediate silicon 
layer (ISL) consists of double-sided sensors in one central 
layer at 23 cm and in two layers at 20 cm and 29 cm for the 
forward regions |51j . The tracking system is completed by 
wire drift chambers with an outer radius of 1.32 m. Their 
length of 3.2 m allows to cover the central region of jr;! < 1. 
The tracking systems provides full coverage and thus pre- 
cise tracking for I77I < 2.0 to measure the momentum of 
charged particles, to find the primary vertex of the col- 
lision as well as to find possible secondary vertices from 
long lived particles like b quark hadrons. 

Outside the tracking system a time of flight system 
(TOF) based on scintillator bars is positioned. This al- 
lows particle identification and is used in identifying b 
hadrons [52] . The solenoid that provides the magnetic field 
for the tracking system is placed between the TOF and the 
calorimetry. 

The central electromagnetic calorimeter is composed 
of alternating layers of lead and scintillator. It covers the 




Fig. 4. Schema of the CDF detector at Run II [49] . 



pseudo-rapidity range of jryj < 3.5. To detect more for- 
wardly produced particles layers of lead and proportional 
chambers are installed at pseudo-rapidities I77I < 4.2. The 
hadronic calorimeter uses lead as absorber material. The 
central region and the end cap wall use scintillator as ac- 
tive material. The forward regions also use gas propor- 
tional chambers. In total the calorimeter covers pseudo- 
rapidities of I77I < 4.2. The calorimeter allows energy mea- 
surements and identification of electrons, photons, jets 
and missing transverse momentum. 

The calorimeter is surrounded by several systems of 
drift chambers to identify muons. The 'central muon sys- 
tem' is placed at a radius of 347 cm. Behind an additional 
3.3 interaction lengths of 60 cm steel of the return joke 
the central muon upgrade system is located. Both sys- 
tems cover pseudo-rapidities of \ri\ < 0.6. The third sys- 
tem, the central muon extension, extends this coverage at 
0.6 < I77I < 1.0. So called barrel muon chambers extend 
the coverage for 1.0 < |ry| < 1.5. Track segments in these 
components are used to identify muons. 

CDF uses a three-level trigger system. The level 1 sys- 
tem has a pipeline for 42 beam-crossings. After the level 

1 trigger the event rate is approximately 10 kHz. At level 

2 trigger processors analyse a substantial fraction of the 
event and further reduce the rate to 200 Hz. The 3rd level 
consists of a cluster of computers which perform an op- 
timised event reconstruction. With this information the 
event rate is reduced to about 40 Hz. The events selected 
by level 3 are stored permanently to tape. 



1.2.3 The D0 Detector 

The D0 (Dzero) detector [551IM] is the second of the two 
detectors recording collisions at the Tevatron. It has been 
significantly upgraded to adapt to the reduced bunch dis- 
tance and increased luminosity of the Tevatron Run II. 
The upgraded D0 detector shall be shortly described here. 
The tracking system at the centre of D0 has been fully 
replaced since Run I and now consists of a silicon micro- 
strip tracker (SMT) and a scintillating-fibre tracker within 
a 2 T solenoidal magnet. The SMT consists of a barrel with 
four layers of single and double sided silicon micro-strip 



Daniel Wicke: Properties of the Top Quark 



detectors with a total length of about 70 cm. The barrel is 
separated into six subsections along the beam-pipe. Each 
subsection is capped with a disk of silicon detectors (F- 
disks). Three additional F-disks are placed on each side 
further outside of the barrel. Larger disks (H-disks) are 
placed at distances of 100 cm and 121cm from the beam 
pipe. The SMT barrel provides excellent r-0-information, 
the disks provide r-z as well as r-cf) measurements. During 
a shutdown period in 2006 the D0 silicon system was ex- 
tended by adding an additional layer of sensors directly on 
top of the beam pipe. This Layer-0 significantly improves 
the ability of vertex reconstruction. D0 measurements are 
usually performed separately for data taken before and 
after the installation of Layer-0. The two run periods are 
commonly denoted as Runlla and Run lib. Outside the 
SMT the central fibre tracker (CFT) is placed. It consists 
of scintillating fibres mounted on eight concentric support 
cylinders at radial distances between 20 and 52 cm cover- 
ing 1 77 1 to about 1.7. At each distance one layer of fibres 
is oriented along the beam axis and a second is mounted 
with stereo angles of of ±3°. 

To the outside of the CFT the solenoidal magnet is 
placed that produces a 2 T magnetic field for the tracking 
components. It is followed by the calorimetry that con- 
sists of a pre-shower detector that is placed in front of the 
cryostat and the sampling calorimeter based primarily on 
uranium/liquid-argon inside the cryostat. The end-caps 
that cover the forward regions have a similar structure, 
i.e. a pre-shower detector and a calorimeter cryostat. The 
calorimeters consists of three regions. The innermost is the 
electromagnetic calorimeter which uses thin (3 — 4 mm) 
plates of depleted uranium as absorber material. It is fol- 
lowed by the fine hadronic calorimeter with 6 mm thick 
plates of uranium- niobium. The outermost subsection, the 
coarse hadronic calorimeters, uses thick (about 47 mm) 
absorber plates of copper (in the central) and stainless 
steel (in the end-caps). The active medium in all three re- 
gions is liquid argon. Additionally, inter-cryostat-detectors 
(ICD) are mounted between the central and the forward 
cryostat to improve on the incomplete coverage of the 
calorimeters at 0.8 < \ri\ < 1.4. 

The D0 muon systems outside the calorimeter consist 
of tree layers of drift chambers, one inside and two out- 
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Fig. 5. Schematic view of the Run II D0 detector [55] 



side a toroidal magnetic field. The central muon chambers 
cover I77I < 1.0 with proportional drift chambers and have 
a magnetic field of 1.9 T. The forward muon chambers use 
mini drift chambers and a toroidal field of 2.0 T. They 
extend the coverage to \j]\ « 2.0. 

Also D0 uses a trigger system with three stages. The 
first level consists of a set of hardware trigger elements 
that provide a trigger accept rate of 2 kHz. In the sec- 
ond level hardware engines and embedded microprocessors 
provide information to a global processor that considers 
individual detector information as well as correlations. It 
reduces the rate by a factor of about 2. The third level 
consists of a farm of computers that reduces the rate to 
50 Hz based on a limited event reconstruction. The ac- 
cepted events are stored to tape for offline analysis. 



1.3 Basic Event Selection 

The selection of events in general, in particular the selec- 
tion of top quark events, is based on the reconstruction 
of a number of different objects: The primary vertex of 
the collision, electrons, muons, jets and transverse miss- 
ing energy. In addition, to reduce the background, often 
methods to identify jets from 6-quarks are applied. In the 
following these objects shall be shortly described in turn. 

The reconstruction of the primary vertex is determined 
by assigning well measured tracks to a common origin in 
the interaction region. The primary vertex is constructed 
event- by- event and is used as reference for some of the 
following objects. 

Electrons are reconstructed by a combination of track- 
ing and calorimeter information. Quality cuts on the 
tracks typically include a -pT threshold of the order of 
10 GeV and the matching energy deposit in the calorime- 
ters should be well contained in the electromagnetic 
subsection. In the absence of bremsstrahlung the en- 
ergy deposit is expected to have a small radius in the 
77-0 plane. Sophisticated algorithms take into account 
bremsstrahlung photons. Discriminating observables in- 
clude the relative sizes of the energy deposit in the 
hadronic and the electromagnetic calorimeter. D0 uses 
a fixed ratio of the electromagnetic to the total energy 
/cm = Ecn\/E < 0.9 |SB] while CDF uses an energy de- 
pendent cut on £;om/^had < 0.055-^0.00045 GeV'^E [57]. 

Muons are identified by the presence of signals in the 
dedicated muon chambers that can be matched to a track 
found in the tracking system. For this purpose CDF and 
D0 extrapolate isolated tracks with standardised quality 
cuts through the calorimeter out to the muon chambers 
to find matching track segments. 

For both lepton objects different quality classes are 
defined, named "loose", "medium", "tight", etc, with in- 
creasingly stricter requirements on the isolation of the 
track and the calorimeter cluster. These different criteria 
can be used to obtain an acceptable signal to background 
ratio. The selection of "loose" leptons excluding those that 
also have a "tighter" identification is often used to define 
sideband samples to extract background estimates from 
data. The triggering and reconstruction efficiencies are 
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usually studied in Z ^ £i events using the tag and probe 
method. 

Jets are reconstructed from all calorimeter towers. 
They usually show a substantial contribution from the 
hadronic subsection and are usually broader than signals 
from electromagnetic particles. At the Tevatron experi- 
ments the "improved legacy" cone algorithms [55] with 
radii of 0.4 and 0.5 are used by CDF and D0, respec- 
tively. Quality cuts typically require that the energy of a 
jet is not contained to more than 90% in a single tower 
and deposits from electron and photon candidates are re- 
moved. The jet energy reconstructed from the calorimeter 
cells needs to be corrected for a number of effects. These 
corrections include imperfections of the calorimeter but 
also energy offsets due to contributions from the under- 
lying event, multiple hadron interactions and noise in the 
electronics. This correction, usually called the Jet Energy 
Scale (JES), is obtained from precisely measured electro- 
magnetic objects by assuming momentum conservation in 
the transverse plane for 7-l-jets events. 

Jets stemming from 6-quarks can be identified due to 
the long lifetime of about 1.5 ps of the B hadrons in such 
jets. At the typical energies in top quark events of 50 to 
100 GeV the mean decay length is of the order of 5 mm. 
This fact is exploited by computing the impact parame- 
ter for tracks or by explicitly reconstructing a secondary 
vertex from the tracks that is displaced from the primary 
vertex. To identify fo-jets CDF uses only tracks that fall 
within the cone radius of the considered jet. A secondary 
vertex is reconstructed in two passes with different track 
requirements [59] . The 2d decay length is computed as the 
distance of primary to the secondary vertex. The signifi- 
cance, i.e. the decay length over its uncertainty, is required 
to be larger than 3. If the direction from the primary to 
the secondary vertex is opposite to the direction of the 
jet, the tag is called a negative tag. These negative tags 
are useful to determine the purity of mis-tag rates. In D0 
first track jets are built with a cone algorithm indepen- 
dent of the calorimeter jets described above. Secondary 
vertices are reconstructed with tracks from a given track 
jet. Calorimeter jets are considered as identified 6-jets, if 
an identified track jet falls within a radius of AR < 0.5. 
The most recent D0 tagging algorithm uses the impact 
parameters of tracks matched to a given jet and infor- 
mation on the secondary vertex mass, the significance of 
displacement, and the number of participating tracks for 
any reconstructed secondary vertex within the cone of the 
given jet. The information is combined in a neural network 
to obtain the output variable, NN^, which tends towards 
one for 6-jets and towards zero for light jets [50] . 

The presence of a neutrino is inferred from the mo- 
mentum imbalance in the transverse plane, that occurs 
because neutrinos are invisible to the detectors. If all ob- 
jects of an event were measured the sum of the transverse 
momenta should vanish. Thus the sum of the transverse 
momenta of all neutrinos can be deduced from the miss- 
ing transverse momentum needed to assure momentum 
conservation. It is derived from the calorimetric measure- 
ments and their direction with respect to the interaction 



region as the negative vector sum of the transverse ener- 
gies and thus usually named Missing Transverse Energy, 
^T- The sum is taken over all calorimeter cells that remain 
after noise suppression. In general corrections for identi- 
fied objects with known energies like electrons, muons or 
jets are applied. 

Triggering and preselection of single top quark events 
and of top quark pair events in the semileptonic and 
dileptonic channels are based on reconstructed lepton, jet 
and ^T objects described above. Typically the leptons 
and the missing transverse energy are required to have 
Pt > 20 GeV. Signal selection in addition requires the 
presence of jets also with a typical momentum require- 
ment Pt > 20 GeV. The number of jets required of course 
depends on the channel under consideration. Some anal- 
yses here improve the signal to background ratio by re- 
quiring one or more identified 6-jets. Others construct a 
likelihood for an event being top quark like based on topo- 
logical quantities to consciously avoid 6-jet identification. 
Additional cuts may be introduced to enhance the sig- 
nal to background ratio or to improve the data to Monte 
Carlo agreement. In the semileptonic selection D0 e.g. re- 
cently requires the leading jet to fulfil pr > 40 GeV and 
avoids events in which missing transverse momentum and 
the selected lepton are aligned. 

Events of the all-hadronic channel have to be selected 
by requiring multijet final states. Due to the overwhelming 
background from multijet events in these analyses usually 
stronger cuts are applied on the transverse momentum of 
the jets. 



2 Top Quark Mass Measurements 

The mass is the only property of the top quark that is 
not fixed within the Standard Model of particle physics. 
The Yukawa coupling responsible for the coupling of the 
top quark to the Higgs boson and thus for the mass of 
the top quark is a free parameter of the Standard Model. 
This already illustrates the importance of measuring the 
top quark mass. 

Moreover Standard Model predictions of electroweak 
precision observables depend on the value of the top quark 
mass via radiative corrections. By correlating the W bo- 
son mass with the top quark mass the mass range for the 
yet undiscovered Higgs boson can be constrained. Mea- 
surements of the top quark mass are thus an important 
preparation for discovering the Higgs boson and will serve 
as a consistency check of the Standard Model after its 
discovery. 

In the following theoretical aspects of fermion mass 
measurements are discussed, before the experimental 
methods used in the various decay channels are explained. 
Then issues of modelling non-perturbative effects in pp 
collisions and their influence on the existing measurements 
are discussed and the combination of the various measure- 
ments to a final result is reviewed. The conclusions contain 
a critical comparison of current results with expectations 
and future prospects. 
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2.1 Theoretical Aspects 

For free particles the physical mass is usually taken to 
correspond to the pole of their propagator, i.e. their value 
of the four-momentum squared, p^ — E"^ —p'^. Because of 
confinement quarks cannot exist as free particles and this 
definition becomes ambiguous. 

The definition of the pole mass is still possible on an 
order by order basis in perturbation theory, but is con- 
sidered to be intrinsically ambiguous on the order of the 
confinement scale 0{Aqcu) [M]. Mass definitions can also 
be obtained following other renormalisation schemes such 
as the MS-mass in the MS-scheme. Other definitions have 
been suggested by [62ll63] . 

The relation between the various mass definitions can 
usually be computed in perturbation theory. For the top 
quark the numerical values of the different definitions may 
differ significantly. In NNLO for example the MS-mass of 
the top quark is about 10 GeV smaller than the pole mass. 

This large difference makes it important to understand 
which definition and order of perturbation theory of the 
top quark mass is measured by the experiments. As we 
shall see below, all direct methods to determine the top 
quark mass use Monte Carlo simulation to either extract 
the mass or calibrate the procedure. Thus the question 
really is, which definition of the top quark mass is used as 
a parameter in these Monte Carlo generators. 

Unfortunately, it is not well understood which field 
theoretic mass definition the currently used generators 
correspond to. Clearly, as the hard matrix elements of 
(most) generators is implemented in leading order, they 
correspondingly use the leading order mass. It is usually 
argued that this corresponds to the pole mass, because in 
the pole mass definition any shifts of the position are to 
be absorbed in the mass definition. Monte Carlo gener- 
ators do not absorb corrections from the parton shower 
or the hadronisation into the mass definition. And it is 
not clear in how far the parton shower and the modelling 
of hadronisation alter the mass definition, nor which ap- 
proximation of QCD this mass parameter corresponds to. 
Partial answers have been given in |63] . but at this point 
a conceptual uncertainty remains that is considered to be 
of the order of 1 GeV. 

Comparisons of the top quark mass from direct mea- 
surements with electroweak precision data to determine 
e.g. the Higgs boson mass currently assume the measured 
top quark mass values corresponds to the pole mass defi- 
nition. They thus have to be interpreted with care. 

Experimentally precise measurements of the top quark 
mass are nevertheless useful. On the one hand they are 
needed to set the mass parameter for simulating top quark 
events, which will be important backgrounds for some 
LHC searches. On the other hand the consistency between 
the various experimental methods and the top quark de- 
cay channels gives confidence in these methods. Finally, it 
seems feasible that a theoretically more precise specifica- 
tion of the top quark mass definition used in the Monte 
Carlo generators can be derived in the future j64]. Any 
bias from the interpretation as a pole mass may then be 
corrected for. 



2.2 Lepton plus Jets Channel 

The lepton plus jets channel is considered to be the golden 
channel in top quark mass measurements. Due to a lepton 
and a neutrino in the final state it has a good signal to 
background ratio. In addition one of the top quark quarks 
is fully measured in the detector. 

Several method have been applied by CDF and D0 to 
measure the top quark mass in this channel. The methods 
exploit the kinematics of the events to different levels and 
make different assumptions about details of the produc- 
tion mechanisms. 



2.2.1 Template Method 

The basic template method was already applied to lepton 
plus jets events in the papers describing the first evidence 
and discovery of the top quark. 

In this method a top quark mass is reconstructed 
for the each of the selected events using the momenta 
measured for lepton and jets and the transverse missing 
energy. The distribution of the reconstructed masses is 
then compared to template distributions from simulation. 
These templates are constructed from signal Monte Carlo 
with varying top quark mass values and contain the ex- 
pected amount of background events. The method thus 
relies on good Monte Carlo modelling of signal and back- 
ground. 



CDF Run II Template Measurement CDF has applied 
an improved template method in up to 5.6 fb^ of data 
combining the lepton plus jets channel and the dilepton 
channel |551IB51I57[I68| . For clarity the two analysis parts 
will be described separately. The dilepton portion of this 
analysis can be found in Section [2.3.1l 

Lepton plus jets events are selected requiring a single 
high pt lepton, large missing transverse momentum and at 
least four jets. Events are separated by the number of jets 
identified as 6-jets based on the transverse decay length of 
track inside the jet [BS]. In case of only a single identified 
jet only events with exactly four jets are considered. For 
events with more than one identified 5-jet more than four 
jets are allowed. 

In each event a top quark mass ml°'^° is fitted using 
a constrained fit. Besides the top quark mass, the mo- 
menta of the top quark decay products (the quarks and 
leptons) are fitted to the observed transverse jet and lep- 
ton momenta, p^^, and the unclustered energy, If^^. The 
unclustered energy is the calorimetric energy in the trans- 
verse direction not associated with any reconstructed ob- 
ject. Jet momenta are corrected to parton level with CDFs 
common jet energy scale correction. The invariant masses 
of the W boson decay products on both sides, Mgg and 
Mil,, are constrained to be consistent with the nominal W 
boson mass within the W boson width. The reconstructed 
top quark mass ml'^'^° is required to be consistent within 
the top quark width with the invariant mass of the top 
quark decay products on both sides, Mbgg and Mte,^. The 
fit x^ thus is written as 



10 



Daniel Wicke: Properties of the Top Quark 




110 



iiO i4(i 160 idO idO liO 246 
mjeco (GeV/c^) 



0.001 


CD 
-100 

3~ 


0.0008 


90 


0.0006 


80 


0.0004 


70 


0.0002 


60 




50 




iiO i40 i6d idO iOd 226 240 

mreco (GeV/C^) 



I 



Fig. 6. Probability densities for signal (left) and background (right) for nominal top quark mass of 170 GeV and no jet energy 
scale shift as obtained for single b-tag events |66) . 
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(11) 



The fitted unclustered energy, U^ , is related to the trans- 
verse momentum of the neutrino used in the computation 
of Mil, and Mmv This y^ is computed for all possible as- 
sociations of quarks to four jets, allowing 6-jets only to be 
matched with b quarks. Only the top quark mass value of 
the association with the best fit x^ is kept. If this best 
X^ > 9.0 the event is rejected. 

To constrain the jet energy scale simultaneously with 
the top quark mass, the dijet mass, rrijj, is measured from 
the non-6-tagged jets among the four leading jets without 
applying the above kinematic fit. This is only unique for 
events with two identified jets. In other events the two jets 
that yield the dijet mass closest to the Vl^-boson mass are 
chosen. 

Thus for each event one top quark mass value, m\'^^° , 
and one dijet mass, rrijj, is entering the following analysis. 

Monte Carlo simulations are used to determine the 
expected behaviour of signal events as well as the back- 
ground contributions. The dominant M^-|-jets background 
is simulated with Alpgen-|-Pythia [TDlfTT] with a nor- 
malisation derived from data. The multijet background 
is modelled with samples containing non-isolated leptons. 
Smaller background from single-top and diboson events 
are taken from Monte Carlo with normalisation to theoret- 
ical cross-sections. All backgrounds are assumed to have 
no dependence on the nominal top quark mass, but are 
allowed to vary with the jet energy scale. Signal samples 
for various nominal top quark mass values, m^, are gen- 
erated using Pythia. All simulations are passed through 
the full CDF detector simulation and reconstruction. 

These simulations are used to generate probability 
density functions, P^^'s, to find a signal event with mea- 



sured values 



and 



'J J 



given a nominal top quark 



mass, mt and jet energy scale shift, Z\jes- Similarly a 



probability density for background events, P°^s^ jg com- 
puted as function of the jet energy scale shift, only. CDF 
uses a Kernel Density Approach, where each simulated 
event contributes not only with its reconstructed values 



and 



'jj' 



but also in the neighbourhood around 



these. The size of this neighbourhood is controlled by 
smoothing parameters. The parameters are chosen dy- 
namically: small near the maximum of a given distribution 
where the statistics require less smoothing and larger in 
the tails. Examples of the resulting density functions are 
shown in Fig. [51 

With these probability densities the likelihood 
£(TOt,Z\jEs) is constructed. The number of signal and 
background events in the single 6-tag and the double h- 
tag sample are used as additional parameters with Gaus- 
sian constraints on the background estimations obtained 
above. Also the jet energy shift, ZijESi is constrained by a 
Gaussian term to be consistent with zero within its nom- 
inal uncertainty, tJc- 



£(mt,Z\jEs) 



exp 



A 



JES 



2a, 



X £1(7714, Z\jEs) X £2(TOt,^JEs) 



£i(mt,Z\jEs) 



exp 






n 



Si + hi 



(12) 



Here i = 1 , 2 indicates the subsamples with one or more 
than one identified 6-jets, respectively. Si and hi are the 
number of signal and background events and h° the back- 
ground expectations in the samples. 

The described likelihood formulae can only be evalu- 
ated at the discrete values of mt and ZijES for which sim- 
ulations were run. To obtain the likelihood for arbitrary 
values of mt and /4jes a quadratic interpolation is used. 
The central result is obtained by maximising the likelihood 
and its uncertainty is quoted as the (largest) mass shift 
corresponding to a likelihood change of Z\log£ — 0.5. 

Systematic uncertainties are evaluated by modifying 
several aspects of the analysis described. The dominat- 
ing uncertainty in top quark mass measurements is the 
jet energy scale. Through the simultaneous fit its contri- 
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Fig. 7. Likelihood contours for the combined lepton plus jet 
and dilepton fit of 5.6 fb"^ CDF data [68]. 



bution is part of the statistical uncertainty in this mea- 
surement. Residual effects remain through uncertainties in 
the pt and 77 dependence. The uncertainty of modelling 
the top quark pair signal events is evaluated by comparing 
pseudo experiments generated with Pythia and Herwig. 
These two systematic uncertainties contribute with about 
0.7 GeV to the uncertainty and are the two largest single 
contributions to the systematic uncertainty. Additional 
contributions in order of decreasing importance include 
the uncertainties on modelling colour reconnection [7^1751 
[74] . the background shape, the parton density functions. 
With 5.6 fb~^ of data using the lepton plus jets channel 
only CDF determines a top quark mass of |68j 



mt - 172.2 ± 1.2stat ± 0.9,yst GeV 



(13) 



For the jet energy scale shift the fit yields Z\jes — (0.10 ± 
0.26)crc, very consistent with the external jet energy scale 
but with a significantly reduced uncertainty. Combining 
the results with an analogous measurement in the dilep- 
ton channel which uses a consistent jet energy scale shift 
(c.f. Section [^Xn yields 



mt = 172.1 ± l.lstat ± 0.9,yst GeV 



(14) 



The likelihood contour of the two dimensional measure- 
ment of mt and Z\jes is shown in Fig. [71 



D0 Run II Template Measurements D0 has produced 
preliminary results for top quark mass measurements us- 
ing the template method [751(76] . The method differs from 
the CDF method as the templates are not smoothened and 
in the use of a binned likelihood. As these results were not 
updated nor published in the last years they shall not be 
described in more detail here. 



2.2.2 Ideogram Method 

The ideogram method was transfered to top quark mass 
measurements from a method used in the DELPHI W 
boson mass measurement |77ll78j . 



Also this method relies on an event by event recon- 
struction of top quark mass values. But in contrast to the 
template method the signal likelihood is constructed event 
by event from the theoretically expected Breit-Wigner dis- 
tribution smeared with the experimental resolution of each 
individual event. Events with a configuration that allows 
a more precise reconstruction thus contribute more than 
events with a configuration that is difficult to reconstruct. 



D0 D0 has published an analysis of 425 pb^ of data 
based on semileptonic top quark pair events [75]. The 
event selection requires an isolated lepton, missing trans- 
verse energy and four or more jets. To identify h quark 
jets the decay length significance of a secondary vertex re- 
construction is used ^80^. No cut is placed on the number 
of identified 6-jets. 

The top quark mass in each event is reconstructed us- 
ing a constrained fit that determines the momenta of the 
top quark pair decay products {(.vhhqq') from the mea- 
sured momenta of the charged lepton, the four leading 
jets and the missing transverse momentum and their un- 
certainties. Constraints are placed on the invariant mass 
of the two light quarks and the two leptons, respectively, 
which are required to be consistent with the W boson 
mass. The reconstructed top and anti-top quark masses 
are required to be equal. The 12 possible assignments of 
jets to quarks that yield different constraints are consid- 
ered. In addition two possible solutions for the neutrino z 
momentum are considered, which results in 24 top quark 
mass values per event. As the common jet energy scale 
of D0 corresponds to the particle level, i.e. what would 
be visible in an ideal detector, the fit uses jet-parton map- 
ping functions determined from Pythia simulation. These 
mappings contain an overall scale factor, /jeSj for the jet 
energy scale that may modify the default jet energy scale. 

The described fit is repeated for various values of the 
jet energy scale factor, so that the resulting mass values 
and the fit x^ &re functions of /jes. Through the con- 
straints to the nominal W boson mass, the fit x^ is ex- 
pected to be the smallest for the correct jet energy scale 
(and top quark mass value). Events with no jet-parton 
combination reaching x^ < 10 at the central jet energy 
scale are rejected at this point. 

To determine the expected performance of the selec- 
tion and the corresponding background contamination, 
top quark pair signal for various nominal top quark mass 
values and M^-fjets background events are generated us- 
ing Alpgen-(-Pythia. The events are passed through the 
full detector simulation and reconstruction. Multijet back- 
ground is modelled with side band data obtained from 
inverted lepton quality cuts. The sample composition is 
determined from a likelihood discriminant that combines 
topological observables, with tracking based jet shape and 
h tagging information, c.f. Fig. [5] The observables are se- 
lected to have low correlation to the top quark mass. 

An event by event likelihood is now constructed to 
observe the top quark masses reconstructed in the kine- 
matic fit, the discriminant and the number of identified 
&-jets given the nominal values of mt and /jes^ 
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Cevt{x;mt,fjES,ftop) = /top ^sig(a;;TOt,/jEs) 

+ {1 - ftop)Pbkg{x; fjEs) ■ (15) 

Here x represents the measured observables for the event 
under consideration and /top is the signal fraction of the 
corresponding sample. The signal and background proba- 
bilities, Psig and Phkg, are proportional to the probabilities 
to find signal or background at the observed discriminant 
value. Because of the selection of observables in the dis- 
criminant, this factor does not depend on the top quark 
mass and the jet energy scale and is factorised as: 

Psig(a;;mt,/jEs) = Psig{D) Psig{x;nit, f:iEs) 

-Pbkg(a;; /jes) = ^bkgC-D) Pbkg(a;; /jes) • (16) 

The top quark mass and jet energy scale dependent por- 
tion of the signal probability, -Psig, is computed as the sum 
over all 24 jet-parton assignments and neutrino solutions. 
Their relative weight, Wi, is computed from the x^ prob- 
ability of the kinematic fit and in presence of identified 
6-jets includes the probability for the quarks associated to 
the fo-jets to produce a tagged jet. Because the result of 
the kinematic fit depends on the jet energy scale factor 
the weights depend on /jes- 

24 

Psig(x; nit, /jes) = X! '^i(/jEs)5'(mi, /jes) 

i=l 
24 

-Pbkg(a;; /,7Es) = ^ Wi{fjEs)B{'mi) (17) 



The background probabilities B{mi) are obtained from 
the simulated H^-|-jets events. The signal probabilities 
'S'(mi,/jEs) are computed as the convolution of a Breit- 
Wigner, BW, that describes the theoretical distribution of 
mass values and a Gaussian, G, that represents the detec- 
tor resolution. As this ansatz is valid only for the correct 
jet-parton assignment and because the wrong pairings do 
contain information about the top quark mass, a second 
contribution to describe wrong pairings is added: 



S{mi,f jes) = fc G(mi,m',cri)BW(m',mt) dm' 

+{1- fc)SwYongimi,mt;nti,g) . (18) 

with fc being the fraction of events in which the weight is 
assigned to the correct jet-parton pairing. For the width 
of the Gaussian the uncertainty, cr^ , of the mass determi- 
nation obtained in the constrained kinematic fit is used. 
Thus events with a well-determined mass from the kine- 
matic fit contribute more than events with a less precise 
fit result. The expected mass value distribution in wrong 
pairings, S'wrong, also contains information about the true 
top quark mass. It is determined from simulation as func- 
tion of the number of identified 6-jets. Also the background 
function, B{mi), is taken from simulation. 

Due to the presence of wrong jet-parton assignments 
and background events the described likelihood does not 
yield unbiased results for the jet energy scale factor, /jes- 
The likelihood is thus corrected with an /jES-dependent 
but mass-independent correction factor. 

The likelihood for the complete sample of observed 
events is then simply the product of the jet energy scale 
corrected likelihood, C^"", for all individual events: 



^("It, /jes, /top) = Jl 



Cl°"{x] mt, /jes, /top) 



(19) 



This likelihood is maximised simultaneously with respect 
to the top quark mass, mt, the jet energy scale factor, /jes 
and the signal fraction /top- 

Before this procedure is applied to data its perfor- 
mance is determined on large numbers of pseudo exper- 
iments with varying nominal parameter values- The bias 
is determined as the mean difference between the nominal 
value and the measured result. The correctness of the fit 
uncertainty is checked from the distribution of pull values, 
i-e- the distribution of differences between the measured 
values and their means normalised to the individual fit 
error. The bias and the width of the pull distributions for 
the fitted top quark mass and jet energy scale factor, /jes, 
are determined as function of the nominal top quark mass 
and jet energy scale factor simultaneously. Linear correc- 
tions are applied to correct for the obtained biases and to 
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Fig. 9. Contours for tlie 2-dimensional likelihood of determined for e+jets (left), /^-fjets and the combined dataset in 
425 pb~^ [79) . The contours correspond to log likelihood differences of 0.5, 2.0 and 4.5. 



correct the uncertainty for any deviation of the pull width 
from the ideal value of one. 

Systematic uncertainties for this measurement are de- 
termined from pseudo experiments with events shifted ac- 
cording to the systematic variation under study. Due to 
the two dimensional fit the uncertainty of the overall jet 
energy scale is contained in the uncertainty obtained from 
the likelihood fit. Residual discrepancies between the data 
and Monte Carlo energy scales still affect the result and 
give the largest contributions to the systematic uncer- 
tainties. D0 evaluates uncertainties due to b quark frag- 
mentation modelling [81; and the calorimeter response to 
6-jets to ±1.30 and ±1.15 GeV, respectively. The uncer- 
tainty from the pT of the jet energy scale yields ±0.45 GeV. 
Another large contribution of ±0.73 GeV comes from the 
signal modelling which is estimated by varying the frac- 
tion of high energy gluons produced in addition to a top 
quark pair. Further uncertainties include the influence of 
uncertainties on the trigger efficiencies, the background 
modelling and the calibration. 

The ideogram method with in-situ jet energy scale cal- 
ibration applied in 425 pb^ of data : t9j yields a top quark 
mass of 



mt 



173.7 ± 4.4s 



i-JES-2'osyst GeV 



(20) 



the jet energy scale factor is determined to be /jes — 
0.989 ± 0.029stat consistent with the nominal value of 
/jes = 1-0, which corresponds to the calibration obtained 
in jet-photon events. The fit probability contour lines for 
the individual lepton channels and the combined results 
are shown in Fig. [SI 



2.2.3 Matrix Element Method 

The Matrix Element method in the lepton plus jets chan- 
nel was pioneered by D0 in Run I [52] and is based on 
ideas described in [55115^ . In this method the likelihood is 
constructed according to the expected distribution from 
the mass dependent top quark pair production matrix el- 
ement. With respect to the previously described methods 
this includes additional information from top quark mass 
dependent kinematics into the measurement. 



D0 Results for the top quark mass using the Matrix Ele- 
ment method have been regularly updated with D0 Run II 
data [85,gBl [571[55] with the latest result including 3.6 fb"^ 
of data. The analyses select events from semileptonic top 
quark pair decays by requiring a single isolated lepton, 
missing transverse momentum and exactly four jets. In the 
recent analyses at least one of the jets needs to be iden- 
tified using D0's neural network b jet identification |60j . 
Vetoes are applied on additional leptons and on events in 
which the lepton and the missing transverse momentum 
are close in the azimuthal direction. 

The top quark mass and an overall jet energy scale 
factor are determined simultaneously from an event by 
event probability that the observed event may occur given 
the assumed values of the top quark mass, mt, the jet 
energy scale factor, /jes and the signal fraction /top: 



Pcvt(a;; wt, /jes, /top) = 
M^) ■ (/top^sig {x;mt, /jes) 



(21) 
(l-/top)Pbkg(2^;/jEs))- 



Here a; represents the measured momenta of the lepton, 
the jets and the missing transverse momentum. A{x) is the 
probability for the configuration x to be accepted in the 
analysis. The probabilities, Psig and -Pbkg, for signal and 
background events to be observed in the configuration x 
given the set of parameters, mt, /jes and /top, are com- 
puted from matrix elements of the dominating processes. 
The signal probability is computed from the ma- 
trix element for top quark pair production and decay 
through quark anti-quark annihilation AAti{y\'mt), con- 
voluted with the transfer function, W{x, y; /jes), that de- 
scribes the probability to observe a parton configuration, 
y, as the measured quantities, x. 



Psigix;mt,fjEs) 



0-obs(TOt,/jEs) 

Y^ dqidq2d$e f{qi)f{q2) 
ivours 

{27r)^\Mtiiy;mtr 



(22) 



qiq2s 



W{x, y;f ]Es) 



The sum is over the possible flavours of the incoming 
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quarks and the integral over their momentum fractions, q\ 
and qi , as well as the 6-body phase space for the outgoing 
particles, ^g- /('Zi) are the parton densities for the incom- 
ing quarks and s is the centre-of-mass energy squared. The 
normalisation factor (Tobsl?^*, /jes) is the cross-section ob- 
servable with the selection used, i.e. it includes effects of 
efficiencies and geometric acceptance. 

For the background probability a corresponding for- 
mula is used with the top quark pair matrix element re- 
placed by the M^-|-jets matrix element, which of course is 
independent of mt- The contribution from multijet events 
is considered to have a similar shape and is not included 
separately. 

Both, Psig and Pbkg contain the same transfer func- 
tion, VF(x, y; /jes)- It is derived from full simulation for 
individual jets and leptons. Only changes of the size of 
the momenta but not of the directions are considered. Be- 
cause it is not known which jet stems from which parton, 
a weighted sum over all 24 possible assignments is made. 
The weight, Wi, reflects the probability of the event's h- 
tags to be consistent with the jet-parton assignment under 
consideration: 






Wix, y; /jes) = WfXxt, yi] /jes) 2^w. 



J|W^j(a;j,2/ij;/jEs)- 

(23) 
The values xg, and yi are the measured and the assumed 
momenta of the lepton, Xj is the measured momentum of 
the jth jet and yi^j the momentum of the matrix element 
parton associated to the jth jet in the jet parton asso- 
ciation number i. We and Wj are the transfer functions 
for leptons and jets, respectively, which are zero when the 
directions do not coincide. 

The likelihood to observe the measured data is com- 
puted as the product of the individual event probabilities, 

-Pevt^ 

'C(r7it,/jEs,/top) = Y[ Pcvt {x ; mt, f^ES, f top) ■ (24) 

events 

For each assumed pair of the nominal top quark mass 
and the jet energy scale factor, rrit and /jeSi the likeli- 
hood is maximised with respect to the top quark fraction. 



/top- In 85, '861187] the top quark mass and jet energy scale 
factor are then determined by maximising the two dimen- 
sional likelihood £(TOt, /jes, /top*(("^t, /jes))- In [SH] the 
jet energy scale factor is constrained with a Gaussian prob- 
ability distribution to its nominal value and its uncertainty 
as obtained from photon plus jet and dijet events. 

Before the method is applied to data, its performance 
is calibrated in pseudo experiments. Random events are 
drawn from a large pool of simulated top quark pair sig- 
nal and M^-|-jets background events with proper fluctua- 
tions of the signal and background contribution such that 
the total number of events corresponds to the number of 
events observed in data. The simulated events were gener- 
ated with Alpgen-I-Pythia and passed through the fuU 
D0 simulation and reconstruction. This procedure is re- 
peated 1000 times for several fixed nominal values of mt, 
/jes and /top- Thus for each of the nominal values the 
signal fraction, the top quark mass and the jet energy fac- 
tor can be measured with the described procedure 1000 
times. The mean result for each set of pseudo experiments 
with fixed nominal m,t , /jes and /top are compared to the 
nominal values in the calibration curves in Fig. [TU] The 
final result is corrected for any deviation of these calibra- 
tion curves from the ideal diagonal. Also the pull width 
is computed and the statistical uncertainty corrected ac- 
cordingly. 

The leading source of systematic uncertainty con- 
tributing ±0.81 GeV (Runllb) stems from the modelling 
of differences in the detector response between light 
quark and b quark jets. The next important contribu- 
tion arises from uncertainties modelling hadronisation and 
underlying event, ft is estimated from the difference be- 
tween Pythia and Herwig and contributes with nearly 
±0.6 GeV. Also the sample dependence of jet energy scale 
corrections in simulation contributes with this size. For 
the first time in [55j this measurement includes an es- 
timate of the uncertainty due to colour reconnection ef- 
fects [7211731174] ■ which contributes 0.4 GeV to the uncer- 
tainty. The squared sum of the individual contributions 
yields a total uncertainty of ±1.4 GeV (Run lib). D0 has 
applied this analysis to their 3.6 fb~^ dataset separated by 
run periods. The 1.0 fb~ Runlla result and the 2.6 fb~ 
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energy scale factor obtained with the Matrix Element method 
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Run lib result, shown in Fig. Illl are combined with the 
BLUE method [89, following the error categories used by 
the Tevatron Electroweak Working Group [9D] (see also 
subsection 12 . 71 below) and yield 



mt = 173.7 ± O.Sstat ± 1.6syst GeV 



(25) 



where in this combined result the uncertainty due to the 
overall jet energy scale factor is contained in the system- 
atic uncertainty. This measurement is currently D0's most 
precise top quark mass result. 



CDF The CDF collaboration is using the concept to mea- 
sure the top quark mass from a likelihood based on the 
production matrix element in the lepton plus channel in 
several variations pTllMligHllMligSl IM] . 

The recent analyses all base on an event selection that 
requires a single isolated lepton, missing transverse mo- 
mentum and exactly four jets, at least one of which is re- 
quired to be identified as b jet. Then the top quark mass 
and an overall jet energy scale factor are determined si- 
multaneously from an event by event probability that the 
observed event may occur given the true values of the top 
quark mass, mt, the jet energy scale factor, /jes and the 
signal fraction /top- The various CDF analyses differ in 
the construction of the likelihood. 

The CDF Matrix Element method (MEM) [9T1I92] fol- 
lows closely the procedure outlined in Section [2.2.31 The 
matrix element used to describe signal is that of the qq — >■ 
if process with its decay. For background the iy-|-4jets ma- 
trix element is employed. Finally, the transfer functions of 
Eq. (P5)) use only parton-jet assignments consistent with 
6-tagging information and assumes the lepton measure- 
ment to be exact. For the background probability, Pbkg, 
the dependence on the jet energy scale is taken from a pa- 
rameterisation of the average likelihood response rather 
than an explicit change of the transfer functions. 



The Dynamical Likelihood method (DLM) 
differs from the Matrix Element method in that 
it bases its likelihood only on the signal contribu- 
tion, Psig(a;;r7T.t, Z\jEs)- In this signal term the squared 
matrix element is factorised into a production ma- 
trix element, the (anti-)top quark propagators and 
their decay matrix elements: \Ai{aia2 -^ tt -^ 
eMqq')\^ = \Ma,a,^ti?VtVi\Vt\''\Vi\\ which removes 
spin-correlations. However, in contrast to previously de- 
scribed measurements it includes gluon diagrams, i.e. 0102 
can be qq or gg. For the treatment of background contribu- 
tions the method fully relies on results obtained in ensem- 
bles of pseudo-datasets similar to the calibration step of 
the other methods. This correction is computed depend- 
ing on the jet energy scale correction and the background 
fraction. 

The Matrix Element Method with Quasi-MC- 
Integration (MTM) [57119811991 uses an again more 
complete matrix element to construct the signal likeli- 
hood. The applied matrix element J 100' includes the qq 
and the gg production channels with full spin-correlations. 
The method treats the background by subtracting the 
expected contribution of the logarithm of the likelihood. 
Equation ((2T|l is thus rewritten as 

logPcvt(a;;TOt,Z\jEs) 

= log Psig (a; ;mt,Z\ jes) - /bkg(g) logPbkg(mt, Z\jes) (26) 

with log Pbkg being the average log Pgig obtained in back- 
ground events. The expected background fraction, /bkg, 
is computed from simulation and applied event by event 
as function of the output 5 of a neural network discrim- 
inant, c.f Fig. 1121 Moreover the analysis removes events 
which are likely to be mismeasured, e.g. due to extra jets, 
misidentification etc. For this events are required to have 
logPovt > 10 for at least some range of rut- 

In all methods the overall likelihood is computed from 
the various event likelihoods following Eq. (|24l) . It is max- 
imised to find the optimal top quark mass and jet energy 
scale parameters. For the MEM also the background frac- 
tion is fitted. The performance is then checked by applying 
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Fig. 13. Results of CDF top quark mass measurements with in-situ jet energy calibration obtained in three variations of the 
Matrix Element method. Left: Matrix Element method a la D0 (MEM) ^, Middle: Dynamical Likelihood method (DLM) [94] . 
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the top quark mass measurement on ensembles of pseudo- 
experiments with various nominal values of the top quark 
mass, nit, and the jet energy scale shifts, /\jes- The re- 
sults are corrected for any observed shifts and the uncer- 
tainties scaled to fit the observed spread of results in the 
ensemble tests. The observed shifts vary between 1.4 GeV 
for the Dynamical Likelihood method [23] ; which relies on 
this calibration to describe the background effects, and 
0.09 GeV for the Matrix Element method [92^, which has 
the most complete background term. 

As all methods apply an in-situ jet energy calibration 
this important systematic is already covered by the un- 
certainty of the fit result. The residual effect of the jet 
energy scale is estimated by varying the pr and rj depen- 
dence. The uncertainty of the jet energy scale of &-jets 
from various sources is considered separately. The uncer- 
tainty of signal modelling is determined from comparison 
of Pythia and Herwig and by varying the amount of 
initial and final state radiation produced in the parton 
shower. Further uncertainties include an estimate of the 
background modelling and treatment in each method. The 
DLM has the uncertainty due to initial and final state ra- 
diation as the most important single contribution, while in 
the MEM and the MTM the residual jet energy scale and 
the difference between generators used for signal simula- 
tion are the two leading contributions. In MEM and MTM 
colour reconnection effects [T^ITSlITi] are estimated and 
yield significant contributions of 0.56 GeV and 0.37 GeV, 
respectively. 

The most recent results of the various Matrix Element 
like methods of CDF were obtained at different integrated 
luminosities. All results yield in-situ jet energy scale cor- 
rections consistent with the default calibration used in 
CDF. 



MEM: 3.2fb~^ [92]: m^ = 172.4 ± 1.4 ± 1.3 GeV 
DLM: 1.7fb"^ [M]: m^ = 171.6 ± 2.0 ± 1.3 GeV 
MTM: 5.6 fb"^ [MIISS]: mt = 173.0 ± 0.9 ± 0.9 GeV 



(27) 



where the first uncertainty is the statistical one and in- 
cludes the overall jet energy scale uncertainty, the second 
uncertainty is the systematic one. The two dimensional 



representation of these results which simultaneously de- 
termine the jet energy scale are shown in Fig. 1131 These 
measurements of course use (partially) the same data, thus 
for the combination only the most precise one is used, tmtil 
their correlation is determined. 



2.2.4 Decay Length and Lepton Momentum Methods 

The top quark mass measurements described so far, aimed 
to use the maximal available information to yield the 
best statistical uncertainty for the given number of events. 
However, all these measurements have significant uncer- 
tainties from the jet energy scale. Alternative observables 
that have little or no dependence on the jet energy scale 
are therefore an interesting complement. 

In the lepton plus jets channel of top quark pair de- 
cays CDF has investigated two observables: the trans- 
verse decay length of b tagged jets and the lepton trans- 
verse momentum. One analysis combines the mean val- 
ues of both observables [101[|102| while the other uses the 
mean and distribution of the lepton transverse momentum 
alone IT03l. 



Two observable mean value method For the method us- 
ing both observables the data selection requires one iso- 
lated lepton, missing transverse momentum and at least 
three jets. For events with exactly three jets two of the jets 
need to be identified as b jets with CDFs secondary vertex 
algorithm [69,. For events with four or more jets only one 
needs to have such a reconstructed secondary vertex. 

The expected sample composition is determined from 
a combination of data and simulation. Top quark pair 
signal is simulated using Pythia with CTEQ5L parton 
densities for various nominal top quark mass values. Also 
single top quark samples are generated with various nom- 
inal top quark mass values. The dominating background 
VF+jets is simulated with Alpgen-|-Pythia, the multijet 
background is modelled from a data sample with modified 
lepton selection criterion. All simulated events are passed 
through the full CDF detector simulation and reconstruc- 
tion. 
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Fig. 14. Left and middle: Calibration curves for the expected mean decay length, {I/2d), and mean lepton momentum, (pt), 
as function of the nominal top quark mass (left and middle). Right: The likelihood for the combined measurement [102) . 



For the top quark mass determination the transverse 
decay length Liu of the h tagged jets is measured with 
the vertex algorithm used for the selection above. The ac- 
tual measurement is performed using the means of the 
observed decay lengths, {Lid), and the transverse lepton 
momentum, (p^). The expected contribution from back- 
ground is taken from the described background model. Be- 
cause the observables chosen are sensitive to the details of 
the event kinematics the signal simulation is corrected to 
yield parton distributions consistent with CTEQ6M [36 . 
In addition the total contribution from gluon fusion is cor- 
rected as function of the top quark mass. The performance 
of the simulation for determining the decay length L2D of 
the 6 tagged jets is calibrated in a dijet control sample. 

With these corrections the expected means of decay 
lengths and transverse lepton momentum and their ex- 
pected statistical spread are determined as function of the 
nominal top quark mass using ensembles of pseudo-data. 
The curves obtained are fitted by quadratic polynomials 
to obtain smooth curves. Figure [U (left, middle) shows 
the mean values (^20) and (p^) and the resulting top 
quark mass ranges from the individual observables. The 
combined result is obtained from the likelihood to find 
the observed mean values at various nominal top quark 
masses, c.f. Fig. [14] (right). 

The dominant systematic uncertainties for the mea- 
surement from (p^) arises from the uncertainty on the 
modelling of initial and final state radiation in the sig- 
nal simulation, from the lepton momentum scale and the 
shape of the background description. In the measurement 
from {L2d) the uncertainty of the data to Monte Carlo 
correction for the decay length dominates the uncertain- 
ties. The jet energy scale gives a non-negligible contribu- 
tion to the measurement from the average decay length 
through effects on the event selection. 

In 1.9 fb~^ of data CDF measures the top quark mass 
from the mean decay length, (^20), and the mean lepton 
momentum, (p^) to be |102| 

mt - 170.7 ± 6.3stat ± 2.6syst GeV , (28) 

where the jet energy scale uncertainty is included (as a 
small) contribution to the systematic uncertainty. The 
larger statistical uncertainty (compared to the Matrix Ele- 
ment methods) prevents a significant weight to the current 
combination of top quark masses |90] . 



Lepton momentum shape method The analysis that 
concentrates on the transverse lepton momentum 



alone |103| is based on events with a single isolated lepton, 
missing transverse energy and at least four jets. At least 
one of the jets is required to be identified as 6-jet. 

Top quark pair signal events are simulated for var- 
ious assumed top quark mass values. Backgrounds due 
to iy-|-jets, .Z-l-jets diboson and single top quark produc- 
tion are described with simulation. Estimates for multi- 
jet background with fake leptons and the normalisation 
for VF-|-light flavour jets are derived from data. These 
signal and background models are used to determine 
parametrised templates for the distribution of the lepton 
Pt spectrum as function of the top quark mass. To ex- 
tract the top quark mass these templates are compared 
to the observed data using an unbinned maximum likeli- 
hood fit. Leptonic Z boson decays are used to calibrate 
the lepton pT ■ 

The uncertainty of the final result is dominated by 
the statistical uncertainty. The systematic uncertainties 
are dominated by effects related to signal and background 
modelling. The uncertainty on the fake lepton description 
contributes with ±1.8 GeV; differences between generators 
for the signal modelling contribute with ±1.8 GeV to the 
systematic uncertainty. 

In 2.7 fb~^ of CDF data the top quark mass deter- 
mined from the shape of the lepton momentum distribu- 
tion is [ins] 



mt 



176.9 ±8.0stat±2.7,yst GeV 



(29) 



This result has been combined with a corresponding study 
of dilepton events [104 , c.f. section B. 3. 31 to yield |105| 



mt 



172.8 ±7.2,tat±2.3,y,t GeV 



(30) 



The large statistical uncertainty of these methods prevents 
a significant contribution to the world average from such 
measurements at the Tevatron. Its independent and low 
systematics will make them more relevant at the LHC. 



2.3 Dilepton Channel 

Due to the small branching fraction the dilepton channel 
has much fewer events, but the two charged leptons also 
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yield a much cleaner signature. For a measurement of the 
top quark mass the dilepton channel has the additional 
complication that the kinematics are under-constrained 
due to the two unmeasured neutrinos. Two of the six miss- 
ing numbers can be recovered from the transverse momen- 
tum balance, two more from requiring that the invariant 
mass of the charged lepton and its neutrino should be con- 
sistent with the W boson mass in the top and the anti-top 
quark decay, a fifth constraint can be obtained by forc- 
ing the top and the anti-top quark masses to be equal. 
Thus for a full reconstruction of an event including the 
top quark mass one constraint is missing. 



2.3.1 Weighting IVlethods 

To completely recover the event kinematics additional as- 
sumptions can be made on a statistical basis, i.e. by as- 
suming the distribution of one or more kinematic quanti- 
ties. For a given event top quark masses corresponding to 
certain values of these kinematic quantities are weighted 
by the probability that these kinematic values occur. The 
top quark mass reconstructed in the given event is then 
chosen such that it corresponds to the largest weight. This 
basic idea leads to various so-called weighting methods 
for measuring the top quark mass in dilepton events that 
mainly differ in the distribution that is assumed a priori. 
In the following the relevant methods recently used by the 
two collaborations are described. 
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Fig. 15. Probability density in the ml'''^°-mT2 plane for a dilep- 
ton signal at mt = 172 GeV (top) and background (bottom) at 
nominal jet energy scale ^7\ . 



CDF Neutrino Weighting One of these weighting meth- 
ods is the neutrino weighting method that is used in CDFs 
combined lepton plus jets and dilepton analysis [65,66,67;, 
[5H] . The lepton plus jets part of this analysis is described 
in Section [^Xn 

Dilepton events are selected by requiring two oppo- 
sitely charged leptons, missing transverse energy and ex- 
actly two energetic jets. In addition the scalar sum of the 
transverse momenta, Ht, is required to be greater than 
200 GeV. Events in which at least one of the jets was iden- 
tified as b jet are kept separately throughout the analysis. 

To reconstruct one top quark mass for each selected 
event the distribution of neutrino pseudo-rapidities is as- 
sumed a priori. The distribution is taken from top quark 
pair simulation and found to be Gaussian with an approx- 
imate width of one. The weight is computed in two steps 
for each assumed top quark mass. 

First, for fixed values of the neutrino pseudo-rapidities, 
rji, and yyp, the event kinematics are reconstructed using 
constraints on the assumed top quark mass and the nomi- 
nal W boson mass (ignoring the measurement of the trans- 
verse momentum). The weight of these 77 choices is con- 
structed as the x^ probability that the sum of the re- 
constructed transverse neutrino momenta agree with the 
measurement of the missing transverse momentum, '^t'- 



w{mt,rjy,rjy) 



^exp 



{i>T-pm-PT{^f 
2crf, 



(31) 



Here, the sum adds the four possible sign choices that 
occur in solving the above constraints and the two possi- 
bilities to assign the two measured jets to the 6 or 6 quark. 
Pxii) Slid Py(j) are the transverse neutrino momenta for 
the choice i and ctt is the experimental resolution of the 
measurement of the missing transverse momentum. 

The second step now folds the weights obtained 
for fixed pseudo-rapidities with the a priori probability, 
.P('7i/,'7p), that these pseudo-rapidities occur: 



W{mt)= Y] P{r]i^,r]i^)wimt,r]^,r]p) 



■Q,^,ri> 



(32) 



The top quark mass that yields the highest W{mt) is the 
reconstructed value, ?7i™'^°, for the event under considera- 
tion. In addition to this top quark mass the variable mT2 is 

used for further analysis. 771x2 — niin 



,« 



/ (1) (2)A 

I 777,^ , 777^ I 



where 77^^, are the transverse mass of the top and the anti- 
top quark, respectively. The minimisation is performed 
over all possible neutrino momenta consistent with the 
missing transverse energy. Previous versions of the anal- 
ysis instead used the scalar sum of the jets and lepton 
transverse momenta and the missing transverse momen- 
tum, Ht- 

Simulation is now used to determine the expected dis- 
tribution of 777,"'^° and 777,T2 ■ The backgrounds in the dilep- 
ton sample stem from fake events with a jet misidentified 
as lepton, from Dell-Yan and from diboson production. 
The dominating fake background is modelled from data. 
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Drell-Yan events are simulated with Alpgen+Pythia, 
diboson events with Pythia. The simulated events are 
passed through the CDF detector simulation and recon- 
struction. 

These simulations are used to generate probability 
density functions for the two observables, m^*^^" and mT2, 
as function of the nominal top quark mass, rrit, and the jet 
energy scale shift, Z\jes, see Fig. [15] From the probability 
densities a likelihood is constructed that is maximised to 
find the final result. The construction of the probability 
densities and the likelihood and the determination of the 
resulting measured top quark mass follows the lepton plus 
jets analysis described in Section [2 .2. II The dependence on 
the jet energy shift is only used for simultaneously fitting 
the lepton plus jets and dilepton events. Systematic uncer- 
tainties for the pure dilepton measurement are dominated 
by jet energy scale uncertainties, followed by the uncer- 
tainty of modelling colour reconnection effects. 

In 5.6 fb~^ of data CDF determines the top quark 
mass [68) in dilepton events using the neutrino weighting 
method and 771x2 to be: 



mt = 170.3 ± 2.0stat ± 3.1,yst GeV 



(33) 

A combined fit including lepton plus jet events and in-situ 
jet energy scale determination the dataset yields 



mt 



172.1 ±l.l,tat±0.9syst GeV 



(34) 



With the combined fit the dilepton result profits from the 
constraint of the jet energy scale. The overall precision is 
still dominated by the lepton plus jets result. 



CDF Neutrino Weighting Methods The neutrino 
weighting method uses the distribution of neutrino az- 
imuthal directions, cj), as an a priori distribution. CDF 
has analysed up to 2.9 fb^ of dilepton events to measure 
the top quark mass [il06j with this method. 

Events are selected by requiring one isolated well-iden- 
tified lepton, one oppositely charged isolated track, miss- 
ing transverse momentum and at least two jets. Vetoes 
are applied on events where the missing transverse mo- 
mentum is close to a jet, on cosmic events, conversions 
and Z boson events. 

For each value pair ((/)^, (f)p) on a grid of azimuthal di- 
rections the top quark mass is reconstructed using a con- 
strained fit that determines the top quark mass, mt, the 
lepton and (anti-)6 quark momenta using the measured 
lepton, track and jet momenta, the missing transverse mo- 
mentum and the constraints on the W boson mass and the 
equality of the top and anti-top quark mass. The x^ is de- 
fined similarly to Eq. pTj) . but uses Breit-Wigner distribu- 
tions rather than Gaussians for the mass constraints. For 
this fit two possibilities exist to assign the two measured 
jets to the b and anti-5 quarks. In addition the quadratic 
nature of the W boson mass constraints gives a fourfold 
ambiguity. Of the corresponding 8 top quark mass results 
at each {4>v,4>i,) pair only the one with the lowest x^ is 
kept. To arrive at a single mass value for the event un- 
der consideration the masses obtained at the various az- 
imuthal directions are averaged weighted with their x" 



probability. Only events with a weight of at least 30% of 
the maximum are considered in the average. 

The top quark mass is finally measured by compar- 
ing the distribution of mass values reconstructed for each 
event in data to templates in a likelihood fit. The tem- 
plates are derived from simulation. For the dilepton signal 



Herwig |107] is used at many nominal top quark mass 
values. The background is simulated with Pythia for 
the Drell-Yan background and Alpgen+Herwig for the 
fakes from VF-fjets production. The diboson background 
was simulated with Pythia and Alpgen+Herwig |"55] . 
Signal and background templates were parametrised to 
obtain smooth templates. The signal templates yield a 
smooth dependence on the nominal top quark mass. 

The likelihood for N events is built as the product 
over the probabilities that each event agrees with the sum 
of the signal and background templates, Pgig and Pbkg, 
weighted according to the signal and background contribu- 
tions, s and b. The likelihood contains a Poissonian term 
on the total number of expected events and a Gaussian 
constraint on the background contribution, 6, to be con- 
sistent with the a priori expectation, b°: 



„ii^^e-(^+^)(s + 6)^ 
L(mt) = e b I 



m 



N 

n 

i=l 



sPsigimj] mt) + bPbkgjmj) 
s + b 



(35) 



In addition a term, >Cparam, that describes the parametri- 
sation uncertainties of the template curves is included into 
the likelihood. Minimising the total likelihood yields the 
top quark mass results, nit, and estimates for the signal 
and background contributions, s and b. The statistical un- 
certainty is obtained by finding the top quark masses cor- 
responding to a log likelihood decrease of 0.5. Due to the 
inclusion of £param in the likelihood this includes the un- 
certainties from the template parametrisation. 

The systematic uncertainties are evaluated by measur- 
ing the top quark mass in ensembles of pseudo data gener- 
ated from simulation with modifications that refiect a one 
sigma change of the systematic under consideration. The 
systematics are dominated by the jet energy scale that 
has to be taken from the external calibration of the CDF 
calorimeter. It contributes ±2.9 GeV. The next to leading 
contributions from the background composition and the 
b jet energy scale contribute less than 20% of this uncer- 
tainty. 

In 2.9 fb~^ of CDF data with the described lepton plus 
track selection the method yields [106] 



m,t = 165.5 



+3.4 

— 3.3 stat 



±3.1,ystGeV 



(36) 



The method shows similar statistical precision as the neu- 
trino weighting method and similar sensitivity to the jet 
energy scale systematics. 



D0 Neutrino Weighting IVIethod D0 has applied the 
Neutrino Weighting method in up to 5.3 fb" of data j llOl 
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Fig. 16. Normalised weight distribution for an example simulated event with rrit — 175 GeV before (dashed) and after (full 
line) smearing lepton and jet momenta [108| (left). Example of the template histograms (middle) and functions (right) of /i,„ 
vs. (Tu, evaluated at rrit = 170 GeV for the efi channel |109j . 



110811109] ■ Events are selected by requiring two isolated, 
oppositely charged, identified leptons or an isolated, iden- 
tified lepton and an isolated track of opposite charge. In 
addition the events are required to have two energetic jets 
and missing transverse energy. For lepton plus track events 
at least one of the jets needs to be identified as 6-jet. Ve- 
toes on the lepton pair (or lepton plus track) invariant 
mass and the scalar sum of transverse momenta are ap- 
plied to reject Z-|-jets and other backgrounds depending 
on the identified lepton types. 

The event kinematics is reconstructed by assuming the 
neutrino rapidity distribution which according to simula- 
tions is expected to be Gaussian. Scanning through the 
range of possible top quark masses an event weight is com- 
puted as function of the top quark mass. 

First, for given values of the neutrino rapidities and 
the top quark mass a constrained fit is performed to de- 
termine the momenta of the top quark decay products, 
hhl'^i~vv. As constraints the W boson mass and the as- 
sumed top quark mass are used, ignoring the measured 
values of the missing transverse momentum. An individ- 
ual weight, w(rnt,r]^,r]jy), is computed from a x^ term that 
compares the sum of neutrino momenta in the transverse 
plane with the measured missing transverse momentum, 
c.f. Eq. (l5lT) . Then the a priori Gaussian distribution of the 
neutrino rapidities are folded into a total weight, Wi{mt), 
by adding the w{mt,r],y,r]j^) at 10 values of the neutrino 
rapidity with appropriate unequal distance. 

To take detector resolution for jet and lepton energies 
into account the determination of Wiirrit) is repeated for 
a number of jet and lepton momenta fluctuated accord- 
ing to their experimental resolution. The weight averaged 
over these fluctuations, W{mt) — {Wi{mt)), shows amuch 
smoother distribution and yields fit results for a wider 
range of top quark masses, see Fig. [12] (left). 

To determine the top quark mass the mean of the 
weight distribution and its variance are computed for each 
individual event: 



mtW{mt) dmt 



m'^W{mt)dmt - ^i^ 



(37) 



Compared to using ^^j alone, including a^ to the following 
extraction of the top quark mass yields an 16% improve- 
ment in the statistical uncertainty. Previous analyses of 
D0 used even more detailed information about W{mt) 
and thus the statistical information was exploited slightly 
better but at the cost of higher complexity |1101I109) . 

The distribution of ^w and a^ values in data is 
now compared to templates derived from simulation. 
Top quark pair signal events were generated with Alp- 
GEN-|-Pythia for various nominal top quark masses. 
Background contributions from Z/7-l-jets are simulated 
using Alpgen-|- Pythia, diboson events with Pythia. 
All simulated events are passed through the full D0 de- 
tector simulation and reconstruction. For Z/7-l-jets events 
with Z/7 — >■ e^e" or /i^/^^ the amount of fake miss- 
ing transverse momenta and of fake isolated leptons or 
tracks is derived from control samples in data and used 
in the normalisation of the samples. Template histograms 
are obtained from the above signal and background es- 
timates. The observed data distributions are compared 
directly to the template histograms or are compared to 
parameterised smooth functions fitted to the templates, 
see Fig. [16] (right) for an example. This yields two meth- 
ods, a histogram based and a function based method. 

The final top quark mass result is extracted by max- 
imising a likelihood as function of the number of signal 
and background events as well as the top quark mass. The 
likelihood describes the agreement of the data with the 
templates, a Gaussian constraint on the expected amount 
of background and a Poisson term for the total number of 
observed events as in Eq. (|35p. 

The performance and precision of the method is tested 
on a large number of pseudo data with known nominal top 
quark mass composed from the signal and background 
models. Calibration curves of the nominal vs. the mea- 
sured top quark mass are obtained by comparing the av- 
erage of the measured top quark masses on many such 
pseudo data with their nominal top quark mass. The ob- 
served offsets of around 1 GeV and deviations of the pull 
distribution from the normal distribution are corrected 
for. 



Daniel Wicke: Properties of the Top Quark 



21 



As for the other dilepton measurements the dominat- 
ing systematic uncertainty stems from the energy scale un- 
certainties for the two b jets in the event. In this analysis 
it contributes about ±1.5 GeV. Sub-leading uncertainties 
arise from limited template statistics and from the dif- 
ference between modelling the top quark pairs with pure 
Pythia and using ALPGEN-t-PYTHiA. 

The histogram based method and the method using 
parametrised templates show a correlation of only 85%. 
They are thus averaged using the BLUE method. Using 
dilepton events selected in Ifb^^ of data D0 determines 
the top quark mass to be |109j 

(38) 



eight values of Wi{mt) may exist. The total event weight 
is computed as the sum of these. 



mt = 176.2 ± 4.8stat ± 2.1syst GeV 



An updated result including up to 5.3 fb^ of data and 
combining with the Matrix Weighting method is described 
at the end of the next section. Despite the large systematic 
uncertainty from the jet energy scale the measurement is 
statistically limited at this luminosity. 



D0 Matrix Weighting Method Another method to mea- 
sure the top quark mass in dilepton events is the Matrix 
Weighting method. This has been applied by D0 in Run II 
to 5.3 fb~ of data [109j . The event selection requires two 
identified oppositely charged leptons, large missing trans- 
verse momentum and at least two jets. As above, vetoes 
on the lepton pair invariant mass and the scalar sum of 
transverse momenta are applied to reject Z-|-jets and other 
backgrounds depending on the identified lepton types. 

To reconstruct the top quark mass in each event this 
method assumes the distribution of lepton energies in the 
top quark rest frame lllj and the parton density func- 
tions. First the kinematics is reconstructed as function of 
the top quark mass. At each assumed top quark mass a 
kinematic fit is applied to reconstruct the top quark decay 
products from the measured quantities. The two leading 
jets are identified with the b quarks, the missing trans- 
verse momentum is required to be consistent with the sum 
of neutrino momenta in the transverse plane. In addition 
kinematic constraints from the W boson mass and the as- 
sumed top quark mass are applied. With the reconstructed 
kinematics a weight as function of the top quark mass is 
computed: 



w,{mt) = fixi)fix2)P{E*,+ \mt)PiE*-\mt) 



(39) 



Here the quark parton densities, f{x), at the quark and 
anti-quark momentum fractions, xi and X2, are explic- 
itly included in the weight. For the neutrino weighting de- 
scribed above they are included implicitly in the expected 
7] distribution. P{E'^± |r7it) are the probabilities that (given 
the hypothetical top quark mass, mt) the reconstructed 
energy of the lepton £^ in the rest frame of the corre- 
sponding top quark is E'^± . The distribution of these en- 
ergies is taken from the top quark pair production and 
decay matrix element [lllj . 

In each event there are up to four possibilities to solve 
the constraint from the W boson mass and two possibili- 
ties to assign the two jets to the two b quarks, thus up to 



Wiimt) = y^Wj(mt) 



(40) 



As for the neutrino weighting method resolution effects 
are included by recomputing the Wi with the measured 
quantities fluctuated according to their resolutions and 
averaging the results: W{mt) = (Wi(mt)). In each event 
the top quark mass which maximises this smeared weight, 
W{mt), is used as the estimator of the top quark mass for 
this event. 

The distribution of the top quark mass estimators ob- 
served in data is compared to templates obtained from 
simulation for a range of nominal top quark mass values. 
The simulations use Pythia to generate top quark pair 
signal and diboson background events. Backgrounds from 
Z-l-jets are generated with Alpgen-|-Pythia. The gener- 
ated events are passed through full D0 detector simula- 
tion and reconstruction. For the comparison and to extract 
the final top quark mass a binned likelihood is used. 

The method is calibrated by using ensembles of pseudo 
data constructed from the simulated events for various 
nominal top quark masses. The obtained bias in the deter- 
mination of the central result and its statistical error are 
corrected for. Ensemble tests are also used to determine 
the effect of systematic uncertainties. As for the other 
dilepton mass measurements the systematic uncertainty 
is dominated by uncertainties related to the jet energy 
scale. They account for ±1.2 GeV. 

In 1.0 fb~ D0 measures the top quark mass using the 



Matrix Weighting method in dilepton events to be [1121 

[inni. 

mt = 173.2 ± 4.9,tat ± 2.0,yst GeV . (41) 

With respect to the D0 neutrino weighting this Matrix 
Weighting method yields a slightly worse statistical pre- 
cision, but has a smaller dependence on the jet energy 
scale. D0 has averaged the two weighting methods de- 
scribed taking correlations into account using the BLUE 
method [109j . The combined result updated to 5.3 fb~^ 
yields [TT^ : 



mt = 173.3 ± 2.4stat ± 2.1syst GeV 



(42) 



As different events contribute to the two methods to a 
different amount, this allows a significant improvement of 
the statistical uncertainty over the individual results. 



2.3.2 Matrix Element Methods 

An alternative to the weighting methods of measuring the 
top quark mass in dilepton events is the Matrix Element 
method described above already for the lepton plus jets 
channel, see Section 12.2.31 Instead of assuming individual 
distributions to add effective constraints, the full knowl- 
edge from the matrix element is used to check the agree- 
ment of the measured events with different top quark mass 
assumptions. 
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CDF Matrix Element Method CDF has applied the Ma- 
trix Element method to the dilepton channel in several 
analyses P~14,115,116|. The analyses require two oppo- 
sitely charged leptons, missing transverse energy and two 
or more energetic jets. 

The probability to observe a given event depends on 
the top quark mass and is computed from folding the ma- 
trix element with parton densities and detector resolu- 
tions along Eqs. (I^TI) and (1^^ . In a dilepton sample of 
course the integration needs to include the unmeasured 
momenta of two neutrinos. The transfer functions that 
connect the measured quantities with the parton momenta 
correspondingly contains terms for only two jets. The mo- 
menta of the charged leptons are kept unsmeared. 

In the 1 fb~ analysis J115I an additional contribution 
to the transfer function was used to connect unclustered 
energy and energy from additional sub-leading jets with 
the transverse momentum of the top quark pair system. 
This more accurate treatment was found to yield better 
statistical power, but also larger systematics. In the recent 



update to 1.8 fb~ |116| such a term is no longer used, 
instead events with more than two jets are rejected. 

The matrix element used for top quark pair events cor- 
responds to the quark annihilation process. Background 
contributions are computed using the matrix elements for 
Z/^+jets, WW+jets and M^-|-3jets. The contributions are 
added according to their expected fractions in the selected 
dataset, depending on the number of identified 6-jets in 
the event. 

The probability for the full sample, P{mt), is the prod- 
uct of the event by event probabilities. As the jet en- 
ergy scale can not be constrained in dilepton events the 
external nominal CDF energy scale is used. The mean, 
j^raw _ JijitP{mt) drrit, is used as the raw measured 
mass. The result is calibrated using ensembles of pseudo 
data with varying nominal top quark masses. The pseudo 
data are modelled using Herwig for signal events, Alp- 
GEN-I-Pythia for Z'/7-|-jets and pure Pythia for diboson 
samples. Misidentified leptons are modelled from data. 
The observed bias and pull width observed by applying 
the above procedure on the ensembles are used to correct 
the raw mass measurement and its statistical uncertainty. 

The systematic uncertainties are again dominated by 
the jet energy scale uncertainty. In the 1.0 fb~^ analy- 
sis [115j the additional inclusion of additional jets is found 
to increase the dependence on the jet energy scale by 15%. 
In the 1.8 fb~^ analysis [116] the jet energy scale accounts 
for a ±2.6 GeV mass uncertainty. The next leading un- 
certainties in this analysis stem from simulation statistics 
for the background description, which of course can be 
addressed, and the uncertainty in the calibration. Differ- 
ences of signal simulation with Pythia vs. Herwig give 
only a marginally smaller contribution. 

In 1.8 fb~^ of dilepton events CDF measured the top 
quark mass with the Matrix Element method as [116] 



mt = 170 A ± 3.1stat ± 3.0,y,t GeV . 



(43) 



D0 Matrix Element Method D0 has applied the Matrix 
Element method to measure the top quark mass in dilep- 
ton events to 3.6 fb^^ [117J. The event selection requires 
one isolated electron and one isolated muon of opposite 
charge and at least two jets. 

The probability to observe a given event is computed 
as function of the top quark mass by folding the matrix el- 
ement with parton densities and detector resolutions along 
Eqs. (PT|) and (f^^ . The integration includes the addi- 
tional unknowns due to the second unseen neutrino and 
the unknown transverse momentum of the top quark pair 
system. The weight function, W{x,y), contains terms to 
describe the detector resolution for two leptons and two 
jets. The two possible jet to (anti-)6-quark assignments are 
summed. For the signal description D0 uses the matrix 
element of quark annihilation to top quark pairs and the 
subsequent decay. The backgrounds are all described with 
the matrix element for Z-|-jets process with Z — >■ t'^t'. 
Special transfer functions were used to relate the measured 
electron and muon momenta from a tau decay to the tau 
momentum. 

The likelihood for the observed data sample as func- 
tion of the assumed top quark mass, C{mt), is the product 
of the event probabilities, c.f. Eq. (|24)) . As dilepton events 
cannot constrain the jet energy scale, the externally de- 
termined nominal jet energy scale is used. The measured 
top quark mass is the one that maximises the sample like- 
lihood C{mt). Its statistical uncertainty is taken from the 
masses that yield an C(mt) decreased by half a unit from 
the maximum. 

The performance of the method is evaluated using en- 
sembles of pseudo data that are created from simulation. 
Top quark pair signal and .Z-l-jets background events were 
generated with Alpgen-|-Pythia, diboson events are gen- 
erated with pure Pythia. The samples were normalised to 
theoretical cross-sections. All simulated events are passed 
through the full D0 detector simulation and reconstruc- 
tion. The bias and widths of the pull distribution observed 
in the ensemble tests are used to correct the above mea- 
sured numbers. 

Systematics are dominated by uncertainties related to 
the jet energy scale of the two b jets in the event. They 
account to about ±2 GeV of the total systematic uncer- 
tainty. 

D0 has applied the Matrix Element method to deter- 
mine the top quark mass to the Runlla (1.1 fb^ ) and 
Run lib (2.5 fb^ ) run periods separately and combined 
the individual results using the BLUE method ^89_ with 
the error categories used by the Tevatron Electroweak 
Working Group [SU] (see also subsection 12.7! below) : 



mt = 174.8 ± 3.3,tat ± 2.6,yst GeV 



(44) 



With the increased statistics the method of J115j that 
yields a statistical improvement at the cost of increased 
systematics is no longer applied. 



Combinations with the other D0 results in the dilepton 
are available in [117) also. 



2.3.3 Lepton Momentum Method 

The so far described methods of determining the top quark 
mass in the dilepton channel, suffer from uncertainties of 
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The combination with the corresponding study of semilep- 
tonic events [103', c.f. section [2.2.41 yields [105) 



fb of data in 



100 150 200 250 300 

Transverse momcnlum (GcV/c) 



Fig. 17. Distribution and fit of lepton transverse momenta 
observed in 2.7 fb~^ of CDF data using dileptonic events. For 
fitted contribution of the background is shown in black. The 
inset shows the corresponding likelihood as function of the top 
quark mass [104) . 



the jet energy scale. This uncertainty dominates here the 
uncertainties much more than in the lepton plus jets chan- 
nel, where the in-situ cahbration allows to constrain the 
jet energy scale. Measurements with lepton based quanti- 
ties thus yield complementary results. 

CDF has performed an analysis of 2 
the dilepton channel using the transverse momentum of 
the leptons to extract the top quark mass |104j . Events 
are selected by requiring two isolated leptons of opposite 
charge and large missing transverse momentum. For the 
signal sample at least two jets are required, one of which 
needs to be identified as 6 jet. 

The background in this selection is dominated by 
fake signals from multijet events, which is estimated 
from events with same charge leptons. Contributions 
from diboson and Z/j events are simulated with Pythia 
passed through full detector simulation and reconstruc- 
tion. Events oi W + 6&-|-jets are generated with Alp- 
gen-I-Pythia. 

To extract the top quark mass from the selected events 
the distribution of the lepton transverse momenta is em- 
ployed. The expected distribution of the leptons' trans- 
verse momenta is parametrised as function of the top 
quark mass. A binned likelihood is used to compare the ob- 
served data to parametrised templates. The observed dis- 
tribution and the best fit templates are shown in Fig. [171 
Systematic uncertainties of this measurement are domi- 
nated by effects related to signal simulation. The choice of 
the generator and the uncertainty in modelling initial and 
final state radiation contribute with 1.5 GeV and 1.3 GeV, 
respectively. The momentum scale contributes with only 
0.7 GeV. 

In this preliminary result based on 2.7 fb^ of data the 
top quark mass determined from the shape of the lepton 
Pt distribution is |104) 



mt = 172.8 ± 7.2,tat ± 2.3syst GeV 



(46) 



mt 



154.6 ± 13.3stat ± 2.3,y,t GeV 



(45) 



The statistical power of these methods is not sufficient to 
yield significant contributions to the average from mea- 
surements at the Tevatron. Due to its independent and 
low systematics it will become more relevant at the LHC, 
though. 



2.4 All Hadronic Channel 



The full hadronic channel is the most difficult of all top 
quark pair decay modes. The channel has the advantage 
of the largest branching fraction and fully reconstructed 
events, but the huge background from multijet production 
makes it difficult to separate signal from background. 

D0 has published a measurement of the top quark 
mass using the all hadronic decay channel in Run I using 
the template method [118) , but not yet published a result 
with Run II. 



2.4.1 Template IVlethod 



CDF has applied a template method with in-situ jet en- 
ergy calibration to 2.9 fb^ of data |119j . The event se- 
lection requires between six and eight jets and vetoes on 
large missing transverse momentum and on identified lep- 
tons in the event. To further suppress the huge multijet 
background without top quarks a neural network is em- 
ployed which uses 13 input variables. These input vari- 
ables include event shape observables like sphericity, the 
minimal and maximal two and three jet invariant masses 
and observables based on the transverse jet energies and 
the jet angles. Only events above a threshold value of the 
neural network output are accepted. The threshold has 
been optimised for the most precise top quark mass result 
and depends on the number of b tagged jets in the events. 
Finally, at least one of the six leading jets needs to be 
identified as b jets using CDFs secondary vertex tagger. 

The background remaining after the selection is esti- 
mated from data. Jet by jet tag rate functions are mea- 
sured in events with exactly four jets and parameterised 
as function of the jet transverse energy, the track multi- 
plicity and number of vertices in the event. These jet tag 
rates are applied to the selected data with six to eight 
jets, but before requiring b tagged jets. Corrections for 
the fact that heavy flavour quarks are produced in pairs 
and for the signal content in the pre-tag samples are ap- 
plied. This background model was verified on data with 
the neural network cut inverted. 
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Fig. 18. Signal templates for various nominal top quark 
masses obtained for the sub-sample of events with exactly one 
h tagged jet. The lines correspond to the fitted functions used 
in the likelihood computation [119) . 



events from top quark pair signal simulation for various 
nominal top quark masses and to the described back- 
ground background estimate. For the background estimate 
the reconstructed top quark and W boson masses are en- 
tering the distributions with weights computed from the 
tagging probabilities and the corrections described above. 
Signal templates are created not only for various top quark 
mass values but also for a range of jet energy scale shifts, 
^.lES- To obtain smooth signal templates the distribu- 
tion obtained at discrete values of the nominal top quark 
mass and jet energy scale shift are fitted by functional 
forms. See Fig. [TS] for examples of such templates and 
their parametrisation. 

Now the results obtained in data are compared to these 
templates with a likelihood that is maximised with respect 
to the top quark mass, m^, the jet energy scale shift, Z\jeS; 
and the number of signal and background events. The like- 
lihood consists of a term describing the agreement with the 
W boson mass templates and of terms for the subsamples 
containing one or more than one h tagged jets, £i and C2'- 



The top quark mass is reconstructed in each selected 
event by performing a constrained kinematic fit to the six 
leading jets. The quark energies are allowed to vary within 
the experimental resolution of the measured jet momenta. 
The two pairs of light quark momenta are constrained 
to the W boson mass within the W boson width. The 
invariant mass of the two triples of quarks from the top 
and anti-top quark decays are constrained to agree with 
the top quark mass within the nominal top quark width, 
which is assumed as Ft = 1.5 GeV. The top quark mass 
and the quark energies are free parameters of the fit. The 
corresponding fit x^ can thus be written as follows: 



^ w 



+ 



(<J- 
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^ w 
2 (Mgj-mr°) 



Pt ) 



<yt 



(47) 



Of the possible jet to parton assignments required for the 
above fit, only those that assign b tagged jets to b quarks 
are used. With six jets this still allows thirty permutations 
for events with one identified b jet and six permutations in 
presence of two tagged jets. The top quark mass obtained 
from the jet-parton assignment that yields the lowest x^ 
is used as the reconstructed mass for the event under con- 
sideration. To obtain a handle on the jet energy scale, also 
the W mass is reconstructed in each event. This is done 
by repeating the above fit with also the W boson mass 
as a free parameter. Again the jet-parton assignment with 
the lowest y^ is chosen for further analysis. 

The distribution of these mass values reconstructed 
in each of the data events is compared to templates for 
various nominal top quark masses. These templates are 
constructed by applying the above fitting procedure to 



£(mf,Z\jEs) = exp 



/\2 
2crn 



£i(mi,Z\jEs)>C2(mt,Z\jEs)- 

The individual terms are constructed very similar to the 
CDF template method lepton plus jets, see Eq. ([T^. 

To estimate the performance of the method, the pro- 
cedure is applied to ensembles of pseudo-data for various 
nominal values of the top quark mass and of the jet en- 
ergy scale shift. The pseudo-data are constructed from 
the simulated signal and the data based background tem- 
plates. Then the nominal top quark mass for an ensemble 
of pseudo-data is compared to the average of the results 
obtained for each of the pseudo-data in that ensemble. 
Similarly, a calibration curve for the jet energy scale shift 
is constructed. These calibration curves show excellent lin- 
earity and only very minor overall offsets. In addition the 



w 2r 

UJ 

-> 

< 1.5 r 

1- 

0.5^ 


-0.5- 



-Ln(L/L ) Contours, 1 + > 2-tag events 



-1.5 




''^ y£ Calibrated Value; 
: — "'-"('-"-m.x) = ''■^ 

-Ln(L/L™^) = 2.0 
■Ln(L/L™'') = 0.5 



lIj 



_Lj 



_Lj 



166 168 170 172 174 176 178 180 182 184 

M,„p[GeV/c=] 

Fig. 19. Top quark mass and jet energy scale shift observed 
with the template method by CDF in 2.9 fb~^ using the all 
hadronic decay channel. The curves correspond to points of 
equal likelihood distance from the optimum [119] . 
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uncertainty estimate is verified with the width of the pull 
distribution and corrected accordingly. 

Systematic uncertainties are estimated as the average 
effect determined on ensembles of pseudo-data with sys- 
tematic effects included. The largest single uncertainty is 
estimated to stem from residual (mass and jet energy de- 
pendent) biases from the template parametrisation, not 
covered by the calibration procedure (to'4 GeV). Residual 
effects of the jet energy scale follow with ±0.5 GeV. This 
analysis also evaluates the effect of changing the under- 
lying event model from a default |120) to a new model 
including colour reconnections [7^[751[7t] and finds an un- 
certainty of 0.4 GeV. 

In 2.9 fb~ of data CDF determines the top quark mass 
with in-situ jet energy calibration from events in the all 
hadronic channel to be |119] 



174. 



1 ± 2.4stat-hJES_i:o syst GcV 



(49) 



The measured two dimensional likelihood is shown in 
Fig. [T21 The jet energy scale shift determined is consis- 
tent with CDFs nominal value, but has a smaller total 
uncertainty. 



2.4.2 Ideogram Method 



Also the ideogram method described for the case of lepton 
plus jet events in Section [2.2.2l is well suited to be applied 
in the all hadronic channel. 

CDF applied this method of measuring the top quark 
mass in up to 1.9 fb~^ of data [12111122] . The first steps of 
the analysis are very similar to the template based anal- 
ysis for the all-hadronic channel. The selection of events 
requires exactly six jets. A neural network is used to fur- 
ther suppress background from multijet production with 
no top quarks. For the signal sample at least two of the 
jets are required to be identified as b jets using CDFs sec- 
ondary vertex algorithm. 

The background contribution after this selection is 
computed by applying jet by jet tag rate functions on the 
data before the b tag requirement. The required tag rate 
functions are derived on events with four jets, where the 
signal contamination is negligible. 

Then a constrained kinematic fit is performed to si- 
multaneously determine the top quark mass and the W 
boson mass for the event under consideration. Transfer 
corrections are applied to correct differences between the 
jet energy inside the chosen cone radius of 0.4 and origi- 
nating quark. These transfer corrections optionally imple- 
ment a jet energy scale shift. The constraints of the fit are 
that the quark momenta yield identical top quark masses 
and identical W boson masses in the two decay chains of 
the event, c.f. Eq. (|47|) . The result for all possible jet to 
parton assignments which assign b tagged jets to (anti-)& 
quarks are kept for the further analysis. 



With these numbers in the ideogram method an event 
by event likelihood is computed that describes the prob- 
ability to observe the reconstructed top quark and W 
boson masses for the various jet to parton assignments 
given the true top quark mass and jet energy scale shift. 
This likelihood consists of probabilities, Pgig and Pbkg, 
corresponding to the signal and the background hypoth- 
esis, c.f. Eq. (IT5|) . The signal and background probabili- 
ties obtained for each of the jet to parton assignment are 
weighted by their x^ probability, Wi , which depends on the 
jet energy scale applied. In addition CDF adds a term to 
the signal probabilities to describe events with no correct 
jet to parton assignment of the selected six jets: 

N 
Psigix;mt,AjEs) = /nm^ Wi(^JEs)'S'(™t,i, Mh/.j, ZijEs) 



+ (1 - fnm)Snmimt,i, Mwa, ^JEs) 



N 



Pbkg(a;;/ijEs) = '^m{A3Es)B{mt,i, Mw,i) 



(50) 



i=l 



N is the number of possible jet parton assignments, i.e. 
6 for doubly tagged events and 18 for events with three b 
tagged jets. Z\jes the jet energy scale shift in units of the 
nominal jet energy scale uncertainty, mt^i and M\Y,i are 
the fit results for the jet parton assignment number i and 
1 — /nm is the fraction of events with no correct jet parton 
assignment. 

The signal probability, S, is written as the convolu- 
tion of the theoretically expected Breit-Wigner, BW, and 
the experimental smearing represented by a Gaussian, G. 
Here two Breit-Wigner functions are needed, one for the 
top quark and one for the W boson, and the Gaussian 
is two dimensional, constructed including correlations be- 
tween the top quark and the W boson mass extracted from 
the kinematic fit: 

S{mt,i,Mw,t,Ams) = 

G{fnt.^,Mw,^,m',M',ai) 

BW{m\mt)BW{M',Mw) dm'dM' . (51) 

The background probability, B, and the probability for the 
signal events with no matching, ^nm, are derived from sim- 
ulation. The final likelihood is then written as the product 
of the event likelihoods times a prior likelihood for the jet 
energy scale shift to describe the external calibration 



C{m.t,AjEs) = exp 



Z\2 



Y\^ Ccvtix;Tnt,AjEs) 



(52) 
and maximised with respect to the top quark mass, rrit, 
and the jet energy scale shift, zijES- 

The performance of this method is evaluated in en- 
semble tests. Ensembles of pseudo data are generated for 
a grid of nominal top quark masses and jet energy scale 
shifts using the Pythia generator and full CDF detec- 
tor simulation and reconstruction. The mean response of 
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Fig. 20. Top quark mass and jet energy scale shift observed 
with the ideogram method by CDF in 1.9 fb~^ using the all 
hadronic decay channel. The curves correspond to points of 
equal likelihood distance from the optimum [122) . 



the analysis on these ensembles is used to correct the top 
quark mass and jet energy scale shift obtained in data. 

Systematic uncertainties are evaluated on ensembles 
of pseudo data with the corresponding systematic vari- 
ation. A variation of initial and final state parameters 
of the simulation yields the dominating contribution of 
±1.2 GeV, followed by the difference obtained from com- 
paring Pythia to HERWiG(±0.8GeV). The residual jet 
energy scale contributes with ±0.7 GeV. 

In 1.9 fb~^ of data CDF measures the top quark mass 
from all hadronic events to be [122] 



mt ^ 165 ± 4.4stat+JES ± l-9syst GeV 



(53) 



Figure [20] shows the result with points of equal likelihood 
distance from the optimum in the top quark mass vs. jet 
energy scale shift plane. 



2.5 Top Quark Mass from Cross-Section 

As discussed in Section 12.11 the definition of quark masses 
is inherently ambiguous. The quoted mass values are often 
considered to be given in the pole mass definition. How- 
ever, this definition usually yields significant higher order 
corrections. For Monte Carlo simulations the order of the 
corrections included in the determination cannot be easily 
computed due to the presence of parton shower cut offs 
and modelling of hadronisation. 

Deriving the top quark mass from the top quark cross- 
section measurements avoids using the simulation for cal- 
ibration and allows the determination of the top quark 
mass using a well defined mass definition in an under- 
standable approximation. 

Both predicted and the measured cross-section depend 
on the true top quark mass. In the theoretical prediction 
this dependence stems from the change in phase space 
with varying top quark mass. For the experimental re- 
sults this dependence is introduced through the selection 
efficiencies which vary with the amount of energy available 
for the decay products. 



D0 has compared their cross-section measurement for 
up-to 1 fb~ of data to theoretical predictions as function 
of the top quark pole mass for various theoretical approx- 
imations [TMlll25[ll2(ij . 

For the extraction of the top quark mass the depen- 
dence of the theoretical predictions and the experimental 
cross-section results are parametrised with a polynomial. 
Both the theoretical and the experimental cross-section 
uncertainties are assumed to be Gaussian to build a like- 
lihood that then allows to find the top quark pole mass 
that yields the best agreement and its uncertainty. This 
assumption is justified also for the theoretical uncertainty 
as it is dominated by PDF uncertainties. 

The measured cross-section is compared to several 
theoretical predictions in Fig. [5TJ A pure NLO predic- 
tion |127j . approximate next-to-next-to leading order [MJ 
155] and a next-to-leading order with resummation of lead- 
ing and next-to- leading soft logarithms T? . All are evalu- 
ated with the CTEQ6.6 parton density distribution [128) . 

In this method statistical and systematic uncertainties 
are already combined. Using the combined cross-section in 
■^±jets, dilepton and r+lepton, D0 derives in 1 fb~^ [126) : 



NLO pTlfnS] 

NNLOapprox [3] 
NNLOapprox H 

NLO+NLL m 



mt = 165.5 ±6.0 GeV 
mt = 169.1 ±5.6 GeV 
mt = 168.2 ±5.6 GeV 
mt = 167.5 ±5.7 GeV 



(54) 



The results on the top quark pole mass agree well with the 
world average of direct measurement, but have a much 
larger uncertainty. As these numbers refer to the same 
data, the differences between the results of about ±2 GeV 
has to be attributed to theoretical differences and may be 
an indication of the sensitivity of the pole-mass definition 
to higher order corrections. 

Recently in a first determination of the MS-mass of the 
top quark has been presented 11291 



Based on the data of 
the same D0 cross-section result [126) a running top quark 
mass of mt = 160.0 ± 3.3 GeV is extracted using the ap- 
proximate NNLO prediction [33] . To compare this number 
to the pole mass result quoted above it has to be converted 
to the pole mass. This conversions strongly depends on the 
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Fig. 21. Comparison of the theoretical top quark pair cross- 
section to the experimental results (^-|-jets, dilepton and 
T-|-lepton) with their top quark mass dependence. The lines 
show parametrised dependence and the error bands [123] . 
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order used. In leading order no change occurs, while us- 
ing the NNLO formula yields rrit = 168.2 ± 3.6 GeV. (For 
simplicity the slightly asymmetric uncertainties were sym- 
metrised by the author.) It should be noted that the MS- 
mass determined from different orders of the MS cross- 
section calculation change by less than 1 GeV between 
leading order and next-to-next-to-leading order, while the 
corresponding pole mass of the same order changes by 
nearly 10 GeV. 



2.6 Modelling of Non-Perturbative Effects 

Modelling non-perturbative effects in simulations is a no- 
toriously difficult task. Beside the effect of hadronisation, 
which is believed to be well understood since LEP times, 
in hadron-hadron collisions multiple parton interactions 
yield a different kind of non-perturbative effects. The un- 
certainties on the top quark mass due to their modelling 
uncertainties have only recently been considered. 

In the context of hadron collisions at the Tevatron sev- 
eral efforts to tune the default model present in Pythia 
to data [120j resulted in several parameter tunes, known 
as Tune A, Tune DW, Tune D6, etc. These tunes required 
large values of parameters describing the colour correla- 
tions between the multiple parton interactions and the 
hard process. This is often interpreted as an implicit hint 
for the need of colour reconnections, i.e. a modification of 
the colour fiow obtained from the simulation of the hard 
and semi-hard process and parton showers. 

Newer versions of Pythia implement an extended 
model to describe the underlying event including 
variants of the Pythia model with explicit colour- 
reconnection [7^175] . These are tuned and their influence 
on the top quark mass measurements is studied with a toy 
top quark mass analysis [73ll74] . Compared to the variants 
of the default PYTHiA-model the new model and its vari- 
ants yields shifts in the extracted top quark mass of up to 
1 GeV. The greater portion, namely 0.7 GeV, of this shift 
can be assigned to differences between the different parton 
showers used in the old and the new underlying event mod- 
els. Less than 0.5 GeV is attributed to non-perturbative 
effects. It should be noted here that additional alternative 
models of colour reconnection suitable for hadron colli- 
sions exists [130.131] , but were not yet studied in the con- 
text of the top quark mass. 

While the smallness of the non-perturbative portion 
of the derived shift on the top quark mass confirms the 
prejudice that non-perturbative effects are expected at the 
order of (a few times) the confinement scale, ^qcd, the 
sizable shift due to parton shower details is very worrying. 
Only very recent top quark mass measurements include an 
uncertainty on the colour reconnection effects (see the pre- 
vious subsections) and confirm the estimated uncertainty 
due to modelling of non-perturbative effects. The greater 
shift due to the parton shower differences is still not in- 
cluded. 

It has recently been realised that the two types of mod- 
els differ mainly in their 6-quark jet energy scale and that 
the shape of 6-jets does not agree with the models using 



the pT-ordered parton shower |1321I133) . The simultane- 
ous determination of the light quark jet energy scale and 
the 6-quark jet energy scale as suggested by [13411135) may 
thus help to resolve the issue. 



2.7 Combination of Top Quark Mass Results 

The Tevatron experiments regularly combine their Run I 
and Run II top quark mass results to take advantage of the 
increase in statistical power. The most recent combination 
was performed in July 2010 [90] . 

For each measurement that enters the combination a 
detailed break down of errors is performed. Uncertainties 
that are believed to have correlations of one measurement 
with any other measurement of the same or the other col- 
laborations are separated so that these correlations can be 
taken into account in the combination. At this point only 
measurements from independent dataset or selections are 
used in the combination. The evaluation of partial correla- 
tions, as they would appear when using results form multi- 
ple methods on the same dataset and channel, are thereby 
avoided. The various uncertainty contributions are consid- 
ered to be either fully correlated or to be uncorrelated. 

The average is then computed with the best linear un- 
biased estimator (BLUE) f^ and yields ^: 



rat = 173.3 ±1.1 GeV . 



(55) 



Uncertainties on this estimation due to the approxima- 
tions of the procedure including a cross-check with re- 
duced correlation coefficients was estimated to be much 
smaller that 0.1 GeV. 



2.8 Mass difference between Top and Anti-top Quark 

According to the CPT-theorem local quantum field theo- 
ries require the mass of any particle to be equal to that 
of its anti-particle. A measurement of the mass difference 
between top and anti-top-quarks is a unique test of the 
validity of the CPT-theorem in the quark sector, as other 
quark masses are more difficult to assess due to hadro- 
nisation. Both Tevatron experiments have extended their 
studies of the top quark mass to determine this mass dif- 
ference. 



2.8.1 Matrix Element Method 

D0 extends the Matrix Element method using 1 fb^ of 
data \XM„- The analysis selects events with one isolated 
lepton, transverse missing energy and exactly four jets. At 
least one of the jets is required to be identified as 6-jet. 

The Matrix Element method used to measure the top 
quark mass as described in Sect. 12.2.'^ is extended by ex- 
plicitly keeping a separate dependence on the top quark 
mass, TOt, and the anti-top quark mass, rai in Eqs. ((2T|) to 
(IM)) . For this the leading order matrix element is rewrit- 
ten to explicitly depend on mt and mi. At the same time 
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the dependence on the overah jet energy scale factor, /jes, 
is dropped. This yields a likelihood C{mt ^ mi, f top) corre- 
sponding to Eq. (IM)) . The likelihood to observe a mass 
difference of Z\ = rrit — mi is obtained by integrating over 
all possible average mass values mavg = {m^t + mi)/ 2. 



C{A) = / Am, 



vg £{mt,mi, fl'°p\mt,mi)) 



(56) 



Here f^^^^{mt, m,i) is the fraction of top quark events fit- 
ted for the considered values of mt and mj . The measured 
mass difference A is the one that maximises this likeli- 
hood. 

As for the mass measurements the method is calibrated 
using pseudo-experiments. The pseudo-data in these tests 
are constructed from ti signal events and VF-|-jets back- 
ground. For the signal various values of mt and mi were 
generated with a modified version of Pythia that allows 
to set mt 7^ mf. It was confirmed that the mass differ- 
ence measured in the various pseudo datasets is very close 
to the generated value. The small deviation in the recon- 
structed value and the derived uncertainty are corrected 
for in the final result. 

Many uncertainties that are important for the mass 
measurement cancel in the determination of the top anti- 
top quark mass difference. This includes the uncertainties 
due to the jet energy scale and the uncertainty due to 
a difference in the detector response between light and 
b quark jets. Instead the dominating systematic uncer- 
tainty in the measurement of the mass difference is that 
of modeling additional jets in events of top quark pairs. 
This signal modeling uncertainty is estimated from data 
and simulation. Using simulation the amount of events 
with more than four jets is varied to agree with data. In 
addition, pseudo experiments are performed in which ap- 
propriate data events with more than four jet are added 
to the standard signal simulation with various amounts. 
In combination D0 finds an uncertainty of 0.85 GeV due 
to signal modelling. Because of the necessity to distin- 
guish top from anti-top quarks also the uncertainty in the 
lepton charge determination and the uncertainty in mod- 
elling the differences of the calorimeter response for b and 
b jets is evaluated. The latter gives a sizable contribution 
of 0.4 GeV in the D0 result, but its determination is lim- 
ited by the statistics of the corresponding simulation. Due 
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Fig. 22. Normalised likelihood distribution as function of the 
top to anti-top quark mass difference, A, as measured by D0 
in the electron and muon plus jets channels, separately [136) . 



to the negligence of multijet background in the likelihood 
computation and its calibration the multijet contamina- 
tion yields another contribution of 0.4 GeV. The total sys- 
tematic uncertainty is 1.2 GeV. 

D0 applies the described method using 1 fb~^ of data 
in the e-|-jets and the /x-|-jets channel separately. The cor- 
responding likelihood functions are shown in Fig. 1221 The 
two channels are then combined by a weighted average 
which yields 

Z\ = 3.8 ±3.7 GeV (57) 

in good agreement with the CPT-theorem expecta- 
tion of zero [136] . As a cross-check the likelihood 
C{m,t,mi,f^°p*-{m,t,mi)) is integrated over all possible 
mass differences and used to determine the standard top 
quark mass. From this the relative difference is obtained 
to be 



A/mt = (2.2 ±2.2)% 



(58) 



2.8.2 Template Method 



For the measurement of the top quark mass differ- 
ence [137) CDF has expanded the template method with 
full reconstruction of the top quark decay products de- 
scribed in Sect. 12.2.11 The event selection requires one 
isolated energetic lepton, missing transverse energy and 
at least four energetic jets in the final state. Events are 
categorised by the charge sign of the lepton and the num- 
ber of identified 6-jets. Only for events with more than one 
identified 6-jets more than four high px jets are allowed. 
In each event the mass difference is reconstructed with 
a fit constrained by the top quark decay kinematics. The 



reco 
t 



X^ for this fit is described by Eq. (fTTj) when replacing m' 
with 172.5 GeV±Z\/2. If the top quark decays hadronically 
and the anti-top quark leptonically, the positive sign ap- 
plies for term including M^qq and the negative sign for the 
term including Mhi^, . If the top quark decays leptonically 
instead, the signs are exchanged so that A represents the 
reconstructed mass difference in all cases: A — mt — mi. 

The x^ is minimised for all possible associations of 
jets to the four quarks requiring that identified 6-jets are 
only associated to &-quarks. For further analysis the values 
of Z\^^) and Z\'-^-' corresponding to the association with 
the smallest and the second smallest x^ are used. Events 
where the lowest x^ is larger than 3.0 (9.0) are rejected 
for events without (with) 6-tagged jets. 

Simulations are used to compute the expected signal 
contributions as function of the two reconstructed mass 
difference values for various nominal mass differences us- 
ing MadGraph+Pythia. The contribution of multijet 
background is modelled using data with loosened lep- 
ton identification criteria. The distribution of mass dif- 
ferences reconstructed in VF+jets background is modeled 
with Alpgen+Pythia. Its normalisation is derived from 
data. Further smaller backgrounds are taken from simula- 
tion normalised to the NLO cross-sections. 
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As for the mass determination (Sect. 12.2.11) the signal 
simulations are used to derive a probability density, P®'^, 
to reconstruct a pair A'^^\A^^'> for a signal event with 
a given nominal mass difference. These two dimensional 
functions depend on the lepton charge, which is the reason 
to separate samples of different lepton charge. Similarly a 
probability density for background events, P^^s^ jg com- 
puted. Then a likelihood for the six subsamples, i = 1...6, 
is computed as 



£,; = exp — 



2< 



n 



p 



6,:-P, 



bkg 



(59) 



with Si and bi being the number of signal and background 
events in the corresponding sample; 6f and ab^ are the 
corresponding background expectations and its error. 

The final likelihood, i.e. the product of all Ci, is avail- 
able only at discrete values of the nominal mass differ- 
ence. For the final maximisation of the likelihood it is in- 
terpolated with polynomial smoothing. CDF found their 
method to be bias free in pseudo experiments. A small cor- 
rection to the derived statistical uncertainty is applied. 

Systematic uncertainties are dominated by the signal 
modeling uncertainty of 0.7 GeV, which is derived by com- 
paring MadGraph to pure Pythia and by comparing the 
results of using the Herwig to the Pythia parton shower. 
The next to dominating uncertainties stem from the de- 
scription of multi hadron interactions and differences in 
the detector response for b and b quarks. The former is 
derived by verifying the agreement of the simulation with 
data as function of the number of primary vertices. The 
latter is determined by comparing the pt balance in di-6- 
jet events. These effects yield systematic uncertainties of 
0.4 GeV and 0.3 GeV, respectively. 

Using 5.6 fb~^ of data CDF finds the top quark to anti- 
top quark mass difference to be [137] 



Z\ = -3.3±1.4tat±1.0syst 



(60) 



which agrees with the CPT-theorem at a little less than 
two standard deviations. The analysis assumes explicitly 
an average (anti-)top quark mass of 172.5 GeV. 



2.9 Conclusions and Outlook to LHC 

In the Standard Model of particle physics the top quark 
mass is an a priori unknown parameter. The Tevatron ex- 
periments have employed many different techniques to de- 
termine its value. Emphasis has been put on the reduction 
of the experimental uncertainties. The Ideogram and the 
Matrix Element methods aim at the maximal utilisation 
of the available experimental information. The simultane- 
ous 'in situ' determination of the jet energy scale along 
with the top quark mass addresses the dominant system- 
atic uncertainty. Also the next dominant uncertainty, the 
6-jet energy scale, can in principle be determined in situ, 
however, more statistics is needed to achieve an improve- 
ment over the current uncertainty. 



These results constrain the Standard Model prediction 
for processes that involve top quarks either as real par- 
ticles or in virtual loops. Besides the real production of 
top quarks, electroweak precision data are sensitive to the 
value of the top quark mass. Even before the discovery of 
the top quark the electroweak precision data constrained 
the top quark mass 138J. Now such indirect values are 
compared with the direct measurements at the Tevatron 
to verify the consistency of the Standard Model [U01I139) . 
Figure [23] shows the direct measurement of the W bo- 
son and the top quark mass compared to the indirect top 
quark mass from electroweak precision measurements. It 
indicates the Standard Model expectation as function of 
the Higgs mass. Also the Standard Model agreement with 
the global electroweak results is usually computed. Fig- 
ure [M] shows the x^ fit as function of the Higgs boson 
mass which also determines the indirect constraints on 
the allowed Higgs boson mass. 

Implicitly, such comparisons assume that the exper- 
imentally measured top quark mass values are given in 
the pole mass scheme. This assumption, however, is not 
confirmed to the level of the experimental precision of 
the combined result from direct top quark mass measure- 
ments. Therefore conclusions based on such comparisons 
need to be drawn with great care. Hopefully, in the fu- 
ture it will be possible to evaluate the exact top quark 
mass definition used in the Monte Carlos and used in the 
calibration of top quark mass measurements. Any bias in 
the current usage of the value can then be corrected for 
and the quoted experimental uncertainties can be used in 
these comparisons unchanged. Until the mass definition of 
the simulations used for calibration can be determined, it 
will be an important task to develop top quark mass mea- 
surements that are based on well defined mass schemes. 
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Fig. 23. Mass of the W boson vs. the top quark mass. The 
smaller (blue) contour indicates the one 68% C.L. of the direct 
measurements, the larger (red) contour the indirect results. 
The expectation in the Standard Model forms diagonal lines 
that depend on the Higgs boson mass [90| . 
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Fig. 24. Quality of the agreement of precision measurements 
with the Standard Model prediction as function of the assumed 
Higgs boson mass after a global fit of the parameters to data. 
The yellow region indicates the direct exclusion limit. The cen- 
tral result is shown as a back line with blue band. Alternative 
fits are shown in other colours [90j . 



The LHC has started data taking in 2010 at an ini- 
tial centre-of-mass energy of 7TeV. First results show a 
tremendous quality of the data taken by the experiments 
and top quark pair production was already seen in about 
3pb"^ 140,141 . With the planned LHC luminosity of 

1 fb~ the LHC experiments will become competitive to 
the Tevatron experiments within 2011. Systematic uncer- 
tainties on the top quark mass using the methods devel- 
oped essentially at the Tevatron were expected to be below 

2 GeV even before the initial data taking [Ti^lTiHllTii] . 
With the performance shown in the first year of running 
a result comparable to the current Tevatron result can be 
expected. The statistical uncertainty in these results will 
still be noticeable. 

A truly negligible statistical uncertainty can only be 
expected after a successful LHC data taking at the de- 
sign centre-of-mass energy of 14 TeV or close. This en- 
ergy can be expected only after a long year shutdown. 
Only then the correspondingly enhanced production cross- 
section will open the door to methods applicable only at 
the LHC. The top quark mass determination involving 
leptonic J/tp decays from the 6-quark jet |1451I146) is ex- 
pected to be less dependent on the jet energy scale allow- 
ing to cross-check the in situ methods. In production with 
very high transverse top quark momentum the top de- 
cay products are clearly separated from the antitop quark 
decay products. These events are theoretically better un- 
derstood and allow a mass measurement in a well defined 
mass definition (6311147) . At a design LHC this might re- 
solve the puzzle of interpreting the current experimental 
top quark mass. 



3 Interaction Properties 

The Standard Model fixes the properties of top quark for 
all three interactions considered in the Standard Model. 
To establish that the top quark discovered at the Teva- 
tron is in fact the Standard Model top quark it is im- 
portant to verify the expected properties experimentally 
and to set limits on possible deviations. This subsection 
will consider interaction properties, i.e. measurements of 
top quark properties and possible deviations from the 
Standard Model that do not assume explicit presence of 
non-Standard Model particles. First measurements of top 
quark properties regarding the weak interaction shall be 
described. Then the verification of the electrical charge is 
summarised followed by properties of the top quark pro- 
duction through the strong force. Measurements that in- 
volve non-Standard Model particles are covered in Sec- 
tion H 



3.1 W Boson Helicity 

One of the first properties of the electroweak interaction 
with top quarks that was measured is that of the W boson 
helicity states in top quark decays. Within the Standard 
Model top quarks decay into a W boson and a b quark. To 
check the expected V — A structure of this weak decay the 
W boson helicity is investigated. Only left-handed parti- 
cles are expected to couple to the W boson and thus the 
W boson can be either left handed (— ) or longitudinal (0). 
For the known 6 and t quark masses the fractions should 
be /_ = 0.3 and /o = 0.7, respectively. The right handed 
(+) contribution is expected to be negligible. 

Depending on the W boson helicity (—,0,-1-) the 
charged lepton in the W boson decay prefers to align 
with the h quark direction, stay orthogonal or escape in 
the opposite direction. Several observables are sensitive to 
the helicity: the transverse momentum of the lepton, the 
lepton-6-quark invariant mass, M/(,, and the angle between 
the lepton and the h quark directions, cos 0* . For best sen- 
sitivity at Tevatron energies cos 9* is measured in the W 
boson rest-frame. 

All of these observables have been used to measure 
the W boson helicity fractions at the Tevatron. The most 
recent and thus most precise results use cos 9* as the ob- 
servable. 



CDF 



Analysis based on Mf^^ The lepton-6-quark invariant 
mass, A/;fc, was used in an analysis of approximately 
700 pb^^ [148|. Events with a lepton plus jets signature 
containing an isolated lepton, missing transverse energy 
and at least 3 jets with identified 6-jets were studied sep- 
arately for one and two identified jets. In addition events 
with a dilepton signature containing two identified leptons 
with opposite electrical charge, missing transverse energy 
and at least 2 jets are investigated. 
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Fig. 25. Distribution of the squared lepton-6-quark invariant mass measured by CDF compared to best fit expectation (left). 
Likelihood distributions as function of the V+A fraction fv+A for the individual and the combined channels. The shaded yellow 
region indicates the allowed 95% confidence range (right) [148] . 



For each selected event the squared invariant lepton-6- 
quark mass, M^f^, is computed. In lepton plus jets events 
with a single identified &-jet the computation is un- 
ambiguous, however, the identified 6-jet and the lepton 
are from the same top quark in only half of the cases, 
see Fig. [25] (left). For events with two identified 5-jets a 
two dimensional distribution of Afj| is constructed. One 
dimension being M;'^ computed with the higher energetic 
6-jet, the other with the lower energetic 6-jet. In the dilep- 
ton events the two dimensional histogram is filled twice, 
i.e. using each of the leptons. 

The contribution of background events in the lepton 
plus jets samples is modelled by Alpgen W + bb events 
and by multijet events obtained from a control data sam- 
ple. In the singly tagged sample multijet events contribute 
15%, in the doubly tagged sample they are neglected. 
For the dilepton sample the background is described by 
about 50% Z-l-jets from Alpgen, 30% W-l-jets with one 
jet misidentified as lepton using a single lepton sample 
and applying the misidentification rate. The final 20% of 
background are from diboson samples, WW and WZ. 

The signal expectation is simulated with Alpgen and 
Pythia assuming a top mass of rrit — 175 GeV for both 
V — A and V + A coupling to the W boson. A binned log 
likelihood fit is used to extract the fraction oi V + A cou- 
pling in the top decay, /+. The likelihood uses Poisson 
probabilities to describe the expected number of events in 
each bin of the M;-| distributions. The parameters of the 
Poisson distributions are taken from the described simu- 
lation and are smeared with nuisance parameters for top 
pair cross-section and the total background contribution 
in each of three samples. These nuisance parameters are 
constrained to their nominal values with Gaussian proba- 
bility distribution. 

The procedure was verified using a large number of 
pseudo experiments with different nominal contributions 
from V + A decays. The fit was found to be stable and 
unbiased. Systematic uncertainties are dominated by the 
uncertainty on the jet energy scale, followed by uncertain- 
ties of the background shape and normalisation as well as 
the limited Monte Carlo statistics. 



The likelihood distributions obtained with this pro- 
cedure are shown in Fig. [51] (right). The left-handed W 
boson fraction in the top quark decay is measured to be 
/+ = — 0.02 ± 0.07 in agreement with the Standard Model 
expectation of a negligible contribution. The upper limit 
is computed using Bayesian statistics with a flat prior for 
/+ between and 1 and yields /+ < 0.09 at 95% CL. Re- 
sults for individual samples differ by at most 1.8 standard 
deviations. 



Analyses based on cos 9* The angle, cos6'*, measured 
in the W boson rest-frame between the lepton and the b 
quark direction yields calculable distributions of this angle 
for each of the possible W boson helicities (—,0,-1-): 



do(r 



= - (1 - cos^ 
4 ^ 



, d±(r) = -(i±cosr). 

(61) 

Measuring cos 9* thus allows to reconstruct the contribu- 
tion of each of these helicities in top quark decays. 

CDF has used this observable in several analyses to 
measure the W helicity fractions in top decays |149lll50j . 
In [149) events with an isolated lepton, missing transverse 
energy and at least four jets including at least one identi- 
fied 5-jet of 1.9 fb~ are investigated. Two different meth- 
ods of reconstructing the full top pair kinematics are used 
to then compute the W boson rest-frame and cos 6** for 
each individual event. 

One method recovers the unmeasured neutrino mo- 
mentum from the missing transverse energy and from solv- 
ing the quadratic equation following from the W ^ iv 
decay kinematics when using the nominal W boson mass. 
This part of the analysis (for reasons that will become 
clear below) is called the "convolution" method. The other 
method (which will be called "template" method) uses 
a constrained kinematic fit to determine the lepton and 
parton momenta, where these momenta are allowed to 
float within the experimental uncertainties of the mea- 
sured quantities. Constraints are built from the W boson 
mass and the equality of the top and anti-top quark masses 
constructed from the fitted lepton and parton momenta. 



32 



Daniel Wicke: Properties of the Top Quark 



I 150 
100- 



50 





data 

f„ = 0.59 (f^ = 0.00) 
f^ = -0.04 (f„ = 0.70) 
f„ = 0.65 , f, = -0.03 



^--S^ 



a. 



• unfolded data 
f^ = 0.01 (f„ = 0.70) 

■ ■ ■ longitudinal 

■ - ■ left handed 
right handed 



-z^. 



+^-7+ 



- -3=^ 



■4r- 



:^ 



Oi^_, — , — , — , — I — , — , — , — , — I — , — , — , — , — I — , — , — , — , — I o' "-"' •• ■ " ■ — I — ■ — ■ — ■ — ' — I — ■ — ■ — ■ — ■ — I — - 

-1 -0.5 0.5 1 -1 -0.5 0.5 1 

cos e* cos 9* 
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data measured for the "template" method compared to the template prediction with the best fit parameters. The right plot 
shows the de-convoluted data of the "convolution" method compared to theoretical curves for individual helicity types and the 
best fit combination. 



Both methods require an association of the measured jets 
to the partons of the top quark pair topology and use the 
quality of the constrained fit to select this assignment. 

The reconstructed cos d* distribution is used in two 
different likelihood fits to determine the helicity fractions 
/o and /+ either individually, fixing the other one to the 
Standard Model value, or simultaneously. Both methods 
require simulation to construct their log likelihood func- 
tion. Pythia and Herwig are used to simulate Standard 
Model top quark pair production. Modified versions of 
Herwig and Madevent+Pythia [151 are used to gen- 
erate samples with varied helicity fractions. The W^-|-jets 
background is simulated with Alpgen-|-Herwig and nor- 
malised to the amount of data before 6-tagging and after 
removing all other background and ii signal events. Multi- 
jet background is taken from a control sample. Additional 
minor backgrounds from diboson, Z-|-jets and single top 
quark production are taken from simulation normalised 
according to their theoretical cross-sections. 

In the "template" method the background expecta- 
tions are combined with signal templates for the three 
different helicity states. The helicity fractions are taken 
from an unbinned likelihood fit with proper correction 
for acceptance effects. For this method an additional cut 
on the scalar sum of all transverse energy was required, 
Ht > 250 GeV. 

The "convolution" method uses the signal simulation 
to derive acceptance functions which are then convoluted 
with the theoretically predicted number of events in each 
bin of cos 9* . The helicity fractions are then taken from 
a binned likelihood fit after subtracting the background 
estimation from data. 

Systematic uncertainties for both methods were deter- 
mined from pseudo datasets with templates modified ac- 
cording to the systematic effect under consideration. The 
jet energy scale uncertainty is among the dominating sys- 
tematic effects in both methods. Only for the /o in the 
two dimensional fits the initial/final state radiation uncer- 
tainties are more important. In general the "convolution" 
technique yields slightly larger systematic uncertainties. 



Data are compared to the estimated results for the optimal 
fit parameters in Fig. [26l For the "convolution" method 
(right) the data shown are de-convoluted and compared 
to the pure theory prediction. 

The results of the methods are combined using the 
BLUE method [55] with the statistical correlation be- 
tween the two methods determined in pseudo experiments 
and the systematic uncertainties considered completely 
correlated. The combined results yields an improvement 
of about 10% compared to the individual methods. The 
model independent 2d fit yields 

/o = 0.66 ± 0.16(stat) ± 0.05(syst) 

/+ = -0.03 ± 0.06(stat) ± 0.03(syst) (62) 

with a correlation of —82% between /o and /+. Upper 
limits on the positive helicity fraction, /+ are not set. 



In |150j CDF investigates the dilepton channel in 
4.8 fb~^ to determine the W boson helicity fractions from 
cos 9*. The selection requires two isolated leptons and at 
least two jets one of which is required to be identified as 
6-jet. 

To reconstruct cos 9* the neutrino momenta are de- 
duced by assuming the top quark pair decay kinematics 
with an intermediate W boson. Requiring that the mo- 
menta of the corresponding decay products are compat- 
ible with the W boson mass and a top quark mass of 
175 GeV and using the missing transverse energy allows 
to determine the neutrino momenta up to an 8-fold ambi- 
guity. To account for the experimental resolutions this de- 
termination is repeated with appropriately smeared input 
momenta. Of the eight (smeared) solutions the one that 
yields the smallest invariant tt mass is used. The distribu- 
tion of measured reconstructed cos 9* values is compared 
to templates with a binned likelihood function. Templates 
for the signal are obtained for various values of /o and /+ 
using a customised Herwig generator. In addition tem- 
plates for background processes of diboson, Drell-Yan and 
fake leptons are added to obtain the complete expectation. 
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Significant systematic uncertainties come from the jet 
energy scale and from the uncertainty of the modehing of 
initial and final state radiation. This special study is also 
suffering from limited template statistics. 

In dilepton events of 4.8 fb~ CDF determines the he- 
licity fractions to be |150l 



/o = 0.78 ± 0.20(stat) ± 0.06(syst) 
/+ = -0.11 ± O.lO(stat) ± 0.04(syst) 



(63) 



Additional improvements are expected from relaxing the 
requirement of identified 6-jets. 



Analysis using the Matrix Element Method In addition 
to the two analyses with explicit reconstruction of the top 
quark kinematics described above, CDF has performed 
an analysis that uses the Matrix Element technique to 
measure the longitudinal W boson helicity fraction, /o, 
in 1.9 fb~ of data [152,153!. The analysis selects events 
with one isolated lepton, large missing transverse energy 
and at least four energetic jets including at least one jet 
identified as 6-jet. 

For each selected event a likelihood, L(/o, f+\{j}), to 
observe the measured quantities, {j} (the lepton and jet 
momenta and the missing transverse energy) , is computed 
as function of the W helicity fractions. This likelihood 
consists of a signal term which depends on /o and /+ , and 
a background term which is independent of these: 

L^ihJ+\{J}) = 

ftopPtiilJ}; /o, /+) + (1 - /top) Pw+MiJ}) (64) 

The signal and background probabilities, Ptt and Pw+jcts, 
are computed by integrating the differential parton level 
signal and background cross-sections according to the 
leading order cross-sections for qq — > tt and W + 4jet pro- 
duction, respectively. These integrations account for the 
experimental resolutions. The product of likelihoods for 
all selected events is evaluated for discrete values of /o 
and /+. At each point the signal fraction, /top, is chosen 
to minimise the total likelihood. The result for the lon- 
gitudinal helicity fraction and its statistical uncertainty 
are taken from the minimum of the log likelihood and its 
change by 0.5 units. 

The method is validated and calibrated using simula- 
tion for top pair signal and various background processes. 
Samples of various nominal /o, /+ values were created by 
reweighting the top pair events according to the expected 
cos 9* distribution. Applying the above method to various 
pseudo datasets with various nominal /o, /+ values yields 
a calibration curve. The observed slope of less than one 
is explained by the incomplete description of signal and 
background each with only a single leading order matrix 
element. The final results are corrected with this calibra- 
tion curve. 

The largest systematic uncertainty for this method 
stems from the uncertainty in simulation. It is estimated 
by checking the difference between Pythia and Herwig. 



In addition, uncertainties due to initial and final state ra- 
diation, different PDF sets, the jet energy scale uncer- 
tainty and various other experimental effects are consid- 
ered. 

Using the described Matrix Element method for 
2.7 fb^^ and assuming nit = 175 GeV CDF finds 



/o = 0.88 ± O.ll(stat) ± 0.06(syst) 
/+ = -0.15 ± 0.07(stat) ± 0.06(syst) 



(65) 



with a correlation of —59%. A variation of the assumed 
top quark mass shifts the results for /o by 0.017 and for 
/+ by —0.010 per ±1 GeV of shift in the top quark mass 
from the central value il53i. 



D0 



The D0 collaboration has investigated a total of 5.3 fb^ 
to determine the W boson helicity fractions based on the 
reconstruction of the decay angle cos 9* in lepton plus jets 
and in dilepton events |1541I155) . In the lepton plus jets 
events the hadronic decay is utilised to measure |cos0*| 
which helps to measure the longitudinal fraction, /q. 

Events are selected by requiring an isolated lepton, 
missing transverse energy and at least 4 jets. No sec- 
ond lepton is allowed in the event. Dilepton events are 
selected with two isolated charged leptons with opposite 
charge, large missing transverse energy and at least two 
jets. Some additional cuts are applied to suppress Z ^ U 
events in the ee and fj,fi channels and to assure a mini- 
mal transverse energy in the efj, channel. All channels use 
a multivariate likelihood discriminant based on kinematic 
observables and the neural network 6-jet identification to 
improve the purity of the top quark pair signal. 

In lepton plus jets events the decay kinematics of the 
top quark pairs is reconstructed using a constrained fit 
that determines the momenta of top quark and the W 
boson decay products from the measured jets and lepton 
momenta as well as from the missing transverse energy. 
Only the four leading jets in pT are used. The fit requires 
the momenta of W boson decay products to be consistent 
with the nominal W boson mass and the momenta of the 
top quark decay products to yield a top quark mass of 
172.5 GeV. Of the 12 possible jet-parton assignments the 
one with the highest probability of being the correct one is 
chosen. This probability for each association is computed 
using the fit x^ as well as the output value of the neural 
network 6-tagger for the four jets and its consistency with 
the light or heavy quark assignment under consideration. 
From the chosen solution cos 9* is computed from the lep- 
tonic side and a second measurement of absolute value 
from the hadronic side, see Fig. [57] (top and middle). 

The kinematics of dilepton events can be solved as- 
suming the top quark mass with a fourfold ambiguity. In 
addition the two possible assignments of the two leading 
jets in pt to the b quarks are considered. For each of these 
solutions the decay angle cos6'* is computed. To explore 
the full phase space consistent with the measurements. 
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Fig. 27. W decay angle distributions as measured in 4.3 fb~^ of 
D0 Run lib data. The top plot shows the distribution obtained 
from the leptonic decay in the lepton+jet events, the middle 
plot the one obtained from the hadronic decay. The bottom 
plot shows the distribution obtained from dilepton events. In 
addition to the data shown with error bars all plots contain the 
best fit prediction shown as full line and the Standard Model 
as a dashed line histogram. The shaded area represents the 
background contribution [155] . 



the measured jet and charged lepton momenta are fluc- 
tuated according to the detector resolution and cos 61* is 
computed for each fluctuation. The average over all solu- 
tions and all fluctuations is computed for each jet to find 
two cos 9* values per event. The resulting distribution for 
Run lib data is shown in Fig. [27] (bottom) . 

The expected distribution for signal top pair events 
is simulated using Alpgen-|-Pythia with various V — A 
to V + A ratios. These samples are then reweighted to 
form samples corresponding to the three W boson helicity 
states. Important backgrounds in the selected event sam- 
ple are M^-l-jets and multijct events in the lepton plus jets 
channel, 11^14^ -|-jets and Z-|-jets in the dilepton samples. 
Backgrounds of a single weak boson with jets are simu- 
lated with Alpgen-I-Pythia, diboson samples are gener- 
ated with Pythia. The multijet contribution is estimated 
from data for each bin of the cos 0* distribution. 

From the data and the estimated signal and back- 
ground contributions a binned likelihood, £(/o,/+), is 



computed for the observed data to be consistent with the 
sum of the backgrounds and the estimates for the three 14^ 
boson helicity states. The normalisation of the background 
is kept as a nuisance parameter with a Gaussian constraint 
to its nominal value. The measured helicity fractions, /o 
and /+, are those that minimise this likelihood. 

Systematic uncertainties are evaluated in ensemble 
tests. Pseudo datasets are drawn from models with sys- 
tematic variations and compared to the standard tem- 
plates to find the resulting shift in the obtained helicity 
fractions. Dominating uncertainties stem from the model- 
ing of ti, which are determined by exchanging the stan- 
dard Alpgen-|-Pythia with pure Pythia, Herwig and 
Mc@NLO. Further significant contributions to the system- 
atic uncertainties stem from background modeling, from 
the limited template statistics, from jet energy scale, jet 
energy calibration and jet reconstruction. An uncertainty 
due to the top quark mass uncertainty of 1.4 GeV is also 
included. 

Combining the datasets of Runlla and Run lib, which 
are analysed separately, D0 finds 



/o = 0.669 ± 0.078(stat) ± 0.065(syst) 
/+ = 0.023 ± 0.041(stat) ± 0.034(syst) 



(66) 



The correlation between the two numbers is about —0.8. 
The result shows a good consistency between the Run Ila 
and Run lib datasets as well as between the dilepton and 
lepton plus jet channel. As for the other good agreement 
with the Standard Model expectation. 



3.2 The CKM Element Vtb 

Another aspect of the weak coupling is that of flavour 
changing charged currents. The Standard Model explains 
these through the CKM matrix that needs to be deter- 
mined from experiment. The elements Vtd and Vts of this 
matrix have been determined from experiment assuming 
the Standard Model and allow to infer 0.9990 < \Vtb\ < 
0.09992 ^6iil56i using unitarity of a 3 x 3 CKM matrix. 
Physics beyond the Standard Model may invalidate these 
assumptions and leave | Vtb \ unconstrained [156^ . The value 
of I Vtb I directly influences the single top quark production 
cross-section and the ratio between top quark decays to 
6-quarks {t — >• bW) and to light-quarks {t — >■ qW with 
q — d,s). Both effects have been studied by the Tevatron 
experiments to constrain the CKM elements related to the 
top quark. 



3.2.1 Single Top Quark 

The single top quark cross-section from Standard Model 
sources is proportional to |Vtb| . Both experiments have 
used this relation to convert their cross-section measure- 
ments to determinations of the CKM element |1571I158| 
I159J . In addition to the uncertainty of the single top 
quark cross-section measurement, theoretical uncertain- 
ties on the Standard Model cross-section need to be taken 
into account. 
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The combination of the CDF and D0 measurements 
yields a cross-section for single top quark production of 



2-76lQ;47pb [ 159; assuming a top quark mass of 170 GcV 
in modelling the signal efficiencies. To extract the CKM 
matrix element it is assumed that | Vtb \ is much larger than 
I Vtd I and I Vts I , so that the top quark decay is dominated by 
the decay to Wb and no significant production through d 
and s quarks in the initial state takes place. No assumption 
about the unitarity of the CKM matrix is made. 

With theses assumptions the analyses performed to 
determine the cross-section can remain unchanged for the 
determination of \Vtb\. Attributing the full deviations of 
the experimental result from the SM prediction to the 
value of I Vtb I the combined CDF and DO single top quark 
production yields |159| : 

|V"th| =0.88±0.07 or \Vtb\ > 0.77 at 95%C.L. (67) 



3.2.2 Top Quark Pairs 

In top quark pair decays the number of identified b jets is 
used to measure the branching fraction for t — >■ bW. This 
fraction can be expressed in terms of the CKM matrix 
elements 



Rb 



\Vtb\' 



(68) 



\Vtd\ +\Vtsr + \Vtbr 

assuming that the top quark decay is restricted to Stan- 
dard Model quarks. 

CDF CDF has investigated 160 pb^^ of data using both 



events with lepton plus jets and dilepton events |160] . For 
the lepton plus jets sample events are selected by requir- 
ing an isolated lepton, missing transverse momentum and 
four jets. Dilepton events consist of two charged leptons, 
missing transverse momentum and two jets. Both sam- 
ples are classified according to the number of jets that are 
identified as b jets. 

The background in the lepton plus jets sample is dom- 
inated by W-|-jets events and multijet events with fake 
electrons. The multijet background is estimated from data 
using a range of control samples. The M^-|-jets background 
is simulated using the Alpgen-|-Herwig generator. In the 
subsamples with one or more than one identified &-jet, it 
is normalised using data before ^-tagging. The fraction of 
heavy flavour in these samples is scaled according to the 
Monte Carlo to data ratio in a control sample of inclusive 
jet events. For the description of VF-|-jets events without 
identified b jets, a neural network based on kinematic ob- 
servables is used, which enriches Vt^-|-jets background at 
low values and top quark pair signal at high values. The 
neural network was trained on Alpgen-|-Herwig sam- 
ples. The distribution of neural network output values 
measured in data is compared to the simulation of lF-|-jets 
and tt to fit the signal and background contribution. The 
shape of the multijet background is included at the rate 
determined above. 

For the dilepton sample the main backgrounds stem 
from Drell-Yan, diboson and from M^-t-jets events with 



fake leptons. The Drell-Yan background for ee and fifi is 
simulated using Pythia normalised to the number of Z 
bosons in a mass window around Mz- Other electroweak 
backgrounds are fully taken from simulation. The VF-|-jets 
backgrounds are taken from the Alpgen-|-Herwig simu- 
lation applying lepton fake rates, that are determined in 
a complementary jet sample. A tag rate probability for 
generic QCD jets is used to find the contribution of fake 
lepton events to the various &-tag subsamples. 

The distribution of the number of 6-tags for top quark 
pair events depends on the branching fraction, Rb. It is de- 
termined from events generated with Pythia and passed 
through full CDF detector simulation (as all simulations 
above) . 

Finally, a Poisson likelihood for the observed data to 
agree with the expectation is constructed as function of 
Rb. Gaussian functions with nuisance parameters are used 
to take systematic uncertainties including correlations be- 
tween the samples and the 6-tag bins into account. The 
dominant uncertainty comes from the background esti- 
mate in the 0-tag samples and the b quark identification 
efficiency. 

In the analysed 160 fb~^ CDF finds Rb = 1.12 ± 
0.2(stat)to:l3(syst) ^6U\. The Feldman- Cousins ap- 
proach |161l l5] is used to compute a lower limit of 



Rb > 0.61 or \Vtb\> 0.78 at 95% C.L. 



(69) 



where the conversion to the limit on the CKM element is 
done assuming three generations and unitary of the CKM 
matrix, only. 



D0 D0 has measured the top quark branching fraction, 
Rb, in conjunction with the top quark pair cross-section 
using 0.9 fb~^ J62. Events are selected for the lepton plus 
jets channel requiring an isolated lepton, missing trans- 
verse momentum and at least three jets. In data b jets are 
identified using a neural network tagger. 

Top quark pair signal is simulated with Pythia includ- 
ing samples in which one or both top quarks decay to a 
light quark and a W boson. The dominating VF-|-jets back- 
ground is simulated using Alpgen-I-Pythia. Its heavy 
flavour content of the M^-|-jets background was corrected 
according to a measurement in a control sample. The fake 
lepton background from multijet events is fully estimated 
from data. Additional smaller backgrounds from dibo- 
son, single top and .Z-l-jets are simulated using Pythia, 
SingleTop and Alpgen-|-Pythia, respectively, and nor- 
malised to their NLO cross-sections. All simulations are 
passed through the D0 detector simulation and recon- 
struction. In the simulation tag rate functions, determined 
on control samples in data, are used to describe the prob- 
ability for a given jet to be identified as b jet. Figure E51 
(left) illustrates the probability to have zero, one or more 
identified b jets in top quark pair events as function of the 
top quark branching fraction obtained from simulation. 

The event samples are separated by lepton type, num- 
ber of jets (3 or > 4) and number of identified b jets (0, 1 
or > 2). The 0-&-tag sample with four or more jets is fur- 
ther split in bins of a topological likelihood discriminant to 
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Fig. 28. Left: Probability to have zero, one or more identified b jets in top quark pair events as function of the top quark 
branching fraction, Rt- Middle: Data compared to contribution from various backgrounds as function of the topological likelihood 
discriminant. Right: Observed number of events as function of the number of identified b jets compared to expectations for various 
values of Rb |162| . 



obtain additional separation between W^+jets background 
and top quark pair signal events, c.f. Fig. E51 (middle) . 

To simultaneously determine the top quark pair pro- 
duction cross-section, ati, and the branching fraction, Ri,, 
a binned likelihood is constructed. Poisson distributions 
according to the expected event count as function of a^ 
and Rb is used for each sample and discriminant bin. The 
normalisation for VF-l-jets is fixed globally by subtract- 
ing all other backgrounds as well as top quark pair es- 
timates from data. Systematic uncertainties are included 
using nuisance parameters with Gaussian constraints. 

For the determination of Ri, the measurement is domi- 
nated by the statistical uncertainty. Systematic uncertain- 
ties are dominated by the uncertainty of the b tagging ef- 
ficiency. In contrast to the CDF measurement, due to the 
global determination of the M^-|-jets normalisation, the un- 
certainty on the size of this background is no significant 
source of uncertainty. 

In 0.9 fb~^ of lepton plus jets data D0 obtains Rt — 
0.97lQ'o8(total) (162j consistent with the expectation of 
the Standard Model. In Fig. [55] (right) the observed num- 
ber of events is compared to expectations for various val- 
ues of Rb- Limits on Rh obtained using the Feldman- 
Cousins procedure yield 

Rb > 0.79 at 95% C.L. (70) 

This limit is converted to a limit on the ratio of |Vtb| to 
the off-diagonal elements: 



Wt 



tb\ 



\Vtd\' + \Vtsf 



> 3.8 at 95% C.L. 



(71) 



The only assumption entering this limit is that top quarks 
cannot decay to quarks other than the known Standard 
Model quarks. Thus it is valid even in presence of an ad- 
ditional generation of quarks as long as the b' quark is 
heavy enough. 



3.3 Flavour Changing Neutral Currents 

Flavour changing neutral currents (FCNC) do not appear 
in the SM at tree level and are suppressed in quantum 
loops [163lll641ll65lll66j . However, anomalous couplings 
could lead to enhancements of FCNC in the top quark 
sector and their observation would be a clear sign of new 
physics [THTlfTHH] . 

The Tevatron experiments have looked for FCNC both 
in top quark decays |169j and in the production of (single) 
top quarks 170, 17T1 |172| . Limits on the single top produc- 
tion through anomalous couplings were also set with LEP 
and HERA data [T751[T71[T751[T751[T771[T751[T7gi[T5n] . 



3.3.1 Top Quark Decay through Z Bosons 

In an investigation of data with a total luminosity of 
1.9 fb^^ CDF looks for top quark pairs that show a flavour 
changing neutral current decay through a Z boson [169j . 
The analysis aims to identify events in which the Z bo- 
son decays leptonically and the second top quark decays 
through a W boson into hadrons. The event selection thus 
looks for a pair of leptons and at least four jets. The 
leptons need to be of the same flavour and have oppo- 
site charge. Their invariant mass is required to be within 
15 GeV of the Z boson mass. The cuts on the total trans- 
verse mass and the transverse energy of the leading and 
sub-leading jets were optimised in simulation. Events fail- 
ing these cuts are used as control sample, events passing 
these cuts are split into events without any identified b jet 
and events with at least one identified b jet. 

The dominant background with this selection stems 
from Standard Model Z-|-jets production, which is sim- 
ulated using Alpgen. Further but much smaller back- 
ground contributions stem from Standard Model top 
quark pair production, and diboson production. The sig- 
nal of top pairs with FCNC decay was simulated using 
Pythia. The events are reweighted to yield helicities of 
65% longitudinal and 35% left-handed Z bosons. 
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To separate signal from background the mass of the 
W boson is reconstructed from two jets, the top quark 
mass is reconstructed by adding a third jet and a second 
top quark mass is reconstructed from the Z boson with 
the fourth jet. A x^ variable is built from the differences 
of the reconstructed masses to the nominal W boson and 
top quark masses, respectively. The x^ of the jet parton 
assignment that yields the lowest x^ is used to build a 
distribution of x^ values. 

The estimated shapes of the various backgrounds and 
the signal events are used as templates that are fitted to 
the distribution measured in data. The main parameters 
of the fit are the branching fraction, S(i — >■ Zq), and the 
normalisation of the dominating Z+jets background (in 
the control sample). Further parameters describe the dif- 
ference of the background normalisation between the sig- 
nal and the control samples (with a Gaussian constraint), 
the b quark identification fraction and the jet energy scale 
shift. The latter is considered to cover all shape changing 
affects. 

The distribution of observed x^ values is shown with 
the best fit of the signal and background templates in 
Fig. [ini Data agree well with the Standard Model tem- 
plates and thus limits on the branching fraction t — >■ Zq 
are set. For this the Feldman-Cousins method is ap- 
plied and yields B{t -^ Zq) < 3.7% at 95% C.L. [1691 The 
Run I result in addition sets a limit on flavour chang- 
ing neutral currents in the photon plus jet mode of 
B{t -^ 7u) + B{t -^ 7c) < 3.2% [IM] . 



3.3.2 Anomalous Single Top Quark Production 

While the above study of top quark decays addresses a 
flavour changing neutral current through the Z boson, in- 
vestigations of the production of single top quark events 
can be used to restrict anomalous gluon couplings. 



CDF In an analysis of 2.2 fb"^ CDF looks for the produc- 
tion of single top quarks without additional jets, u{c) + 
g -^ t. To select such events with a leptonic top quark 
decay, one isolated lepton, transverse missing energy and 
exactly one hadronic jet are required. The jet must be 
identified as 6 jet. Additional cuts are used to reduce the 
backgrounds without a W boson as in the single top quark 



analyses |170j . 

To describe the expected background from Standard 
Model processes diboson and top quark pair events are 
simulated with Pythia and normalised to the NLO 
cross-sections. Single top quark events are simulated us- 
ing MadGraph-I-Pythia 182,183 . Finally, processes of 
weak vector bosons are simulated with Alpgen-|-Pythia. 
In these samples the heavy flavour contribution is en- 
hanced according to the findings in a control sample. The 
total normalisation of the W^-l-jets samples is taken from 
sideband data. The signal of FCNC production of single 



top quark is simulated using TopRex-|-Pythia [184) . 

Due to the large background from W^-l-ljet data, a neu- 
ral network is employed to differentiate between FCNC 
and Standard Model production. Fourteen observables, 
which each allow a significance of more than 3(7 in dis- 
criminating signal and background, were chosen as inputs 
to the neural network. They utilise kinematical properties 
of the measured quantities and the reconstructed W boson 
as well as the output of a special flavour separation neural 
network. The neural network is trained on samples with 
equal amount of signal and background. It is then applied 
to the individual signal and background samples to obtain 
templates for all simulated physics processes considered, 
see Fig. [20] (left). 

To determine the possible contribution of FCNC single 
top quark production the background templates are added 
according to their expected contribution and a binned 
maximum likelihood fit is used to measure the contribu- 
tion due to FCNC production. Systematic uncertainties 
are parametrised in the likelihood function with Gaussian 
constraints. They are dominated by uncertainties on the 
cross-sections of the background samples normalised to 
NLO and the selection efficiency for signal events. 

CDF finds no significant contribution of FCNC sin- 
gle top quark production in 2.2 fb^ of data. The limit 
on the allowed production cross-section crf^^^ is set us- 
ing Bayesian statistics with a flat prior for positive cross- 
sections and yields af'^^'^ < 1.8 pb at 95% C.L. [T^ . 
This cross-section limit is converted to limits on FCNC 
top quark-gluon coupling constants following J186II187) . 
Assuming that only one of the couplings differs from 
the Standard Model expectation CDF finds Kgtu/A < 
0.018 TcV"^ or Kgtc/A < 0.069 TeV~^ Expressed in terms 
of the top quark branching fraction through this processes 
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Fig. 30. Neural network output compared to 2.2 fb ^ of CDF data (left). Upper limits on the anomalous branching fractions 
derived from the limits on the anomalous couplings to u-quarks (middle) and c-quarks (right) |170II185] . 



these limits correspond to B{t -^ u + g) < 3.9 • lO^"' and 
B{t ^ c + g) < 5.7 • 10"^ as shown in Fig. [301 (middle 
and right). These small limits justify the approximation 
of pure Standard Model decays made in simulating signal 
samples above. 



D0 D0 has set limits on the FCNC anomalous couplings 
of the top quark in up to 2.3 fb"^ of data [TTTlfrT^] . The 
analyses investigate the singly production of a top quark in 
association with at least one additional jet. The event se- 
lection requires an isolated charged lepton, missing trans- 
verse momentum and at least two jets. Exactly one of the 
jets has to be identified as b jet. With this the selection 
follows closely the selection used for the measurements of 
the single top quark production 188,189,190,157 . 

In the most recent analysis the single top quark sam- 
ples for the SM and the FCNC signal processes are sim- 
ulated with the SingleTop generator 191,192J. Back- 
ground contributions from top quark pair, Vl^-|-jets and 
Z-|-jets production are simulated using Alpgen-|-Pythia 
and diboson production is simulated by Pythia. All sam- 
ples are passed through Geant to simulate the D0 detec- 
tor and then reconstructed with the standard event recon- 
struction. The SM samples of single top quarks, top quark 
pairs, Z-|-jets and dibosons are normalised to the NLO 
(or better) cross-sections. Events from M^-|-jets produc- 
tion are normalised to data before 6-tagging accounting 
for other simulated backgrounds and multijet background. 
The multijet background is described using data in which 
the lepton candidates fail one of the lepton identification 
cuts. 

To separate FCNC from Standard Model prediction 
the method of a Bayesian neural network (BNN) is em- 
ployed. A large number of variables is used as input to 
the neural network. These variables describe object and 
event kinematics, top quark reconstruction, jet width and 
angular correlations. For training the BNN the two FCNC 
processes including tgu and tgc couplings, respectively, are 
combined into a single training sample. Separate BNNs are 
trained for each lepton flavour and jet multiplicity. The 
distribution of BNN outputs observed in data agrees well 



with the pure SM expectation, thus limits on the allowed 
anomalous gluon couplings are computed. 

These limits are computed using a Bayesian approach. 
A likelihood for the distribution of neural network outputs 
observed in data to occur is computed from the events ex- 
pected in the Standard Model and the FCNC production 
of single top quark as function of the anomalous gluon 
couplings Kgtu/yl and Kgtc/A. The likelihood for each bin 
is based on a Poisson distribution. Systematics are taken 
into account by smearing the Poisson parameters with a 
corresponding Gaussian distribution. The dominant un- 
certainties stem from shape changing effects like those 
from the jet energy scale and the modelling of b quark 
identification. In addition normalisation uncertainties for 
the background simulations and the overall luminosity un- 
certainty give a significant contribution. 

The likelihood is folded with a prior flat in the FCNC 
cross-sections and exclusion contours are computed as 
contours of equal probability that contain 95% of the 
volume. The two dimensional limits on the squared cou- 
plings observed by D0 in 2.3 fb~^ [172] are shown in 
Fig. 1311 One dimensional limits are obtained by integral- 
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ing over one of the two anomalous couplings and yield 
ngtujA < 0.013 TeV"^ and Kgtc/A < 0.057 TeV"\ corre- 
sponding to branching fractions of ;B(t — T' u+g) < 2.0-10^^ 
and B{t^ c + g) < 3.9 • IQ-^. These D0 results are the 
currently most stringent branching fraction limits. 



3.4 Top Quark Charge 

The top quark's electrical properties should be fixed by its 
charge. However, in reconstructing top quarks the charges 
of the selected objects are usually not checked. Thus an ex- 
otic charge value oi\qt\ = 4e/3 is not excluded by standard 
analyses. Furthermore the amount of photon radiation off 
the top quark may a priori differ from the expection based 
on the top quark charge. The Tevatron experiments have 
searched for deviations from the SM expectation in both 
aspects of the top quark charge. 



3.4.1 Exotic Top Quark Charge 

To distinguish between the Standard Model and the exotic 
top quark charge it is necessary to reconstruct the charges 
of the top quark decay products, the W boson and the b 
quark. The W boson charge can be taken from the charge 
of the reconstructed lepton, but finding the charge of the 
b quark is more difficult. 



D0 D0 has performed an analysis of £-|-jets events 
with at least two 6-tagged jets in 370 pb~^ using a 
jet charge technique to determine the charge of the b 
jets |193) . Semileptonic events are selected following the 
cross-section analysis by requiring exactly one isolated lep- 
ton, transverse missing energy and four or more jets. At 
least two of the jets must be identified as b jets using a 
secondary vertex tagging algorithm. 

The charge of a jet can be defined as the sum of the 
charges of all tracks inside the cone of that jet. In this 
analysis the sum has been weighted with the component 
of the track momenta transverse to the jet momentum, 
PX- 



iict 



T.P1, 



with 



0.6 



where the sums run over all tracks, i, within the jet under 
consideration and qi is the charge sign of the track i. Be- 
cause particles may easily escape the jet cone such a jet 
charge fluctuates strongly from event to event, so only sta- 
tistical statements can be made. It is crucial to determine 
the expected distribution of Qjot in the case of 6 or 6 quark 
and, because a significant fraction of charm quarks gets 
flagged by the secondary vertex tagger, also for the c and 
c quarks. These expected distributions, c.f. Fig. |31] (left), 
are derived from dijet data using a tag and probe method. 
To determine the top quark charge an assignment of 
6-jets to the leptonic or hadronic event side is necessary. 
This analysis uses the quality of a fit to the ti hypothesis. 
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which uses the W boson and top quark masses as con- 
straints, to select the best possible assignment. The jet 
charge for the b jet on the leptonic (hadronic) side, gb, 
ilbh) is then combined with the charge of the measured 
lepton qi to define two top quark charge values per event: 
Qicp = \qi+qb^ and Qhad = \-qi+qbJ\- The distribu- 
tion of the measured top quark charges is compared to 
templates simulated for the Standard Model and the ex- 
otic case, where the exotic case has been obtained by in- 
verting the jet charge, see Fig. [32] (right). The top quark 
pair events are simulated using Alpgen-|-Pythia. The 
dominating W^-|-jets background is simulated using Alp- 
GEN-I-Pythia with a normalisation to data. Multijet tem- 
plates are derived from data alone. All simulated events 
are passed through full D0 detector simulation. 

An unbinned likelihood ratio accounting also for re- 
maining background yields a p-value for the exotic case 
of 7.8% and a Bayes factor of Bf = 4.3 favouring the 
Standard Model charge scenario |1931I194) . 



CDF, Jet Charge The corresponding analysis by CDF in- 
vestigates 1.5 fb^ with events from the semileptonic and 
the dileptonic decay channel |T95 . The former are selected 
requiring an isolated charged lepton, missing transverse 
momentum and at least four jets. Two of the jets are re- 
quired to be identified b jets using CDFs secondary vertex 
algorithm. Dilepton events are selected by asking two op- 
positely charged leptons, missing transverse momentum 
and at least two jets, one of which needs be identified as 
& jet. 

Compared to the D0 analysis the jet charge is com- 
puted slightly differently. Instead of the transverse mo- 
mentum the scalar product of the jet and the track mo- 
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mentum is used to weigh the measured charges: 



Q 



jet 



E (Pi ■ Pjct) 



with K ~ 0.5. 



(73) 



Depending on the sign of Qjct the identified b jet is con- 
sidered to stem from the b or the b quark. The purity 
of this assignment is caUbrated on dijet events with two 
identified b jets. One of the jets is required to contain a 
muon, that serves as tag in the tag and probe method. 
The resuhing purity is corrected for effects due to muons 
from secondary decays, for B meson mixing and for hght 
or c-quark jets misidentified as b jets. 

To compute the charge of the top and anti-top quark 
the jets need to be associated to the leptons. In the 
semileptonic channel a kinematic fit with constraints on 
the top quark mass and the VF-boson mass is used. The 
jet-parton association with the lowest x^ of this fit is kept. 
In the dilepton channel the invariant mass of each pair of 
one lepton and one jet, M;|, is computed. The combina- 
tion which does not produce the largest value of M^f^ is 
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used. In both channels cuts on x^ and Af;|, respectively, 
are used to enhance the purity of correct assignments. 

Each event can now be classified as Standard Model 
like or as exotic model like. To obtain a statistical interpre- 
tation a likelihood is computed as function of the fraction 
of Standard Model like signal pairs, /+. Nuisance param- 
eters that represent the number and purity of signal and 
background events are optimised for each value of /+. The 
systematic uncertainties considered include effects from 
the choice of parton density function, the uncertainties in 
the simulation of initial and final state radiation, the jet 
energy scale and the choice of the generator. All system- 
atic uncertainties are included in the statistical treatment 
through their effect on the nuisance parameters. 

Background predictions are obtained as for other CDF 
lepton plus jets and dilepton analyses based on a mixture 
of simulation and data. For this analysis each background 
is checked for a correlation between the charge of the signal 
lepton and the jet charge value of the corresponding b jet. 
Such a correlation could occur for the semileptonic channel 
from the bb background when the lepton from the b decay 
passes the lepton selection criteria and from single top 
quark events. In both cases the correlation is found to be 
small and consistent with zero within uncertainties. 

In 1.5 fb^^ of data CDF finds the most likely value of 
the fraction of Standard Model like signal events to be 
/+ = 0.87. This corresponds to a p- value of 31% |195j . see 
also Fig. [321 Because this is larger than the a priori chosen 
limiting probability of 1% to falsely reject the Standard 
Model hypothesis, CDF claims to confirm the Standard 
Model hypotheses. The confidence limit corresponding to 
this 1% choice is computed as 87%. CDF computes the 
Bayes factor to be 21og(-B/) = 12, which shows that this 
analysis with 1.5 fb~ yields a much stronger exclusion of 
the exotic hypothesis than the D0 analysis of 370 pb^^ 
described above. 
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Fig. 33. Top: Distribution of the obtained fraction of Stan- 
dard Model like pairs, /+, in ensembles with the exotic and 
the Standard Model top quark charge. Indicated is the value 
measured in data: /+ = 0.87. Bottom: Distribution of p-values 
for the Standard Model hypothesis to be true obtained in en- 
sembles of the exotic model. Indicated is the a priori chosen 
hmiting probability of incorrectly rejecting the Standard Model 
if Standard Model is true, 1%, and the corresponding proba- 
bility of rejecting the Standard Model if the exotic model is 
true, 87% il95l. 



CDF, Soft Lepton Tag With more data an alternative 
method to determine the charge of the b quarks is to con- 
centrate on the leptonic decays of the b quarks. CDF ap- 
plies this method on 2.7 fb~^ |196) . 

Events are selected requiring an isolated energetic lep- 
ton, missing transverse energy and at least four jets. Ve- 
toes are applied against additional isolated energetic lep- 
tons, conversion electrons, cosmic muons, and Z bosons. 
In the four jets at least one soft lepton and one secondary 
vertex tag must be reconstructed. The soft lepton tag is 
optimised to select leptonic b quark decays, b — >■ li^X, but 
to suppress cascade decays, 6 — )• c — )• ivX. With this the 
charge of the b quark can be deduced from the measured 
soft lepton charge. 

To resolve the ambiguities of assigning the measured 
jets and energetic leptons to the (anti-)top quark decay 
products a kinematic fitter is applied. For the computa- 
tion of the top quark charge the solution with the best 
fit quality is used. Only solutions which assign the tagged 
jets to b quarks are considered. Events with a bad fit qual- 
ity are rejected, if x^ < 9 or x^ < 27 for events with one 
or two tagged jet, respectively. 
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To estimated the expected performance of the anal- 
ysis, signal events are generated with Pythia using 
EvtGen |197) for the decays. Background from M^+jet 
events are generated using Alpgen+Pythia correcting 
the heavy flavour contribution by a factor derived in 
W + Ijet events. These samples are normalised to the 
data before requiring any b tags. The contribution of di- 
boson single top, Z+jets and Drell-Yan events is consid- 
ered using simulation normalised to the theoretical or ex- 
perimentally measured cross-sections. Contributions from 
multijet events with fake leptons are estimated from data. 
The total background estimation is B = 2.4 ± 0.8 events. 
The purity of the charge determination is determined from 
these simulations and then calibrated in data using pure 
bb events. 

In 2.7 fb~^ of data CDF finds a total of 45 events, 
TVsM — 29 of which are reconstructed as SM like and 

A^XM 



16 with the exotic charge value [196) . The statis- 
tical significance of this result is obtained by considering 
a large number of pseudo experiments. Figure |3H shows 
the expected distribution in terms of the asymmetry 



A 



J_NsM__NxM±BDb 
Ds NsM + NxM-B 



(74) 



where Ds and Di, are the dilution factors obtained in the 
calibration step for signal and background, respectively. 
The SM scenario corresponds to positive values, the ex- 
otic model to negative values. In 69% of the pseudo ex- 
periments assuming the SM the asymmetry turned out 
to be smaller than the measured value. Only in 0.9% of 
the pseudo experiments assuming a fully exotic top quark 
showed an asymmetry larger than the measured value. 
From these studies the exotic charge model can be ex- 
cluded at 95%C.L. 




-2024 

Normalized Asymmetry 



Fig. 34. Distribution of the asymmetry observed in a large 
number of pseudo experiments for the SM and the exotic charge 
case. The measured result is indicated by a vertical line [196J . 



3.4.2 Photon Radiation of Top Quarks 

The top quark charge is also expected to define the amount 
of photon radiation off the top quark. The measurement of 
this process is thus a complementary way of verifying the 
electrical charge of the top quark. The CDF collaboration 
has searched for top quark pairs with an additional photon 
in 1.9 fb~ as part of a more general search (198) . 

For the search of the tt^ final state, events are selected 
that contain an isolated energetic photon, one isolated en- 
ergetic lepton, large missing transverse energy and at least 
three jets. One of the jets is required to be identified as 
&-jet. 

With this selection the dominant backgrounds stem 
from jets misidentified as photons and from light-jets 
misidentified as &-jets. The former is estimated from data 
by measuring the photon fake rate in jet events as func- 
tion of f^T- This fake rate is then applied on events pass- 
ing the selection without the photon requirement. Also the 
amount of fake &-jets is determined from data and then ap- 
plied to the events without the requirement of a 6-jet iden- 
tification. Additional backgrounds including W^7+jet and 
diboson events are estimated using MadGraph+Pythia 
normalised to NLO cross-sections. MadGraph-|-Pythia 
passed through the full detector simulation and the re- 
construction is also used to determine the signal selection 
efficiencies. 

In 1.9 fb^ of data CDF finds a total of 16 events which 
pass the described selection. The expectation including 
SM tij events is 11.2 ± 2.2 events. Attributing the fuU 
difference between the expectation and the data to the 
tt'j process, CDF measures the cross-section for radiation 
off top quark pair events of 



atij = 0.15±0.08pb 



(75) 



with the quoted uncertainty being dominated by the sta- 
tistical uncertainty. The result is in agreement with the 
SM expectation of 0.080 ± 0.011 pb [TM] . 



3.5 Spin Correlations 

At the level of the hard interaction the spins between the 
top and the antitop quark are correlated in top quark pair 
production. Because top quarks decay before they hadro- 
nise these correlations are conserved in the weak decay and 
thus "good" quantum-mechanic observables |199ll200ll201l 
120211203) . The degree of correlation depends on the pro- 
duction and decay processes. Its measurement thus probes 
the details of the production mechanism (when assuming a 
SM weak decay). The amount of correlation also depends 
on the reference axes used to define the top and anti-top 
quark spin states. At the Tevatron the spin-correlations 
expected in the SM are largest in the so-called off-diagonal 
and beam bases. 

Both Tevatron experiments have measured the spin 
correlation of top quark pairs in different channels and 
using different spin bases. 
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3.5.1 Dilepton Channel 

In the dilepton channel both experiments use the nor- 
malised double differential cross-section with respect to 
the angle of flight directions of the two leptons as their 
observable. The angle of flight, cos9±, is measured with 
respect to the spin quantisation axis chosen in the analy- 
sis. This double differential cross-section depends on the 
spin correlation coefficient k: 



dV 



1 



a dcos^+dcos0_ 



= -(l + Kcos6'+cos0_) (76) 



In the off-diagonal and beam bases the SM predicts a cor- 
relation coefficient k of about 0.8. 



CDF The CDF analysis of the dilepton channel is based 
on 2.8 fb~^ |204j . The event selection requires two leptons 
of opposite charge, large missing transverse energy and at 
least two energetic jets. A veto is applied on Drell-Yan 
events and on special multijet fake configurations. 

The analysis selects the off-diagonal basis for the defi- 
nition of the flight directions. In this basis the spin quan- 
tisation axis is defined in the tt centre-of-mass frame. It is 
defined as the direction which deviates from the direction 
of flight of the (anti-)top quark by the angle ^ in clock- 
wise direction. The value of ^ is computed as function of 
the top quark velocity, (3, and the angle between the top 
quark and proton flight direction, 9*: 



tan ^ == v^l -/32 tan ( 



(77) 



The directions of flight with respect to the direction de- 
fined by ^ are measured in the rest frames of the top and 
the anti-top quark. The analysis uses the angles obtained 
for the leptons (cos0+,cos0_) and the (anti-)5-quarks 
{cos 9b, cos 9^). 

To reconstruct these angles in an individual event a 
full kinematic reconstruction of the top quark pair and its 
decay to l^l^vvbb is necessary. In the dilepton channel 
with its two neutrinos this requires the six constraints. 
The reconstructed W boson and the top quark masses 
need to be consistent with their nominal values (two con- 
straints each). In addition the sum of the transverse mo- 
menta of the reconstructed neutrinos are required to agree 
with the measured missing transverse energy, ^t- Due to 
the quadratic nature of the corresponding equations these 
constraints yield up to four solutions for the unmeasured 
neutrino momenta, Pi, and pp. 

To further improve the reconstruction a kinematic fit 
is applied that varies the reconstructed (anti-)6 quark en- 
ergies within the experimental resolution of the b jet mea- 
surements and the reconstructed sum of transverse neu- 
trino momenta within the resolution of ^t- Besides these 
resolutions the kinematic fit includes the probability den- 
sities for the distribution of the Pz , the pt and the invari- 
ant mass of the top quark pair in its likelihood. Of the 
initially up to four solutions for the neutrino momenta, 
Pv and Pv, the one which gives the best likelihood is used 
to compute the flight directions in each event. 



Simulations are used to determine the expected out- 
come of this measurement as function of the spin correla- 
tion. For the tt signal simulation Pythia is used. Because 
the Pythia samples do not contain spin correlations be- 
tween the generated top quarks, the generated events are 
weighted proportional to 1 + k cos 0+ cos 9- according to 
the generated values of cos 6*+ and cos0_. Background of 
dibosons and Drell-Yan are simulated using Pythia and 
Alpgen-|-Pythia, respectively. Background due to fake 
leptons is simulated from data with a single energetic lep- 
ton and jets. In these events one jet is artiflcially inter- 
preted as a fake lepton. Both signal and background tem- 
plates are smoothened by a polynomial function, where 
for the signal template the dependence on k, is kept. 

The measured two-dimensional distributions of the re- 
constructed (cos0+,cos0_) and (cos 6*;,, cos 6*^) values are 
now compared to the templates to determine the spin cor- 
relation coefficient k, using a likelihood fit. The measured 
value is corrected slightly according to a calibration per- 
formed on ensembles of pseudo experiments. 

In 2.8 fb~ of dilepton top quark pair events CDF mea- 
sures the spin correlation coefficient in the off-diagonal 
basis to be 

K = 0.32l°:^^ . (78) 

The total uncertainty is by far dominated by the statistical 
error. The leading systematics stems from the uncertainty 
on relative contribution of signal and background events. 
At the available statistics this results poses only a very 
minor constraint on the real spin correlations. 

D0 The measurement of the spin correlation by D0 uses 
dilepton events of up to 4.2 fb~ of data [205 . Events are 
selected requiring two oppositely charged leptons, large 
missing transverse energy and at least two jets. Drell-Yan 
events are vetoed near the Z boson resonance. 

In this study the fiight directions in Eq. (j76p are com- 
puted using the beam basis. The spin quantisation axes 
for the beam basis are the directions of the proton and 
anti-proton in the tt rest frame. The flight directions with 
respect to the quantisation axes are measured in the top 
and anti-top quark rest frame, respectively. D0 uses the 
direction of fiight angles for the two measured leptons, 
cos 9+ and cos 0_ . 

For the full reconstruction of the top quark pair kine- 
matics D0 applies the Neutrino Weighting Method de- 
scribed in the context of top quark mass measurement 
in Sect. 12.3.11 For the spin correlation measurement the 
weight, W , is computed for a fixed top quark mass of 
175 GeV, but scanning possible values of the flight angles 
to find its dependence on the product cos0+cos0_, i.e. 
for each event W = W{cos9j^ cos 9-). For further analysis 
the mean, /i, of the weight function in each event is used: 



H = xW{x) dx 



(79) 



The distribution of /i values found in the selected 
events is now compared to templates obtained from simu- 
lation. Signal templates for a range of k values are ob- 
tained by reweighting Pythia. Templates for the Zf-f 
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events is generated with Alpgen+Pythia, while for di- 
boson processes Pythia is used alone. Both background 
types are normalised to the theoretical cross-sections. To 
enhance the bad simulation of the Z boson -pT distribu- 
tion, a reweighting of these events has been applied. The 
measured spin correlation parameter k is obtained by a 
binned likelihood fit of the templates to the measured 
data. 

D0 has studied the performance of this method in en- 
sembles of pseudo experiments using various nominal val- 
ues of K. This calibration is used in the Feldnian- Cousins 
procedure p,61,6j to obtain the final results. Using dilep- 
ton events D0 determines the spin correlation coefficient 
in the beam axis 



K - -n 1 7+0-64 



(80) 



The uncertainty is dominated by statistics. The leading 
systematic uncertainty is found to come from signal mod- 
elling. This includes a dependence on the assumed top 
quark mass as well as differences observed when replacing 
Pythia by Alpgen+Pythia or Mc@nlo. The observed 
spin correlation is in slight tension with the SM expecta- 
tion of about 0.8. 



by reweighting Pythia events according to the generated 
decay angles. In addition templates for the expected back- 
ground are produced, including multijet, M^-fjets and di- 
boson events. 

To determine the spin correlation the two dimensional 
distribution of cos 9i cos 6d vs. cos 6'^ cos 0b is fitted with 
the described templates using a binned likelihood fit. The 
contribution of the same spin and the opposite spin tem- 
plates are allowed to float freely, the background contribu- 
tion is allowed to float but constrained to the expectation 
within its errors. The determined fractional contribution 
of the opposite spin contribution, /o, is converted to the 
spin correlation coefficient using n — 2fo — 1. 

In 5.3 fb~ CDF determines the spin correlation in the 
helicity and in the beam bases to be 



Khclicity = 0.48 ± 0.48stat ± 0.22, 
Kbcam = 0.72 ± 0.64,tat ± 0.26, 



^syst 
*syst 



(81) 



The measurement is clearly dominated by the limited 
statistical uncertainty. Of the systematic uncertainties the 
uncertainty on the signal modelling is by far dominating. 



3.5.2 Lepton plus Jets Channel 

CDF has also determined top quark spin correlations 
in the semileptonic decay channel [20611207) . In 5.3 fb~ 
events with one energetic lepton, large missing transverse 
momentum and at least 4 jet are selected; at least one of 
the jets is required to be identified as 6-quark jet. 

This analysis uses the the beam and the helicity bases. 
The spin quantisation axes are defined in the ii rest frame. 
For the beam basis the beam direction is used, for the 
helicity basis the (anti-)top quark directions are taken. 
Then the angles between the fiight directions of the top 
quark decay products, bbli^qq' , and the quantisation axis 
are reconstructed in the (anti-)top quark rest frame. The 
analysis considers the lepton direction, cos 9i, the bottom 
quark from the hadronic top quark decay, cos 9b , and the 
down- like quark, cos 9d- 

For the determination of these angles the full event 
kinematics of the top quark pair decay is reconstructed 
using a kinematic fitter which constrains the top quark 
mass to 172.5 GeV. Only events with a reasonable fit qual- 
ity are considered in the analysis. In the following analy- 
sis the distributions of the cos of the reconstructed angles 
are compared to templates. For the reconstructed value of 
cos6'd the jet which is closest to the 6-jet in the W rest 
frame is used. This is found to be the down-type quark in 
about 60% of the cases. 

For the analysis of the angles in the helicity basis, sam- 
ples of the four possible helicity combinations of the top 
and anti-top quark are generated with a modified version 
of Herwig. From these four combinations, equal helicity 
and opposite helicity samples are obtained assuming par- 
ity and CP conservation. In the analysis of the angles with 
respect to the beam basis the signal samples are obtained 



3.6 Charge Forward Backward Asymmetry 

At the Tevatron the initial state of proton anti-proton is 
not an eigenstate under charge conjugation. Thus in prin- 
ciple also the final state may change under this operation. 
In QCD, however, such a charge asymmetry appears only 
at next-to-leading order and arises mainly from interfer- 
ence between contributions symmetric and anti-symmetric 
under the exchange of top and anti-top quarks |2081I209| 
i210u211U212j . 

Experimentally, CDF and D0 investigated forward 
backward asymmetries pTH^THI^TSlI^ 



A 



FB 



Nf-Nb 
Nf + Nb 



(82) 



where A^f and A'^b are the number of events observed in 
the forward and backward direction, respectively. The for- 
ward and backward directions are either defined in the 
laboratory frame, i.e. according to the sign of the rapidity 
of the top quark, yt , or can be defined in the frame where 
the top quark pair system rests along the beam axis, i.e. 
according to the sign of the rapidity difference between 
top and anti-top quark. Ay = yt — yt- The two different 
definitions of forward and backward yield two different 
asymmetries that are labelled A^q and Ap^ according to 
their rest frame of definition. In the Standard Model at 
NLO asymmetries are expected to be 0.05 and 0.08, re- 
spectively |217j . but at NNLO significant corrections are 
predicted for the contributions from ti + X |218j . 

The smallness of the asymmetries expected within the 
Standard Model make them a sensitive probe for new 
physics. 
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Fig. 35. Distribution of the top quark rapidity (left) and rapidity difference (middle) as measured by CDF in 5.3 fb ^ compared 
to the Standard Model prediction. Right: The functional dependence of the rapidity difference of the top quark pair invariant 
mass, Mtt EH]. 



CDF 

The CDF collaboration has investigated up to 5.3 fb~ of 
data and measures both charge asymmetries defined above 
from top quark pairs with semileptonic decay |214[|216] . 
The event selection requires an isolated lepton, missing 
transverse energy and at least four hadronic jets, one of 
which must be identified as b jet. 

The top and anti-top quark kinematics are recon- 
structed from the jet momenta, the lepton momentum and 
the missing transverse momentum using mass constraints 
from the W boson and the top quark. The reconstructed 
values of these masses are constrained by the nominal val- 
ues of Mw = 80.4 GeV and to* = 172.5 GeV and b tagged 



jets are assigned to b quarks only [216j . Of the possible jet 
parton assignments the one with the best fit probability 
is taken. The rapidity of the hadronically decayed (anti-) 
top quark, yh, is multiplied by minus the charge, Qe, of 
the lepton to obtain the top quark rapidity: yt = —QeVh- 
For the tt frame asymmetry, A^^ the rapidity difference is 
computed Ay ^ yt- yt as Ay = Qi {yi -yh), with Qi and 
yh as above and yi being the rapidity of the leptonically 
decayed (anti-)top quark. 

The NLO Standard Model expectation of top quark 
pair production is done with MCFM [219] and the next- 
to-leading order generator Mc@NLO |220| which contain a 
small asymmetry. Leading order signal simulation without 
asymmetry from Pythia is used to check for any detec- 
tor or selection asymmetry. The dominating background 
events of VF-l-jets are simulated with Alpgen-|-Pythia, 
diboson backgrounds and single top quark events are sim- 
ulated with Pythia and Madevent, respectively. The 
normalisation of the VF-l-jets background and contribu- 
tions from misreconstructed multijet events are estimated 
from data. 

The uncorrected rapidity and rapidity difference dis- 
tributions measured in data are compared to the expecta- 
tions in Fig. [33] These distributions differ from the true 
particle level shape due to acceptance and reconstruc- 
tion effects. After background subtraction CDF derives 
the particle level distributions inverting the acceptance 
efficiencies and migration probability matrices as derived 
from Pythia simulation with zero asymmetry using a re- 



duced number of only four bins. The final asymmetries are 
computed from these unfolded distributions. 

The dominating systematic uncertainties are the back- 
ground normalisation and shape. For v4pg the amount of 
initial and final state radiation contributes significantly, 
while for Ap^ the jet energy scale is the next leading un- 
certainty. Further uncertainties from the parton distribu- 
tion functions, due to colour reconnection and the MC 
generator are considered. 

The final asymmetries measured in the update to 
5.3 fb"^ of CDF data are WB: 



(83) 



A% = 0.150 ± 0.050,tat ± 0.024,yst 
^^B = 0.158 ± 0.072stat ± 0.017,yst 



These values are somewhat larger than the 0.038 and 0.058 
expected in the Standard Model at NLO, respectively, but 
agree within two standard deviations. 

CDF also investigated the dependence of the asymme- 
try on several other topological and kinematic properties. 
Considering the two ranges of Ay yields: 

A*^Q{\Ay\ < 1.0) = 0.03 ± 0.12 (SM: 0.39) 
A*^Q{\Ay\>l.0)^0M±0.26 (SM: 0.123). (84) 

Clearly, the deviation from the expectation is driven by 
the effect at large rapidity differences. 

The functional dependence of the asymmetry on 
the invariant mass of the top quark pairs is shown in 
Fig. [551 (right). Separating this result in two bins yields 



^FB(^tt<450GeV) 
^FB(^tt>450GeV) 



-0.12 ±0.15 (SM: 0.04) 
-0.48 ±0.11 (SM: 0.09). (85) 



The events with high invariant mass show a deviation 3.4 
standard deviations from the NLO prediction obtained 
with MCFM. 

CDF completes their study by verifying that the asym- 
metries are consistent with CP conservation by separately 
considering events with positively and negatively charged 
leptons. The slightly enhanced asymmetry observed in the 
inclusive measurement thus seems to stem from effects at 
large rapidity difference and high invariant mass of the 
top quark pair system. 
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D0 

The D0 collaboration has investigated up to 4.3 fb'^ 
of data to measure Apg in semileptonic top quark pair 
events |2131I215) . The event selection requires exactly one 
isolated lepton, missing transverse momentum and at least 
four jets, the hardest of which must have pr > 40 GeV. 
At least one of the jets is required to be identified as b jet 
with D0's neural network tagger. 

The top quark pair kinematics is reconstructed by fit- 
ting the momenta of the top quark decay products to 
the measured jet and lepton momenta and the missing 
transverse energy with constraints on the reconstructed 
W boson and top quark mass to 80.4 GeV and 172.5 GeV, 
respectively. Only the b jets and the remaining three lead- 
ing jets are used. The possible jet parton assignments are 
reduced by assigning identified b jets only to b quarks. In 
the final analysis only the assignment with the best fit 
probability is used. The rapidity difference with the cor- 
rect sign is determined from the rapidities reconstructed 
of the leptonic and the hadronic side (y£ and yh) and the 
lepton charge, Qf. Ay = Qg {ye - yh). 

To estimate the dominant background of VF-|-jets pro- 
duction a set of observables well described by the simu- 
lation is used to construct a likelihood discriminant that 
does not depend on Ay. The expected shape of top quark 
pair signal and the T^-t-jets background in the distribution 
of the discriminant and on the asymmetry is determined 
from Mc@NLO and ALPGEN-f Pythia simulation, respec- 
tively, passed through D0 detector simulation and recon- 
struction. The effect of multijets events on the asymmetry 
and the discriminant is determined from data that fail the 
lepton identification. Other backgrounds were checked to 
have negligible effects. 

The final reconstructed asymmetry, Apg, in signal 
events is determined by maximising the combined like- 
lihood of the observed discriminant distribution and the 
distribution of the sign of Ay as function of the signal and 
background contributions and of the signal asymmetry. 
The dominant systematic uncertainties on the asymmetry 
are the jet energy calibration and the asymmetry recon- 
structed in VF-|-jets events. All of them are much smaller 
than the statistical uncertainty. 

In 4.3 fb~^ of data D0 finds a final observed asymme- 
try Em of 



A 



tt obs 
FB 



= 0.08 ± 0.04,tat ± O.OI 



syst 



(86) 



To keep the result model independent and in contrast to 
the CDF results this number is not corrected for accep- 
tance and resolution effects. Instead it needs to be com- 
pared to a theory prediction for the phase space region 
accepted in this analysis which is corrected for dilution 
effects. For NLO QCD and the cuts used in this analysis 
D0 evaluates A% = O.Oll[J;^i. Thus as for CDF this re- 
sult corresponds to an asymmetry that is slightly higher 
than expected in NLO QCD, but not by more than two 
standard deviations. 

In addition to the QCD expectation in the published 
analysis [213) D0 provides a parameterised procedure to 
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Fig. 36. Limits on a possible fraction, /, of resonant top quark 
pair production through a Z' boson obtained from the mea- 
surement of the forward backward asymmetry in D0 [213] . 



compute the asymmetry expected for an arbitrary model 
of new physics. As an example the measurement's sen- 
sitivity to top quark pair production via a heavy neutral 
boson, Z' , with couplings proportional to that of the Stan- 
dard Model Z boson is studied. Pythia is used to obtain 
a prediction of this kind of top quark pair production and 
due to the parity violating decay yields large observable 
asymmetries of 13 to 35% depending on the assumed Z' 
boson mass. Limits on the possible fraction of heavy Z' 
production are determined as function of the Z' boson 
mass using the Feldman- Cousins approach. These limits 
are shown in Fig. [36] and can be applied to wide Z' reso- 
nance by averaging the appropriate mass range. 



3.7 Differential Cross-Section 

Measurements of the differential cross-sections of top 
quark pair production can be used to verify the produc- 
tion mechanism assumed in the Standard Model. Due to 
the required unfolding these measurements are especially 
cumbersome. The CDF collaboration has measured the 
differential cross-section with respect to the invariant top 
quark pair mass, ^jffziMu), using 2.7fb"^ of data P^ . 
The event selection requests a lepton with high transverse 
momentum, large missing transverse momentum and at 
least four jets. At least one of the jets needs to be identi- 
fied as &-jet. 

The invariant mass of the top quark pairs is recon- 
structed from the four-momenta of the four leading jets 
in pt, the four momentum of the lepton and the missing 
transverse energy. The z-component of the neutrino is not 
reconstructed but used as if it was zero |222| . 

The dominating background in this selection stems 
from VK-f jets production. Its kinematics simulated with 
Alpgen-|-Pythia correcting heavy flavour contribution 
for differences between data and Monte Carlo. The re- 
quired normalisation is measured in data before applying 
the &-jet requirement [223 . Multijet background is extrap- 
olated from data with low missing transverse momentum. 
The smaller backgrounds of diboson, Z-|-jets and single 
top quark is fully taken from simulation using Pythia, 
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Fig. 37. Differential top quark pair production cross-section 
measured by CDF in 2.7 fb~^ of data using the semileptonic 
decay mode. Indicated are the total uncertainties for each bin, 
excluding the overall luminosity uncertainty of 6% [221| . 



Alpgen+Pythia and MadGraph, respectively. All sim- 
ulated events are passed through the CDF detector simu- 
lation and reconstruction. 

To obtain the differential cross-section from the back- 
ground subtracted distribution of observed Mtt values, 
acceptance effects and smearing effects from the recon- 
struction need to be corrected for. The required accep- 
tance correction is computed from signal simulation with 
Pythia. Factors to correct for differences between data 
and Monte Carlo observed in control samples are applied 
for the lepton identification and &-jet identification rates. 
The distortions of the reconstructed distribution are un- 
folded using using the singular value decomposition 224] 
of the response matrix that is obtained from simulations. 

Relevant systematic uncertainties arise from the back- 
ground normalisation, the acceptance, parton density dis- 
tributions, the used Monte Carlo generator and jet en- 
ergy scale. The relative contributions of the uncertainties 
strongly depend on M^f. To reduce the uncertainty on the 
jet energy scale an in-situ calibration of the jet energy 
scale is performed. This uses the invariant dijet mass re- 
constructed from the two non-6-tagged jets most consis- 
tent with Mw- 

The differential cross-section obtained in 2.7 fb^^ of 
data using the semileptonic decay mode is shown in 
Fig. [321 Ell]- The consistency with the Standard Model 
expectation is computed using Anderson-Darling statis- 
tics |225j . The observed p- value is 0.28, showing good 
agreement with the Standard Model. 

Finding no evidence for physics beyond the Stan- 
dard Model limits on gravitons in a Randall-Sundrum 
model |226j decaying to top quarks are set using the 
CL^ method [HTUHSj. Signal is modelled with Mad- 
Graph-|-Pythia assuming a first resonance with a mass 
of 600 GeV. The Anderson-Darling statistics is used as test 
statistics in the CL^ method. For the ratio of the warping 
parameter over the Planck mass CDF finds k/Mp\ < 0.16 
at 95% C.L, see Fig. |3H1 
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Fig. 38. Expected and observed limit on k/M-pi in a Randall- 
Sundrum model obtained from dau/AMu |221ll222j . 



The invariant top quark pair mass was used in further 
analyses by both CDF and D0 to search for new physics. 
These results are described in Section ITU 



3.8 Gluon Production vs. Quark Production 

Top pair production at the Tevatron with ^/s — 1.96 TeV 
takes place either through quark anti-quark annihilation 
or through gluon fusion. The former is expected to domi- 
nate with the gluon fusion contributing about 15%. Due to 
the large uncertainties of the large-x gluon density in the 
proton the exact size of the gluon contribution is rather 
uncertain [gOllHTlIM] . 

Two properties of the two production processes allow 
to separate them and to measure their relative contribu- 
tions. Close to threshold the spin states of the gluon fu- 
sion are J = 0, J^ = 0, while the qq annihilation yields 
J = 1, Jz = ±1 [^229]. This yields angular correlations 
between the charged leptons in the dilepton channel. Al- 
ternatively one can exploit the difference in the amount 
of gluon radiation from quarks and gluons: The gluon fu- 
sion processes are expected to contain more particles from 
initial state radiation. CDF has used both features to mea- 
sure the gluon fraction of top quark pair production. 



3.8.1 Angular Correlation Methods 

Dilepton Channel CDF has investigated 2.0 fb~^ of data 
with an event signature of top quark pair dilepton 
events [230j . The selection requires two oppositely charged 
leptons, at least one of which must be isolated, and at least 
two jets. The scalar sum of the lepton and jets transverse 
energies must exceed 200 GeV. Additional cuts are placed 
to reject cosmic particles, leptons from photon conversion 
and Z boson events. 

The azimuthal angle between the two leptons is mea- 
sured in each event. Then a template method is used to 
measure the fractional contribution of the different pro- 
duction mechanisms. The expected behaviour of signal 
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Fig. 39. Distribution of azimuthal angle between the two lep- 
tons, A(j), observed in 2.0 fb~^ of CDF data, compared to the 
best fit template curves 



events is simulated using Herwig with the top quark mass 
set to rrit = 175 GeV and the CTEQ5L parton distribution 
function. Pythia and Mc@nlo are used in systematic 
studies. Backgrounds are dominated by diboson produc- 
tion and Z boson events with tauonic decay. These are 
simulated by Pythia. In addition the background from 
events with only one true lepton and a jet misidentified 
as lepton are described using data. All simulated events 
are passed through the full CDF detector simulation and 
reconstruction. 

The angular distributions obtained for events pro- 
duced by qq annihilation and gg fusion and the sum of 
backgrounds are separately fitted with smooth functions, 
that then serve as signal and background templates. 

The measured fraction of top quark pairs produced 
through gluon fusion is obtained from an unbinned 
likelihood fit of these templates to the observed data, 
c.f. Fig. [39] (left). Systematic uncertainties include uncer- 
tainties on the template shapes, the acceptance differences 
between qq annihilation and gluon fusion, the used matrix 
element, initial and final state radiation and PDF uncer- 
tainties. The uncertainties are determined as function of 
the nominal gluon production fraction and several of them 
may contribute up to 10%. All systematic uncertainties 
are included in the determination of the Feldman- Cousins 
band, see Fig. |39] (right), which is used to obtain the final 
result with errors. 

In the investigated 2.0 fb^ of dilcpton events CDF 



obtains a gluon fusion fraction of 0.53 ± 0.37 [230) . The 
total uncertainty is dominated by statistical uncertainties 
and is not yet able to restrict the theoretical uncertainties 
on the gluon fusion production. 



Lepton plus Jets Channel Angular correlations are also 
used by CDF in an analysis of 0.96 fb~^ with lepton plus 
jets events (23 1| . The events are required to contain one 
energetic lepton, large missing transverse energy and at 
least four jets. One of the jets is required to be identified 
as b jet through the presence of a secondary vertex. 



In each event the decay chain of the top quark pair 
decay is reconstructed from the four leading jets using a 
kinematic fit with constraints on the W boson and the top 
quark mass, c.f. Section [2. 2. II Only jet parton assignments 
which associate the tagged b jet(s) with a b quark are 
considered. The one with the best fit quality is used for 
further analysis. 

From the reconstructed top quark pair decay, eight ob- 
servables are used to feed a neural network that is trained 
to distinguish between qq and gg production. The observ- 
ables are the (cosine of the) angle between the top quark 
momentum and the beam direction of the incoming pro- 
ton, the top quark velocity and the (cosines of the) six 
angles of the top quark decay product as defined in the 
off-diagonal spin basis, c.f. Section 13.51 The first two ob- 
servables are contributing about one third of the total sen- 
sitivity, each. The six decay angles yield the remaining 
third. 

The neural network is trained separately for events 
with one and with more than one identified b jet. Simu- 
lation of the two signal production processes is done with 
Herwig, the dominant background of VF-hjets is gener- 
ated with Alpgen-I-Herwig. The generated events are 
passed through the full detector simulation and recon- 
struction chain of CDF. To obtain the measured gg pro- 
duction fraction, fgg, templates for the neural network 
output are constructed from the simulation as function 
of fgg and a likelihood to observe the measured data as 
function of fgg is maximised. The Feldman- Cousins ap- 
proach |161| IS] is applied to restrict the final result to the 
physically allowed range. 

To determine the systematic uncertainties various 
pseudo experiments with systematically varied signal 
and/or background templates are studied. The deviation 
of fgg obtained in these samples from the standard tem- 
plate result is considered as systematic uncertainty. The 
dominant uncertainty is found to stem from background 
shape and composition as well as from the differences be- 
tween leading and next-to-leading order simulation of the 
signal. 

In 0.96 fb~^ of lepton plus jets events CDF is able to 
set limits on the gg production fraction of fgg < 0.61 at 
95%C.L. [231) . This method yields independent informa- 
tion, but is not as sensitive as the following method. 



3.8.2 Soft Track Method 

Another method that CDF applies to measure the fraction 
of top quark pair production through gluon fusion relies on 
differences that occur because gluons have a higher prob- 



ability to radiate than quarks |232j . The analysis is based 
on top quark pairs with semileptonic decays in 0.96 fb~ 
of CDF data. 

As sensitive observable the average number of soft 
tracks per event, (iVtrk), with 0.9 GeV < pr < 2.9 GeV 
in the central detector region \r]\ < 1.1 is used. In simu- 
lation this is shown to have linear relation to the average 
number of gluons, (Ng), in the hard process. This relation 
is calibrated on two samples with different gluon content: 
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W + Ojets events for low gluon content and dijet for high 
gluon content. W + Ijet events are used as cross-check. 

The dijet event sample for calibration is selected re- 
quiring a leading jet with transverse momentum between 
80 and 100 GeV and a second recoiling jet with |Z\(/)| > 
2.53. Vetoes are applied on lepton candidates and missing 
transverse energy. The W^-|-jets sample is selected requir- 
ing an isolated lepton, large missing transverse energy. For 
the signal top quark pair sample in addition at least four 
jets are required. At least one of these jets must be iden- 
tified as 6-jet. In the W^-f jets and top quark pair samples 
vetoes on additional lepton candidates and on leptons con- 
sistent with photon conversion or cosmic rays are applied. 

After calibration of the relation between (A'trk) and 
{Ng), the fraction of events with a high gluon content, /g, 
is determined using a binned likelihood fit. The fit result 
is corrected according to the expected background contri- 
bution. The fraction of events with high gluon content in 
the background, fg^^, is extrapolated from the VF-|-jets 
sample with up to three jets to the four or more jet sam- 
ple. The expected amount of background in the selected 
signal sample is determined following the neural network 
based method in (223) . The obtained high gluon content 
in top quark pair events, /**, in a last step is corrected 
for the differences in acceptance between the gluon fusion 
and the qq annihilation processes. 

The systematic uncertainties of this measurement are 
dominated by uncertainties of the calibration procedure 
and were determined by varying the corresponding pa- 
rameters in the analysis. 

In the dataset of 0.96 fb~ CDF determines a gluon fu- 
sion fraction of top quark pair production in semileptonic 
events of 0.07± 0.14stat ± O.OTsyst UM- This number cor- 
responds to an upper limit of 0.33 at 95% C.L, well in 
agreement with the Standard Model expectations. Also 
this measurement is statistically limited. A combination 
of this result with the result of the Angular Correlation 
Method described in the previous section yields an about 
10% improvement on the upper limit '23 1| . 



3.9 Top Quark Width and Lifetime 

The top quark width and its lifetime are related by Heisen- 
berg uncertainty principle. In the Standard Model the top 
quark width is expected to be 1.34 GcV corresponding to 
a very short lifetime of about 5 • 10^^^ s. Experimentally, 
these predictions have been challenged for deviations in 
very different analyses. CDF constrains the lifetime from 
the distribution of reconstructed top quark mass values 
and from the distribution of lepton track impact parame- 
ters. D0 combines the measured t — ^ Wb branching frac- 
tion and the single top ^-channel cross-section to measure 
the top quark width. 



values from top quark pairs decaying to lepton plus jets 



in up to 4.3 fb^ of data [23311234) . After selecting events 
with one lepton, large missing transverse momentum and 
at least four jets. One of the jets is required to be identified 
as 6-jet. 

In these events the top quark mass is reconstructed 
using a kinematic fit that determines the four-momenta 
of the top quark decay productions from the measured jet 
and lepton momenta and the transverse missing energy. 
The fit uses constraints that force the W boson decay 
products to build the W boson mass within the width of 
the W boson and the reconstructed top and anti-top quark 
masses to be equal within the top quark width. In the am- 
biguous association of jets to partons identified &-jets are 
only associated to 6-quarks. Of the remaining associations 
and the two solutions for the neutrino 2-momentum, the 
one with the best x^ is used. It was checked that the use 
of the constraint of the equality of the top quark masses 
width does not destroy the sensitivity to the true width. 

To find the measured value of the top quark width the 
distribution of top quark masses reconstructed with the 
best association in each event is compared to parametrised 
templates with varying nominal width. Templates for top 
quark pair signal events were generated using Pythia 
with rrit — 172.5 GeV |234) . Background contributions of 
M^-l-jets are modelled Alpgen-|-Pythia. Multijet contri- 
butions from data with non-isolated leptons. Single top 
quark and diboson events are simulated with MadGraph. 
The template distributions for discrete values of the nom- 
inal top quark width are parametrised to obtain smooth 
template functions that can now be interpreted as prob- 
ability densities. The measured top quark width is deter- 
mined in an unbinned likelihood fit. 

Recently, in addition to the reconstructed top quark 
mass, the invariant mass of the jets assigned to the 
hadronic W decay is considered as an observable. This 
allows a simultaneous fit of the top quark width and the 
jet energy scale [234) . 

The Feldman-Cousins approach |161l l6) is used to de- 
termine the final result excluding the unphysical values 
of negative widths that may occur in the fit. The jet res- 
olution followed by colour reconnection effects yield the 
biggest single contribution to the systematic uncertain- 
ties, are propagated to the final Feldman-Cousins band by 
convoluting their effects with the fitted width function. 

Including all systematics this CDF analysis of 4.3fb^ 
yields an upper limit of the top quark width Ft < 7.5 GeV 
at 95% C.L. which corresponds to n > 8.7 • lO^^^^s. At 
68%C.L. the top quark width is determined as 0.4 GeV < 
n < 4.4 GeV [2M1. 



3.9.2 Lepton Impact Parameter 



3.9.1 Top Quark Mass Distribution 

The limit on the top quark width was obtained by CDF 
from the the distribution of reconstructed top quark mass 



The limit from the lepton track impact parameter distri- 
bution was obtained by CDF using lepton plus jets events 



in 318 pb^ of data [235) . Events are selected requiring 
one isolated lepton, missing transverse energy and at least 
three jets. At least one of the jets has to be identified as 
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6-jet. The lepton track needs to be reconstructed with at 
least three R-(j> positions in the CDF sihcon tracker. 

The lepton impact parameter, do, chosen as observable 
in this measurement is defined as the smallest distance 
between the collision point and the lepton track in the 
transverse projection. The collision point is computed as 
the position of the beam line in the transverse plane at 
the reconstructed z position of the primary vertex. 

The distribution of lepton impact parameters expected 
in an ideal detector for various top quark lifetimes is sim- 
ulated with Pythia. The resolution of the CDF detector 
is measured in Drell-Yan data near the Z boson resonance 
and used to derive the templates for the real detector ex- 
pectation. The dominant backgrounds like M^-fjets con- 
sist of prompt leptons, which are described by the zero 
lifetime template. But the distribution of multijet events, 
backgrounds with r leptons and electrons from photon 
conversions need to be modelled. Multijets and electron 
conversions are modelled from control samples in data. 
Backgrounds with r leptons are modelled using Herwig. 

From these templates a likelihood as function of CTt 
is built. The maximal likelihood is obtained from the 
fim template. Systematic uncertainties on the signal and 
background systematics are computed with correspond- 
ingly varied templates and are dominated from the un- 
certainty on the detector resolution for prompt leptons. 
The Feldman- Cousins approach is used to determine the 
observed limit of cTt < 52.5 /im at 95% C.L. 



3.9.3 Branching fraction and single top cross-section 

The total width of the top quark can be written as the 
ratio of the partial width for the decay t — ^ Wb and the 
corresponding branching fraction: 



Ft = r{t -> Wb)/B{t ^ Wb) 



(87) 



D0 measures the top quark width by relating the par- 
tial width and branching fraction to the results of two 



independent analyses |236j . The branching fraction in the 
denominator is equal to the branching fraction ratio, i?;,, 
measured from top quark pairs (c.f. Section [3. 2. 2p . This 
assumes that the top quark always decays to a W^-boson 
plus a quark. The partial width is derived from the sin- 
gle top production cross-section which is proportional to 
r{t —> Wb). When considering only the t-channel cross- 
section, (Tt-channcb this proportionality is valid also in the 
presence of anomalous couplings in the iVF6-vertex. Thus 
the partial width is determined as 



r(t -^ Wb) = <Tt_channcl 



n^jt^Wb) 

crSM 
i — channel 



where the superscript SM indicates the Standard Model 
expectations. As the total width is by definition larger 
than the partial width, a lower bound of ri > 1.21 GeV at 
95% C.L. can be set from the i-channel production cross- 
section jl90) alone. In the combination of the partial width 
with the branching fraction systematics are classified and 



treated as either fully correlated or uncorrelated. They are 
dominated by uncertainties in the background description 
and the description of the 6-jet identification in the input 
analyses. 

With the combination of the t-channel production 
cross-section measured in 2.3 fb~^ [190] and the branching 
fraction ratio measured in 1 fb~ |162j D0 obtains a top 
quark width of Ft = 1-99^^55 GeV corresponding to a top 
quark lifetime of n = (3.3J:oJ) • IQ-^s s ^55] . 

3.10 Outlook to LHC 

Within the SM the interaction properties of the top quark 
are fully defined once the top quark mass is known. The 
verification of the predicted properties establishes the top 
quark as the particle expected in the SM. Huge progress 
has been made in the last years in experimentally verifying 
the expectations. Measurements of the W boson helicity, 
the CKM element Vtb and searches for FCNC confirm the 
expected weak interaction properties. The electric charge 
has been challenged and the strong interaction properties 
have been verified in differential cross-section measure- 
ments, in a verification of the contribution of gluon fusion 
processes and also in the forward backward charge asym- 
metry. First tests of the spin structure and determinations 
of the top quark width complete the current picture. All 
current results are compatible with the SM expectations, a 
single deviation of 3.4(T appears in the charge asymmetry 
at high Mti. 

Despite the great progress, the precision of the ex- 
perimental knowledge of interaction properties of the top 
quark is still limited. Many results have not yet reached 
a precision of 10%. As a consequence of the experimental 
progress even the more precise results, like the W boson 
helicity and the measurement of Vtb, sue statistically lim- 
ited. 

The near future will see updates of the results with 
the yet unanalysed Tevatron data to a total 10 fb~ . In 
addition the LHC schedule envisions to collect Ifb^^ of 
proton-proton collisions at ^/s = 7 TeV. The LHC experi- 
ments will then have about twice as many top quark pairs 
as the Tevatron and even three times as many single top 
quarks. For the top quark pairs the production through 
gluon fusion and for the single top quark the t-channel 
production will dominate. When systematic uncertainties 
can be controlled similarly well as at the Tevatron, the 
LHC results on this dataset will still be limited by statis- 
tics. 

Only the update of the LHC to the design centre-of- 
mass energy of 14 TeV will provide sufficiently many top 
quarks to enter an area of precision measurements in the 
verification of top quark interactions. 

An important exception to these statements is the 
measurement of the forward-backward asymmetry, which 
cannot be measured in pp collisions. Instead related asym- 
metries need to be studied [25711^1111^17 . 238 , 239 . These 
measurements will probably require large luminosity as 
only quark annihilation diagrams contribute to the asym- 
metries. 
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4 New Particles in Top Quark Events 

The phenomenology of the top quark may also be altered 
by particles that are not expected within the Standard 
Model, but in one of the many models of new physics. Such 
particles beyond the Standard Model may occur in the top 
quark production or its decay, depending on the specific 
model or its parameters. Some models of new physics also 
contain new particles with signatures that are very similar 
to the Standard Model top quark. The Tevatron experi- 
ments have checked for all these different extensions of the 
Standard Model in the top quark sector. 

This subsection will actually start with a process that 
is expected in the Standard Model though at very low 
rate: associated Higgs production. In some models of new 
physics this process is expected to be enhanced. Then par- 
ticles beyond the Standard Model in the top quark decay 
will be discussed. Finally, searches for the production of 
particles that look like the top quark but are not are de- 
scribed. 



4.1 Associated Higgs Boson Production, ttH 

Top quark pair production may be associated by the pro- 
duction of a Higgs boson. For parameters where Higgs 
bosons dominantly decay to bottom quark pairs, i.e. low 
Higgs masses, this associated production is a possibility to 
measure the top quark Yukawa coupling. While the corre- 
sponding cross-section in the Standard Model is too low to 
allow a Higgs discovery in this channel alone, it still con- 
tributes to the combination of the Standard Model Higgs 
searches. In some models including new physics an en- 
hancement of ttH production is expected |240ll241|l242j . 



D0 



D0 performed an analysis searching for associated Higgs 
production in events with a lepton (e or /i) missing trans- 
verse energy and at least four jets [243] . The analysis uses 
the scalar sum of transverse momenta, Ht, the number 
of jets and the number of jets identified as 5-jets to dis- 
criminate the Standard Model backgrounds and top pair 
production containing no Higgs from the signal. 

Signal events are simulated using Pythia. For tt pro- 
duction pure Pythia simulation was compared to Alp- 
GEN-I-Pythia simulation. Due to the difference between 
the two simulations a 50% uncertainty was assigned to the 
contribution of ttbb through QCD processes. Background 
from W-|-jets events is simulated with Alpgen-|-Pythia 
and normalised to data. Multijet background was com- 
pletely estimated from data. Smaller backgrounds are 
taken from simulation normalised to NLO cross-sections. 

In the investigation of 2.1 fb~ the observed data agree 
with the Standard Model expectations within statistical 
and systematic uncertainties [243) . To compute limits sig- 
nal and background contributions are fitted to the data for 
a background only assumption and for a signal plus back- 
ground assumption. Limits on a{ttH)-B{H — >■ bb) are then 



derived using the CL^. method |227U228j for Higgs masses 
between 105 and 155 GeV. 

For Mh = 115 GeV the cross-section limit corresponds 
to about 60 times the Standard Model value. While this 
allows to exclude unexpectedly large Higgs boson produc- 
tion in association with the top quark, its contribution to 
the Standard Model Higgs search remains small. 



4.2 Charged Higgs Boson 

Particles beyond the Standard Model in the final state of 
top quark pair events may alter the branching fractions 
of the various top quark decay channels and modify the 
kinematic properties of the final state. 

Charged Higgs bosons appear in many extentions of 
the Standard Model due to the need for an additional 
Higgs doublet with a separate vacuum expectation value. 
These models are characterised by the ratio of the vac- 
uum expectation values of the two Higgs doublets, tan/3. 
A charged Higgs boson can replace the W boson in top 
quark decays. Because charged Higgs bosons have different 
branching fractions than W bosons this alters the branch- 
ing fractions to the various top quark pair decay channels. 
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Fig. 40. Charged Higgs boson branching fraction 
as function of tan/? [244] . 
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Fig. 41. Results from recasting the CDF top quark cross-section measurements [244j . Left: Exclusion region in the MSSM 
mn-tan /3-plane for an example benchmark scenario corresponding to the parameters indicated below the plot. Right: Upper 
limits on B{t — >■ H^h) derived without assumptions on the charged Higgs boson branching fraction. 



If its mass is different from the Mw it also modifies the 
kinematic properties of the top quark pair final state. 

In the Minimal Supersymmetric Standard Model 
(MSSM) |245) the decay at low tan/3 is dominated by 
hadronic decay to cs at low Higgs boson masses and to t*b 
for Higgs boson masses above about 130 GeV. For tan/3 
larger than about 1 a leptonic decay to tv dominates, 
c.f. Fig. 001 The figure also shows the expected branching 
fraction of t — > H^h which is especially large for very low 
and very high tan/3 and rather small in the intermediate 
range. 



CDF 

CDF has performed two analyses with different ap- 
proaches. An analysis based on 0.2 fb~ uses the CDF 
tt cross-section measurements in various channels and re- 
casts the interpretation to obtain limits on the charged 
Higgs boson production. A second more recent analysis 
on 2.2 fb~ investigates the kinematic differences between 
lepton plus jets events from top quark pair production 
with Standard Model decay and those including charged 
Higgs boson decays. 



Recast of Top Quark Pair Cross-Section To obtain lim- 
its on a possible charged Higgs boson contribution in 
top quark decay CDF utilises cross-section measurements 
performed in the lepton plus jets channel (with exactly 
one 6-tag or two or more h tags), the dilepton channel 
and the r plus lepton channel [244] . Care is taken to 
avoid overlap between the various channels. Beside the 
Standard Model decays of a top quark through a W bo- 
son, four decay modes through the charged Higgs boson 
are considered: H'^ — >■_ fv, 7?+ — >■ cs, H'^ -^ t*b and 
i7+ -^ W^h^ — > W^bb. The latter has a non- negligible 
contribution at intermediate values of tan /3. 



Selection efficiencies are taken from simulation of top 
quark pair events for various masses of the top quark, 
the charged and neutral Higgs boson, h^ . The simula- 
tion takes the dependence of the width of the top quark 
and the charged Higgs boson into account. The produc- 
tion cross-section is kept at its Standard Model value for 
rat = 175 GeV: crtt = 6.7 ± 0.9pb. 

The event counts observed in data in the four chan- 
nels are compared to the expectations in three different 
ways. For specific benchmarks of the MSSM a Bayesian 
approach is used to set limits on tan/3. This analysis uses 
a fiat prior on log tan /3 within the theoretically allowed 
range. These limits are computed for various values of 
the charged Higgs boson mass and five different parame- 
ter benchmarks. Figure Hi] (left) shows the results for one 
specific benchmark. 

For the high tan/3 region 77+ — )■ fv dominates in a 
large fraction of the MSSM parameter space. Setting the 
branching fraction of i?+ — > fv to 100%, limits on the 
charged Higgs contribution to top quark decays are set 
using Bayesian statistics. A fiat prior for B{t — >■ H'^b) 
between and 1 is used. For charged Higgs boson masses 
between 80 GeV and 160 GeV CDF can exclude B{t -^ 
H+b) > 0.4 at 95%C.L. 

Finally, a more model independent limit is computed 
by scanning the full range of possible charged Higgs bo- 
son decays. For all five H^ decay modes considered the 
branching fraction is scanned in 21 steps, assuring that 
the sum of branching fractions adds to one. Limits on 
B{t ^ H^b) are computed for each combination. The 
least restrictive limit is quoted. Also this analysis is re- 
peated for various charged Higgs boson masses. The lim- 
its obtained in this more general approach, shown in 
Fig. UT] (right), exclude only very high contributions of 
charged Higgs bosons to top quark decays of above ap- 
proximately 0.8 to 0.9, depending on the charged Higgs 
boson mass. 
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Fig. 42. CDF dijet mass distribution with 120 GeV Higgs bo- 
son events assuming B{t — > H^b) = 0.1. The size of Higgs 
boson signal corresponds to the expected upper limit branch- 
ing ratio at 95% C.L. for 120 GeV [246] . 



Investigation of Kinematic Differences At low tan/3, 
where the charged Higgs boson can also decay to cs, CDF 
used the invariant dijet mass to search for a possible H^ 



contribution in top quark pair events |246) . Lepton plus 
jet events are selected requiring at least two of the four 
leading jets in px to be 6-tagged. The four leading jets 
are used in a kinematic fit that requires consistency of the 
fitted lepton and neutrino momenta with the W boson 
mass and reconstructed top quark masses to be 175 GeV. 
The dijet mass of the hadronic W decay remains uncon- 
strained. The jet parton assignment with the best x^ is 
used and the charged Higgs boson mass is reconstructed 
from the non-6-taggcd of the 4 leading jets. For events 
with more than 4 jets the 5th jet is added to its closest 
neighbour if their AR < 1.0 to improve the dijet mass 
resolution. 

Background events are dominated by top quark pair 
production with Standard Model decay. Further processes 
included are Vl^-|-jets, Z'-l-jets, diboson, single top quark 
and multijet events. Except for multijet events the back- 
grounds are estimated from simulation. The normalisation 
for iy-|-jets taken from data, for the others it is taken from 
theory. The multijet background is fully determined from 
data. 

To determine a possible contribution of charged Higgs 
boson in the decay of top quark pair production a binned 
likelihood fit is performed. The likelihood is constructed 
with templates for the backgrounds and using the branch- 
ing fraction of top quark to charged Higgs bosons, the 
number of top quark pair events and the number of back- 
ground events as parameters. The number of background 
events is constrained within the uncertainty to the expec- 
tation. The observed dijet mass distribution and the fitted 
background composition is shown in Fig. |42] including a 
charged Higgs boson contribution of 10%. 

Systematic uncertainties are computed by fitting 
pseudo data created from systematically varied templates 
with the standard unshifted templates. The change of the 
branching ratio due to the systematic variation is taken 
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Fig. 43. CDF limits B{t -^ H+b) observed in 2.2 fb"^ data 
(red dots) compared with the expected limit assuming the 
Standard Model (black line with uncertainty bands) 



as systematic uncertainty for each variation considered. 
These uncertainties are included to a final likelihood by 
convoluting a Gaussian distribution to the original likeli- 
hood. Systematic uncertainties are dominated by the sig- 
nal modeling that is derived from replacing the default 
Pythia sample of ti events by a Herwig sample. For Higgs 
boson masses close to the W boson mass the jet energy 
scale uncertainty becomes dominant. 

For various assumed Higgs boson masses 95%CL limits 
on the branching fraction are determined by integrating 
the likelihood distribution to 95% of its total. As shown 
in Fig. US limits between about 10% and 30% can be set 
depending on the mass of the charged Higgs, consistent 
with the expected limits for pure Standard Model top 
quark decays. This result is less model dependent than 
the above CDF limits, but not as strict within the models 
used above. 



D0 



The D0 collaboration searched for a light charged Higgs 
boson contribution in top quark decay by reinterpreting 
the cross-section measurements in various decay channels 
and for heavy charged Higgs boson contributing to single 
top quark production. 



Charged Higgs Boson in Top Quark Decay In the search 
for a light charged Higgs boson in top quark decay D0 uses 
the cross-section analyses for lepton plus jets, dilepton and 
lepton plus tau decay channels, with lepton refering to e 
and /i only [126.248,247] . The channels are kept disjoint 
and further separated into subsamples depending on the 
number of jets, the number of h tags and the lepton type. 
The number of expected events for each of the subsamples 
is computed from ti simulation with Pythia including the 
Standard Model decays and decays of the top quark to a 
leptophobic or a tauonic charged Higgs boson and charged 
Higgs boson masses between 80 and 155 GeV. 
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Fig. 44. Limits on the contribution of a charged Higgs boson in top quark decays for a leptophobic model (left) and a tauonic 
model (right) obtained in 1 fb^^ of D0 data assuming the SM ti production cross-section [247| . 



A likelihood for the observed data is built each of the 
two models given the number of expected events as func- 
tion of the branching fraction B{t -^ H^b). The observed 
B{t — )■ H^b) is extracted by maximising the likelihood. 

In a first iteration the production cross-section is fixed 
at a value of au — 7.48 pb. Limits are set according to the 
Feldman-Cousins procedure including systematic uncer- 
tainties. The systematic uncertainties in this method are 
dominated by uncertainties due to the assumed top quark 
pair cross-section, the luminosity and 6-jet identification. 
The resulting limits obtained with lfb~^ of data exclude 
a branching fraction above around 20% for the pure lepto- 
phobic model and above 15 — 20% for the tauonic model, 
c.f. Fig.iilMZ]- 

In a second iteration the tt production cross-section 
is treated as a free parameter and determined simul- 
taneously with the limits on charged Higgs production, 
see Fig. HD This reduces the assumptions made in the de- 
termination of the limit. In addition the result is much less 
sensitive to the luminosity. In this method the description 
of multijet background becomes the largest systematic un- 
certainty. Such a two dimensional fit is only possible for 
the tauonic decay model where this analysis includes chan- 



nels that get enhanced in the presence of a charged Higgs 
and others that get depleted. In this channel the sensi- 
tivity enhances by more than 20% for low Higgs boson 
masses I247L 



D0 further discusses the implication for various su- 
persymmetric models and sets exclusion limits on these 
models in the plane of model parameters tan f3 and Mh± ■ 



Charged Higgs Boson in Single Top Quark Production 

Heavy charged Higgs bosons can not only occur in the 
decay of top quarks but may contribute to single top quark 
production. Their signature is identical to Standard Model 
s-channel single top quark production, but may have a 
resonant structure in the invariant mass distribution of 
its decay products, the top and the bottom quarks. 

Following their single top quark analysis, D0 selects 
events with an isolated lepton, missing transverse energy 
and exactly two jets, one of which is required to be identi- 
fied as &-jet |249| . Background estimation for W-|-jets and 
tt production is simulated using Alpgen. Standard Model 
single top quark production is modelled using Single- 
Top |191II192) . Charged Higgs boson signal events are sim- 
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Fig. 45. Limits on the contribution of a charged Higgs boson 
in top quark decays for the tauonic model obtained in 1 fb~^ of 
D0 data in a simultaneous fit of the SM tt production cross- 
section 
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Fig. 46. Observed and expected invariant mass of the W boson 
and the two leading jets in 0.9 fb"^ of D0 data [249| . 
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Fig. 47. Exclusion areas derived from 0.9 fb ^ of D0 data for 
Type I two Higgs Doublet Model (2HDM) [249]. 

ulated with a narrow width for the charged Higgs boson 
using CompHEP. Three types of two Higgs doublet mod- 
els (2HDM) are considered. In the Type I 2HDM one dou- 
blet gives mass to all fermions; in the Type H model one 
doublet gives mass to the u-type quarks (and neutrinos), 
the other to the d-type quarks and charged leptons. This 
model is realised in the MSSM. In the Type HI 2HDM 
both doublets contribute to the masses of all fermions. 
Due to the different couplings the cross-section of single 
top quark production in these 3 models is quite different. 
Standard Model and charged Higgs boson production 
of single top quarks is separated by reconstructing the 
invariant mass of the two jets and the W boson. This 
distribution shows good agreement between data and the 
Standard Model expectation, see Fig. [311 Bayesian statis- 
tics is used to set limits on the allowed cross-section for 
single top quark production through a charged Higgs bo- 
son. For the Type I 2HDM some region in tan /3 vs. M^± 
can be excluded, c.f. Fig.|47l a significant fraction of phase 
space is not accessible by the analysis in its current form 
due to the restriction to small H^ decay widths 



4.3 Heavy Charged Vector Boson, W' 

New charged gauge bosons, W, are expected in extensions 
of the Standard Model with additional gauge symmetries 
and in supersymmetric models, see e.g. [6^j 245] . Its cou- 
plings may be to left-handed fermions, like for the Stan- 
dard Model W boson, or include right-handed fermions. In 
general a mixture of these two options is possible. If the 
W' boson has left handed couplings, it will have a size- 
able interference with the SM W^ boson [250j . For purely 
right-handed couplings, a leptonic decay may only occur 
when the right-handed neutrinos are lighter than the W 
boson. In this case the decay to a top and bottom quark is 
an interesting channel to perform direct searches for such 
W' bosons. 

Both CDF and D0 search for various types of W' 
bosons decaying to tb pairs in conjunction with their single 
top quark analyses. The main discriminating observable 
is the reconstructed invariant mass of the decay products, 
which was also utilised to search for a heavy charged Higgs 
boson, c.f. Section [¥^ 



D0 



D0 has published a search for a heavy W' boson with 



decay to top and bottom quarks using 0.9 fb~ |251j . The 
event selection follows the single top quark analysis and 
requires one isolated lepton, missing transverse momen- 
tum and two or three jets, one of which must be identified 
as 6-jet. 

The invariant mass, vs, of the bottom and the top 
quark decay products is computed from the measured 
four-momenta of the leading two jets, the charged lep- 
ton and the neutrino. The transverse momentum of the 
neutrino is identified with the transverse missing momen- 
tum, its z-component infered by solving M^ = (pi + p„) 
choosing the solution with the smaller \pl\. 

The distributions expected within the Standard Model 
from a combination of simulation and data. Single top 
quark and top quark pair production is generated with 
SingleTop and Alpgen-|-Pythia normalised to their 
theoretical cross-sections. W^-|-jets background is gener- 
ated with Alpgen-(-Pythia and normalised to data be- 
fore 6-tagging in a way that it includes diboson back- 
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Fig. 48. D0 results on a search for W' boson decaying to top 
and bottom quark using 0.9 fb^^ of data. Top: Expected and 
observed limits on a left-handed W' production cross-section 
times branching fraction of the decay to top and bottom quark 
as function of Mw' compared to the theory prediction. Bottom: 
Same but for right-handed W' production [251) . 
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grounds. Also the W^+heavy flavour fraction is derived 
from data. Multijet background is fully taken from data. 
Samples of W boson events with masses up to 900 GcV are 
generated in conjunction with the single top quark sam- 
ples taking interferences with the W boson into account 
that are present for the left-handed W bosons. Because 
of this interference the Standard Model single top quark 
production in the s-channel is treated as part of the signal 
in the search for ]¥[ bosons. 

The distribution of reconstructed v^ measured by D0 
agrees with the expectation from the Standard Model. 
Limits on a possible contribution from W' bosons decaying 
to top and bottom quarks are derived as a function of Mw' 
assuming couplings like in the Standard Model, though 
possibly to right handed fermions. D0 uses the Bayesian 
approach with a flat non-negative prior on the cross- 
section times branching fraction. Expected and observed 
results are shown in Fig. 1481 Comparing upper limits on 
the W boson cross-section times branching fraction to top 
and bottom quark to the NLO theory predictions 252 
excludes left-handed W bosons with Mw < 731 GeV. If 
only hadronic decays are allowed the right handed W^ bo- 
son is excluded for M^y' < 768 GeV, when leptonic decays 
are also possible the limit is 739 GeV. 

Without assuming the coupling strength the Bayesian 
approach is used to determine a limit on the size of this 
coupling relative to the Standard Model, see Fig.|49l These 
limits assume no interference between the Standard Model 
W and the W bosons. 

In computing the above limits systematic uncertainties 
are included. They include effects due to uncertainties on 
the integrated luminosity, the theoretical cross-sections, 
branchings fraction, object identification efficiencies, trig- 
ger efficiencies, fragmentation models, jet energy scale and 
heavy flavour simulation. 



CDF 



In an investigation of 1.9 fb^ of data |253j CDF selects 
W^+jets events requiring one Icpton (e, ^) isolated from 
jets, missing transverse energy and two or three energetic 
jets. At least one of the jets must be tagged as 6-jet. 

In these events the neutrino momentum, p,^, is infered 
from the missing transverse momentum and by solving 
^w — (Pi +Pi') foi' the longitudinal component of the 
neutrino. The W boson mass, Mw, is set to its nominal 
value, pi is the measured lepton momentum. In case of 
complex solutions CDF assigns the real part of the solu- 
tion to the longitudinal neutrino momentum. The invari- 
ant mass of the lepton, the neutrino and the two leading 
jets, Mwjj, is then used as a discriminating observable. 

The distribution expected within the Standard Model 
is computed from a combination of simulation and data. 
The contribution of events containing a real W boson is 
taken from simulation. M^-|-jets samples are normalised 
to data before 5-tagging using a scale factor to correct 
the heavy flavour contribution to fit the observation in 
VF-|-ljet data. The other samples are normalised to their 
theoretical expectation. The identification of heavy fiavour 
jets is estimated from simulation and corrected with a 
scale factor. Misidentification of light flavour jets is com- 
puted from mis-tag rate functions. The contribution from 
events without real W bosons are estimated from events 
with electrons that pass only a subset of the full electron 
identification and thus are enriched with jets misidentified 
as electrons. 

W boson signal events are simulated using Pythia for 
W boson masses between 300 and 950 GeV with fermion 
couplings identical to the Standard Model W boson. When 
the right-handed W boson is heavier than the right- 
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Fig. 49. D0 results on a search for W' boson decaying to 
top and bottom quark using 0.9 fb^^ of data. Limits on the 
W' boson couphng relative to the Standard Model Vy-boson 
coupling [251) . 
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Fig. 50. Limits on W^ cross-section times branching fraction 
to tb as function Mjy/ compared to theory obtained in the 
CDF search for Wr using 1.9 ffa-^ of data [253] . 
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Fig. 51. Limits on the W^ couphng strength relative the 
Standard Model coupling as function M^/ . obtained in the 
CDF search for Wr using 1.9 fb"^ of data [253] , 



handed neutrinos, the branching fraction to iv is corrected 
according to the additional decay modes. 

Limits are constructed according to the CL^ 
method _227.:228^ . Probabilities are computed from pseudo 
experiments which are generated including variations due 
to systematic uncertainties. The dominating systematic 
uncertainties are the jet energy scale and the scale factor 
used to account for differences between simulation and 
data in the 6-tagging algorithm of CDF. 

Limits on the WJ^ boson production cross-section are 
set as a function of M^y assuming the Standard Model 
coupling strength. These are converted to mass limits by 
comparison to the corresponding theoretical expectation 
and yield A/^y > 800 GeV for W^ bosons which decay 
leptonically arfd Af^/^ > 825 GeV for A/^,^ > Mwj^- For 
the more general case that the W^ coupling is a priori un- 
known the W^ coupling strength, g', relative to the Stan- 
dard Model coupling, gw, is constrained. Limits are com- 
puted from the above analysis as function of the assumed 
Mw'- The observed and expected limits derived by CDF 
for Mw and g'/gw are shown in Fig. [501 and ICT 



4.4 Resonant Top Quark Pair Production 

Due to the fast decay of the top quark, no resonant pro- 
duction of top quark pairs is expected within the Stan- 
dard Model. However, unknown heavy resonances decay- 
ing to top quark pairs may add a resonant part to the 
Standard Model production mechanism. Resonant pro- 
duction is possible for massive Z-like bosons in extended 



gauge theories [254j , Kaluza- Klein states of the gluon or Z 
boson PSSII^ . axigluons gST], Topcolor pSSlI^ . and 
other theories beyond the Standard Model. Independent 



of the exact model, such resonant production could be 
visible in the reconstructed tt invariant mass. 



CDF 



CDF has employed several different techniques to search 
for resonances in the tt invariant mass distribution. All 
analyses use a very similar event selection: an isolated lep- 
ton, missing transverse energy and four or more jets. One 
analysis (the Matrix Element plus Template method) is 
also applied to the full hadronic channel selecting six or 
seven jets. Their main difference is the method to recon- 
struct the ti invariant mass distribution. 



Constrained Fit plus Template The method that uses 
the least assumptions reconstructs the invariant mass us- 
ing a constrained fit 260 and is performed requiring at 
least one identified b jet. In the fit the final state lepton 
and quark momenta are determined from the measured 
lepton momentum, the missing transverse energy (which 
is assumed to stem from the unseen neutrino) and the 
measured jet momenta. Constraints are imposed that re- 
quire the sum of neutrino and lepton momenta as well as 
the two light quark momenta to be consistent with the W 
boson mass. In addition these pairs in combination with 
one b quark need be consistent with the top quark mass of 
175 GeV. The fitted momenta may vary within the exper- 
imental resolution of their assigned measurement and the 
mass constraints are varied within the natural widths of 
the W boson and the top quark, respectively. The fit thus 
assumes the lepton plus jets decay topology of top quark 
pair events. 

Of the multiple jet parton assignments the one with the 
best x^ from the fit is used to compute the top quark pair 
invariant mass for each event. The expected distribution of 
this observable is dominated by Standard Model top quark 
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Fig. 52. Distribution of the invariant top quark pair mass re- 
constructed with a constrained fit in 1 fb~^ of CDF data 12601. 
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Fig. 53. Limits on the cross-section times branching fraction 
for resonant top quark pair production observed in lfb~^ of 
CDF data [260) . Theoretical curves are shown for various model 
and used to set mass limit on the corresponding resonance. 



pair events which are simulated using Herwig. Further 
backgrounds like VF+jets, misidentified multijet events, 
diboson and single top quark events are modelled with 
combination of simulation and control data. 

Templates for a resonant production of top quark pairs 
are simulated using Pythia for resonance masses between 
450 and 900 GeV. The resonance couplings are propor- 
tional to those of a Standard Model Z boson. The width of 
this Z' boson was kept at 0.012^/^'. The observed invari- 
ant top quark pair mass distribution is compared to the 
distribution expected from the Standard Model in Fig. W2\ 

From the expected distributions CDF constructs a like- 
lihood for the expected bin content of the distribution 
of the invariant top quark pair mass as function of the 
resonant production cross-section time branching ratio, 
(JxB, the number of Standard Model top quark pairs and 
the number of non-ii events. Nuisance parameters with 
Gaussian constraints are used to implement the effect of 
systematic uncertainties. These include uncertainties that 
affect the relative background normalisation and the lu- 
minosity, the uncertainty of the jet energy scale and the 
shape change due to the top quark mass uncertainty. Mi- 
nor contributions come from varying the PDFs between 
CTEQ6M [55] and MRST [IHI] parametrisations and the 
uncertainty on the strength of initial and final state radi- 
ation. 

To find the upper limits the maxima of the likelihood 
as function of axB is integrated to the point where the 
integral reaches 95% of its area. This is done for each as- 
sumed resonance mass between 450 and 900 GeV in 50 GeV 
steps. Expected and observed limits are shown in Fig. [531 
At high resonance masses the observed limits exclude a 
resonant top quark pair production with axB > 0.55 pb 
at 95% C.L. A production through leptophobic Topcolor 
assisted Technicolor is excluded for resonance masses up 
to 720 GeV. 



Matrix Element plus Template The resolution of the 
reconstructed invariant top quark pair mass can be im- 
proved by assuming additional information. In an analy- 
sis of 680 pb~ CDF employed the Matrix Element tech- 
nique to reconstruct the invariant mass distribution that is 
used to search for resonant production in lepton plus jets 
events [262) . Following the mass analysis ^3] described 
in Section 12.2.31 for each event a probability density, 
P({p}|{j}), is computed to find the momenta of the top 
quark pair decay products (4 quarks, charged lepton and 
neutrino, {p}) given the observed quantities, {ji}. This 
probability is computed from the parton density functions, 
the theoretical Matrix Element for Standard Model top 
quark pair production and decay and jet transfer func- 
tions that fold in the detector resolution. It is converted 
to a probability density for the top quark pair invariant 
mass using 

Pf{Mu\{j}) = Jd{p} Pi{p}\{j}) 5{Mu - m{{p})) . 

(89) 
The result for all possible jet-parton assignments is 
summed, before the mean value is computed and taken 
as the reconstructed invariant top quark pair mass for the 
event under consideration. Here the 6-tagging information 
is used to reduce the number of allowed jet parton assign- 
ments. The events are not required to contain 6-tagged 
jets in the event selection. 

Template distributions for Standard Model and reso- 
nant tt processes are derived from Pythia, W^-|-jets events 
from Alpgen-|-Herwig including full detector simula- 
tion. As a resonance signal a Z' boson with a width of 
1.2%AIz' was generated. The multijet background tem- 
plate was taken from data. The Standard Model tt sam- 
ples and diboson samples are normalised to the theoret- 
ical cross-section, while the sum of multijet and W^-|-jet 
samples are scaled to fit the observed data and depend 
on the assumed signal contribution. The resulting ex- 
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Fig. 54. Distribution of the invariant top quark pair mass 
reconstructed with the Matrix Element technique in 680 pb~^ 
of CDF data using semileptonic events [262| . 
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Fig. 55. Expected and observed limits on ax ■ B{X — )■ tt) 
obtained by CDF in 680 pb~^ using the Matrix Element tech- 
nique in Icpton plus jets events 
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Fig. 56. Expected and observed upper limits on ax-B{X — >■ ti) 
obtained by CDF with a Matrix Element technique in 2.8 fb~^ 
of data using all hadronic events [25^ . 



pected distribution compared to the observed data is 
shown in Fig. 1541 

The possible contribution from a resonant top quark 
pair production is computed using Bayesian statistics. The 
posterior probability density is build from a likelihood 
that implements Poissonian expectations in each bin. The 
parameters of the Poisson distribution are smeared with 
Gaussians according to the systematic uncertainties. Fi- 
nally, the posterior probability density is convoluted with 
a flat prior in az' ■ B{Z' — >■ tt). The uncertainties of the 
Standard Model top quark pair production cross-section, 
the jet energy scale and the variation of initial and final 
state gluon radiation have the largest impact on the re- 
sulting limits. 

The events observed by CDF in 680 pb^^ of data show 
no evidence for resonant top quark pair production and 
upper limits are derived for az' ■ B{Z' — > tt) as shown 
in Fig. 1551 A comparison to the leptophobic Topcolor as- 
sisted Technicolor model yield an exclusion of this model 
for Mz' < 725 GeV at 95%C.L. Using additional assump- 
tions about the kinematics through the matrix element 
thus allows to exclude slightly higher Z' boson masses de- 
spite using less data. 

The Matrix Element plus Template method was also 
applied in the all hadronic decay channel using 2.8 fb~^ 
of CDF data with six or seven jets [263] . Top quark pair 
events are enriched with a neural net event selection and 
b jet identification. The dominant multijet background is 
described with a data driven method from events without 
the b identification requirement. 

No evidence for resonant top quark pair production is 
found. CDF computed upper limits on the resonant pro- 
duction cross-section time branching fraction, ax ■B{X — >■ 
tt), as shown in Fig. [56] In the leptophobic Topcolor 
assisted Technicolor model resonance masses of Mz' < 
805 GeV are excluded at 95%C.L. 



Dynamical Likelihood Method with Massive Gluon in- 
terpretation Another analysis of CDF is based on the dy- 
namical likelihood method (DLM, see also Section [23]3]). 
It investigates 1.9 fb^^ of data with an isolated lepton, 
missing transverse energy and exactly four reconstructed 



jets [264j . In contrast to the previously described analysis 
the invariant top quark pair mass is reconstructed with- 
out using the production part of the matrix element in 
the construction of the probability densities in Eq. ((89|) . 
This avoids a bias towards the Standard Model produc- 
tion mechanism. The resulting distribution and the corre- 
sponding Standard Model expectation is shown in Fig. [571 

The distribution is then used to search for a new 
colour-octet particle, called a massive gluon. An unbinned 
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Fig. 58. Limits on the coupling strength of a massive gluon, G, 
contribution deduced for various masses and two widths i264i. 



likelihood fit based on the production matrix elements 
with and without massive gluon contribution is used to 
extract the possible coupling strengths of such a massive 
gluon contributing to the top quark pair production. The 
likelihood is computed for various masses and widths of 
the massive gluon. The systematic uncertainties are incor- 
porated in the likelihood calculation. Jet energy scale and 
top mass uncertainties are the largest contribution to the 
total uncertainty on the fitted coupling, A. 

The observed data agree with the Standard Model ex- 
pectation within ^ 1.7(7. This agreement is cross-checked 
by reconstructing the top quark pT-distribution, which is 
also found to be in agreement with the Standard Model 
expectation. Limits on the possible coupling strength of 
a massive gluon, G, contributing to the top quark pair 
production are set at 95% C.L. for various values of the 
width, Fc as function of the mass, Mg- Figure [551 shows 
the expected and observed limits for two choices of the 
massive gluon width. 



D0 

D0 investigated the invariant mass distribution of top 



quark pairs in up to 3.6 fb^ of lepton plus jets events [2651 
[266,267 . The event selection requires an isolated lepton, 
transverse missing momentum and at least three jets. At 
least one of the jets needs to be identified as b jet. Signal 
simulation is created for various resonance masses between 
350 and 1000 GeV. Separated resonance samples were gen- 
erated with couplings proportional to Standard Model Z 
boson couplings, with pure vector couplings and with pure 
axial-vector couplings. The width of the resonances was 
chosen to be 1.2% of their mass, which is much smaller 
than the detector resolution. 

The top quark pair invariant mass, Mtt, is recon- 
structed directly from the reconstructed physics objects. 
A constrained kinematic fit is not applied. Instead the mo- 
mentum of the neutrino is reconstructed from the trans- 
verse missing energy, fx, which is identified with the 
transverse momentum of the neutrino and by solving 
^w ~ (P^ + Pi^)^ -^'-"" *^^^ z-component of the neutrino 
momentum, pi and Pi, are the four-momenta of the lepton 
and the neutrino, respectively. 

The tt invariant mass can then be computed without 
any assumptions about a jet-parton assignment that is 



needed in constrained fits. Compared to the constrained 
fit reconstruction, applied in an earlier analysis, this gives 
better performance at high resonance masses and in addi- 
tion allows the inclusion of ^-1-3 jets events. The expected 
distribution of Standard Model processes and the mea- 
sured data is shown in Fig. 1591 For comparison a resonance 
with a mass of 650 GeV is shown at the cross-section ex- 
pected in the Topcolor assisted Technicolor model used 
for reference. 

Cross-sections for resonant production are evaluated 
with Bayesian statistics using a non-zero flat prior (for 
positive values) of the resonant top quark pair cross- 
section time branching fraction, ax-B{X — )• it). A Poisson 
distribution is assumed for the number of events observed 
in each bin of the likelihood. The prior for the combined 
signal acceptance and background yields is a multivariate 
Gaussian with uncertainties and correlations described by 
a covariance matrix. The measured ax • B(X — >■ tt) corre- 
spond to the maximum of the Bayesian posterior proba- 
bility density, limits are set at the point where the integral 
of the posterior probability density from zero reaches 95% 
of its total. Expected limits are obtained by applying the 
procedure when assuming that the observed result cor- 
responded to the Standard Model expectation. The limits 
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Fig. 59. Expected and observed ii invariant mass distribution 
for the combined (a) ^ -I- 3 jets and (b) ^ -I- 4 or more jets chan- 
nels, with at least one identified b jet. Superimposed as white 
area is the expected signal for a Topcolor assisted Technicolor 
Z' boson with Mz' = 650 GeV [267] . 
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Fig. 60. Left: Expected limits ax-B{X — )■ ii) vs. tlie assumed resonance mass for 0.9 fb~^. From bottom to top tlie lines represent 
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except mt and complete systematics 268 . Right: Expected upper limits with complete systematics in 3.6 fb~^ compared to the 
observed cross-section and exclusion limits at 95%CL 12671. 



obtained for the Z-like, veetor and axial-veetor resonances 
are found to agree with each other, thus these limits are 
valid for a general (colour neutral) narrow resonance. 

These expected limits were used to optimise major 
analysis cuts and the 6-tag working point. In Fig. l60l (left) 
the expected limits are used to visualise the effect of 
the various systematics by including one after another. 
The lowest curve corresponds to a purely statistical limit. 
Adding the jet energy uncertainty shows that this uncer- 
tainty mainly contributes at medium resonance masses. 
The various object identification efficiencies and the lu- 
minosity are added. They essentially scale like the back- 
ground shape. Finally the effect of the top quark mass 
is included and it is most important at low resonance 
masses 



In 3.6 fb~^ of D0 data the observed cross-sections are 
close to zero for all considered resonance masses, as shown 
in Fig. |60] (right). The largest deviation (around 650 GeV) 
corresponds to a little more than 1.5 standard deviations. 
Thus limits are set on the ax ■ B{X —>■ tt) as function 
of the assumed resonance mass, Mx- The excluded values 
range from about 1 pb for low mass resonances to less then 
0.2 pb for the highest considered resonance mass of 1 TeV. 
The benchmark Topcolor assisted Technicolor model can 
be excluded for resonance masses of Mz' < 820 GeV at 
95%CL. 



4.5 Admixture of Stop Quarks 

A final very fundamental question that may be asked in 
the context of top quark physics beyond the Standard 
Model is, whether the events that are considered to be 
top quarks actually are all top quarks or whether some 
additional unknown new particle is hiding in the selected 
data. The top quark's supersymmetric partners, the stop 
quarks ti and t2, are possible candidates in such a sce- 



4.5.1 D0, Lepton plus Jets 

The stop quark decay modes to neutralino and top quark, 
Xit, or through chargino and b quark, Xi b, both yield 
a final state with a neutralino, a b quark and a W bo- 
son, XibW. The neutralino is the lightest supersymmetric 
particle in many models and is stable if i?-parity is con- 
served. Then it escapes the detector and the experimental 
signature of stop quark pair production differs from that 
of semileptonic top quark pair events only by the addi- 
tional contribution to the missing transverse energy from 
the neutralino. 

D0 has searched for a contribution of such stop quark 
pair production in the semileptonic channel in data with 
0.9 fb~^ |2691l270j . Semileptonic events were selected fol- 
lowing the corresponding ti cross-section analysis by look- 
ing for isolated leptons (e and /i), missing transverse en- 
ergy and four or more jets. At least one of the jets was 
required to be identified as 6-jet using D0's neural net- 
work algorithm. 

To describe the Standard Model expectation a mix- 
ture of data and simulation is employed. The description 
of top quark pair production (and of further minor back- 
grounds) is taken fully from simulation normalised to the 
corresponding theoretical cross-sections. For VF-j-jets the 
kinematics is taken from simulation, but the normalisation 
is taken from data. Multijet background is fully estimated 
from data. As signal the lighter of the two stop quarks, 
<i, is considered. Production of iiii is simulated for var- 
ious combinations of stop quark and chargino masses, 
mj^ , m~± . For the sake of this analysis the stop quark 
mass was assumed to be less or equal to the top quark 
mass. The neutralino mass m^o = 50 GeV was chosen to 
be close to the experimental limit. 

To detect a possible contribution of stop quark pairs 
the differences between stop quark pair events and Stan- 
dard Model top quark pair production kinematic variables 
are combined into a likelihood, C = Pstop/ (-Pstop + Psm)- 
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Fig. 61. Left: Expected and observed distribution for the reconstructed top quark mass. Right: Expected likelihood distribution 
for Standard Model and signal compared to data. In both plots only events with four or more jets are shown. The stop quark 
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Fig. 62. Expected and observed limits on the stop quark pair production cross-section compared to the expectation in the 
various MSSM parameter sets [270| . 



The kinematic variables considered include the transverse 
momentum of the (leading) b jet, distances between lead- 
ing 6 jet and lepton or leading other jet. Additional vari- 
ables were reconstructed by applying a constrained fit. 
In this fit reconstructed physics objects (lepton, missing 
transverse energy and jets) are assigned to the decay prod- 
ucts of an assumed semileptonic top quark pair event and 
the measured quantities are allowed to vary within their 
experimental resolution to fulfil additional constraints. It 
was required that the W boson mass is consistent with the 
invariant mass of the jets assigned to the two light quarks 
as well as with the invariant mass of the lepton with the 
neutrino. The masses of the reconstructed top quark pairs 
were constrained to be equal. Of the possible jet parton 
assignments only the one with the best x^ was used. From 
the constrained fits observables the angle between the b- 
quarks and the beam axis in the bb rest frame, the bb in- 
variant mass, the distances between the 6's and the same- 
side or opposite-side W bosons and the reconstructed top 
quark mass are considered. The likelihood was derived for 
each rrif , "^ -± combination separately and the selection 
of variables used has been optimised each time. Figure I5T] 



shows the separation power of the reconstructed top quark 
mass and the full likelihood for the case of m^ =175 GeV, 
TO -± = 135 GeV and the comparison to the observed data. 



To determine limits on the possible contribution 
of stop quark pair production in the selected channel 
Bayesian statistics is employed using a non-zero flat prior 
for positive values of the stop quark pair cross-section. A 
Poisson distribution is assumed for the number of events 
observed in each bin of the likelihood. The prior for the 
combined signal acceptance and background yields is a 
multivariate Gaussian with uncertainties and correlations 
described by a covariance matrix. The systematic uncer- 
tainty is dominated by the uncertainties on the theoretical 
cross-section of top quark pair production, on the selec- 
tion efficiencies and the luminosity determination. Fig- 
ure [52] shows the expected and observed limits on the 
stop quark pair production cross-section compared to the 
MSSM prediction for various values of m^^ and to -± . The 
theoretically expected stop quark signal cross-section in 
the MSSM is smaller than the experimental limits for all 
parameter points considered. 
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4.5.2 CDF, Dilepton 

CDF has searched for a contribution of stop quarks in the 
dilepton channel using upto 2.7 fb~ of data [271) . The 
dilepton event signature was chosen to cover chargino de- 
cay modes to ^ -I- 1' that do not involve an intermediate W 
boson and thus may not have a corresponding hadronic 
decay to build semileptonic events. 



Reconstructed Stop Mass, B-Tagged Channel 
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CDF collected data using an inclusive high-p^ trigger 
and selects events with 2 leptons, one of which needs to be 
isolated from calorimeter energies not associated to that 
lepton. The events are required to have missing transverse 
energy and at least 2 jets. A Z boson veto is applied for ee 
and pill events. To suppress the leading background, Stan- 
dard Model top quark pair events, from the selected event 
a cut in the plane of Ht vs. A plane is applied, where A 
is the product of the azimuthal angles between the lead- 
ing jets and the two leptons: A — Z\0(jet]^,jet2)Zi0(i?i,^2)- 
This cut reduces top quark pair production by a factor of 
2, but reduces stop quark by approximately 12% only. 

The Standard Model background expectation is mod- 
elled using simulation and control data. Simulation of top 
quark pair production and other minor backgrounds are 
normalised to their NLO cross-section. Z-|-jets samples 
are normalised to control data of low missing transverse 
energy near the Z-pole separately for events with and 
without 6-tags. To model events with faked leptons that 
may stem from other top quark pair decay channels, from 
Ty-|-jets or from multijet events, parametrised lepton fake 
rates are derived from a large sample of generic jets. These 
fake rates are applied to events with lepton plus electron 
or muon like events to find the contribution of fakes in the 
signal region. 

To describe signal events stop quark pair produc- 
tion is simulated with various combinations of neutralino, 
chargino and stop quark masses lighter than the top quark 
mass. Generated samples are interpolated through a tem- 
plate morphing technique to obtain any combination of 
masses within the generated range. 

The reconstructed stop quark mass is used to distin- 
guish a stop quark signal from Standard Model back- 
grounds including top quark pair production. The stop 
quark mass, mi , is determined following an extention of 
the dilepton neutrino weighting technique. 

6-jets are assigned to their proper lepton based on jet- 
lepton invariant mass quantities. A correct assignment is 
reached in 85% to 95% of the cases with both 6-jets being 
the leading 2 jets. Neutralino and neutrino are considered 
as a single, though massive, pseudo particle. For given (j) 
directions of the pseudo particles the particle momenta 
are determined with a fit to the measured quantities us- 
ing constraints on the assumed pseudo particle mass, the 
assumed chargino mass and the equality of the two stop 
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Fig. 63. Distribution of reconstructed stop quark masses in 
data and simulation for two choices of the stop quark, chargino 
and neutralino masses 12711. 



quark masses. The reconstructed stop quark mass is com- 
puted as weighted average of the fitted stop quark masses, 
where the average is computed over all values of (j), with 
weights of e~-^ . The expected and observed distributions 
for two choices of parameters are shown in Fig. [521 

The combination of reconstructed stop quark mass 
templates from the signal and the various background 
components is fitted to data. Systematic uncertainties en- 
ter the fit through nuisance parameters, signal and back- 
ground contributions are allowed to vary within their rate 
uncertainties and the shape may vary according to CDFs 
morphing technique. The ratio of the likelihoods is used to 
do the limit-setting according to the CLg technique |227l 



Depending on the dilepton branching ratio limits are 
set in the stop quark vs. neutralino mass plane. Figure [Ml 
shows the results for two choices of parameters. These 
limits are derived using only very few assumptions, these 
are (a) x? is the LSP and q, I and v are heavy, (b) ra^ < 
rrit and (c) m^± < mi — rrif,. Thus the limits arc valid 
over a large range of SUSY parameter space. 
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Fig. 64. CDF observed 95%CL limits in the stop mass vs. chargino mass plane. The left plot shows limits for a chargino mass 
of 105.8 GeV, the right for 125.8 GeV gn]- 



4.6 Heavy Top-like Quark, t' 

Another class of hypothetical particles that may hide in 
samples usually considered as top quarks are new heavy 
quarks, in this context usually called t' quark. These new 
particles a considered to decay to Wq and thus show a 
signature very similar to that of top quarks. 

Heavy top- like quarks appear in a large number of new 
physics models: A fourth generation of fermions [2 72) . Lit- 
tle Higgs models [2 73) and more named in |274| . Strong 
bounds are placed on such models by electroweak precision 
data, but for special parameters the effects of the fourth 
generation particles on electroweak observables compen- 
sate. Among other settings a small mass splitting between 
the fourth w-type quark, t' , and its isospin partner, 6', is 
preferred, i.e. Wb' -|- Mw > ruf |272| . Especially, when the 
new top-like quark is very heavy, it should be distinguish- 
able from Standard Model top production in kinematic 
observables. A search for such heavy top-like quarks has 
been performed by both Tevatron experiments. 



CDF 



The CDF collaboration has repeatedly analysed their sam- 
ples of lepton plus jets events to search for a new top-like 
quark. The published result uses 760 pb~^ |274) and the 



preliminary updated results 4.6 fb^^ of data [275] . The 
analyses consider pair production of a new quark heavier 
than the top quark and with a subsequent decay to Wq. 
The event selection requires exactly one isolated lepton, 
large missing transverse energy and at least four energetic 
jets. 

The dominant Standard Model processes that pass this 
selection are top quark pair production which is simu- 
lated using Pythia, M^-|-jets events which are simulated 
using Alpgen-I-Herwig or Pythia and normalised to 
data and multijet events which are modelled from data 
with reversed lepton identification. Minor backgrounds 



like Z-t-jets and diboson events are considered to be de- 
scribed by the M^-|-jets simulation. The simulation of t'F 
signal is performed using Pythia. 

As the t' quark decay chain is the same as the top 
quark decay chain a constrained fit to the kinematic prop- 
erties is performed. The momenta of the quarks and lep- 
tons from the t^'^ quark and the following W boson decays 
are fitted to the observed transverse momenta. The con- 
straints require that the W boson decay products form 
the nominal W mass and that the decay products of the 
t^'' quark yield the same mass on the hadronic and the 
leptonic side. Of the 12 different jet-parton assignments 
CDF chooses the one with the best x^ of the fit. The cor- 
responding t^') mass from the fit, Mrcco, is used as one ob- 
servable to separate signal from background. The second 
observable is the total transverse energy, Ht, i-e. the sum 
of transverse energies of the observed jets, the lepton and 
the missing transverse energy. The choices explicitly avoid 
imposing 6-quark tagging requirements. The expected and 
observed individual distributions are shown in Fig. 1651 

The signal and background shapes in the two dimen- 
sional M^cco-Ht plane are used to construct a likelihood 
for the observed data as function of the assumed t' quark 
cross-section, at' ■ Then Baysian statistics is employed to 
compute expected and observed limits on ut' that are 
shown in Fig. 1661 

Systematic uncertainties are implemented through nui- 
sance parameters that are constrained in the fit with a 
Gaussian function to their nominal value within their ex- 
pected uncertainty. The jet energy scale is named as one 
of the largest uncertainties. In addition uncertainties on 
the Q^ scale used in simulating W^-f jets, on initial state 
and final state radiation, on the multijet background de- 
termination, the integrated luminosity, the lepton identi- 
fication, the parton density functions and the expected t' 
quark cross-section as function of the t' quark mass are 
considered. 

The limits on t' quark pair production cross-section. 
It' , determined in this search for a new top-like heavy 
quark are compared to the theoretical prediction 131' . As- 
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Fig. 65. Expected and observed distribution of tlie the total transverse energy, Ht, and the reconstructed r'' mass in 4.6 fb ^ 
of CDF data including a hypothetical signal at rrit' — 450 GeV |275j . 
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Fig. 66. Expected and observed limits on the cross-section of a 
new top- like quark determined from 4.6 fb~^ of CDF data 275 . 
The shaded bands show the expected one and two sigma vari- 
ation on the expected upper limit. 



suming t' — )■ Wq CDF concludes that a t' pair production 
can be excluded for m^ < 335 GeV at 95%CL. However, 
for masses of rrit' ~ 350 GeV the observed limit is viforse 
than the expectation by about two standard deviations, 
which indicates a surplus of data for that range. 



D0 

The D0 collaboration has searched for a heavy top like 
quark with a decay to a W-boson and a quark using 
4.3 fb~^ of data |276| . Events are selected which show 
an isolated lepton, large missing transverse energy and 
at least four jets. 

The dominating backgrounds from Standard Model 
processes are top quark pair and VF-f-jets production. 



which are simulated using Alpgen+Pythia. The tt sim- 
ulation and additional smaller backgrounds from diboson 
and ^-boson production are normalised to the expected 
Standard Model cross-sections. The normalisation of the 
M^-|-jets contribution is determined below. The contribu- 
tion of multijet background faking a lepton is derived from 
a sample with loosened lepton selection. Signal simulation 
is obtained from Pythia for 13 different t' mass values 
using a fixed decay width of 10 GeV. 

Also D0 uses a reconstructed i'-quark invariant mass 
and the scalar sum of transverse momenta to distinguish 
the t'-quark signal from the Standard Model backgrounds. 
The t'-quark invariant mass is derived in a full event re- 
construction assuming the expected decay t'i' -^ ivbqqh. 
In this reconstruction the momenta of the assumed decay 
products are fitted to the reconstructed measured objects 
applying constraints that enforce the mass of the interme- 
diate M^-bosons and equality of the initial t' and t' masses. 
All possible associations of the four leading jets in px to 
the quarks are considered. To select the best association in 
addition to the fit x^ a Ax^ term is computed that prefers 
low transverse momenta for the reconstructed t'-quark. 
The t'-mass reconstructed, mfit, with the association that 
minimises x^ + ^X^ is taken for further analysis. 

Two dimensional distributions of Ht vs. mfit are used 
to fit the data compositions with free Ty-|-jets and signal 
normalisations, c.f. Fig. [57] This is done for all 13 t'-quark 
mass hypotheses and for a background only hypothesis. 
D0 uses the CL^ method [2?7[l^^ with a Poisson likeli- 
hood as the test statistics to determine the cross-section 
limits. D0 includes several systematic uncertainties af- 
fecting the normalisation of the fixed background contri- 
butions as well as shape changing effects in the limit calcu- 
lation. Of the normalisation uncertainties the luminosity 
uncertainty gives the largest single contribution. 
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Fig. 67. Expected and observed distribution of the the total transverse energy, Ht, and the reconstructed t''' mass in 4.3 fb ^ 
of D0 data including a hypothetical signal at m^/ = 300 GeV [276| . 
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Fig. 68. Expected and observed limits for the production of 
a i'-quark obtained in 4.3 fb~^ of D0 data. The shaded bands 
indicate the one and two sigma deviations from the expected 
limit 



The limits on a t'-quark production cross-section ob- 
served by D0 in 4.3 fb^^ [276 and the corresponding ex- 
pectations are shown in Fig. |BS] as function of the hypo- 
thetical t' mass. By comparing to the theoretical predic- 
tion [31j D0 excludes heavy top like quarks with masses 
less than 296 GeV at 95% CL. Over almost the full range 
the observed limit stays behind the expected limit in- 
dicating a surplus of data in the t'-signal range. Corre- 
spondingly the expected mass exclusion of 330 GeV is not 
reached. 



4.7 Outlook to LHC 

Due to the high mass of the top quark it has been spec- 
ulated that the top quark may play a special role in par- 
ticle physics. Many models of new physics involve new 
particles that may occur in top quark production or de- 
cay. The Tevatron experiments have investigated a wide 
range of options for such new particles in top quark events 
and also looked for particles with signatures very similar 
to that of top quarks. So far no significant deviation has 
been found. 



In the near future the LHC is scheduled to deliver at 
least 1 fb"^ at yi = 7TeV. With such a dataset the LHC 
experiments will have roughly twice as many top quark 
pairs and three times as many single top quarks as the 
Tevatron with the expected 10 fb^ . The benefit of this 
data, however, depends on the process in which the new 
particle appears. 

Searches for new particles that occur in the top quark 
decay, like e.g. a charged Higgs boson, will benefit directly 
from the increased dataset and yield results competitive 
to those of the Tevatron. Searches for new particles that 
occur in the production like, Z' or W' bosons or a t' quark, 
will benefit from the larger centre-of-mass energy. The 
cross-sections of such processes, however, scale with the 
corresponding parton luminosity, which increases much 
more for the gluon than for the quark luminosity. Particles 
like the Z' boson that can by produced only by quark- 
antiquark annihilation thus will have less signal events in 
1 fb~ than the Tevatron will have in the expected 10 fb~ . 
Such analyses will still extend the Tevatron results to the 
phase space areas opened with the higher centre-of-mass 
energy. Finally, the search for a boson produced in the 
single top quark s-channel, like the W , will suffer from 
the comparably moderate increase of the corresponding 
cross-section and thus require more statistics at the LHC 
to become competitive to the Tevatron results. 

With the design LHC at i/s = 14 TeV the larger pro- 
duction cross-section and the enlarged phase space will 
again extend the results significantly. An ATLAS study 
assuming collisions at ^/s — 14 TeV and a luminosity of 
lfb~^ 1 144) has considered the potential for discovering 
resonant top quark pair production through a narrow Z' 
boson. A significant degradation of the selection efficiency 
is expected at high top pair invariant masses because the 
top quark decay products join into the same jet more and 
more often. Given the performance visible in the first re- 
sults shown at conferences from LHC experiments searches 
for new physics in top quark events will greatly advance 
at each new energy reached by the LHC after collecting a 
moderate luminosity of e.g. 1 fb~ . 
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5 Conclusions 

More than fifteen years after the discovery of the top quark 
the experimental verification of the properties expected in 
the Standard Model has reached a first level of maturity. 
With the luminosity delivered by the Tevatron accelerator 
in its Run II to the CDF and D0 experiments many top 
quark properties have been challenged and contributions 
for new particles have been searched for. 

The top quark mass measurements, with the elabo- 
rate statistical methods of the Matrix Element technique, 
have now reached a statistical precision of less than 1 GeV 
and experimental systematics of 1 GeV making the top 
quark mass the most precisely known mass of all quarks. 
This precision significantly exceeds the precision goals of 
about 3 GeV of the Tevatron Run II programme set for 
2fb~^ [4711277] . Despite the low branching fraction of the 
dileptonic channel, even results in this channel alone now 
have achieved an experimental precision exceeding this 
goal. In order to reach these small uncertainties, it was 
important to constrain the jet energy scale, leading to the 
dominating uncertainty, in-situ to data. Unfortunately, 
the experimental precision on the top quark mass is cur- 
rently not matched by a corresponding theoretical under- 
standing. The uncertainty of the mass definition used in 
the simulations that the experiments apply to calibrate 
their measurements is not known better than to the or- 
der of 1 GeV. This complicates the comparison of the top 
quark mass results, e.g. with other electroweak precision 
data, and is a fundamental problem also for the LHC pro- 
gramme. Discussions on these issues have started between 
the experimental, the generator and the theoretical com- 
munity in order to collect the information needed to over- 
come this issue. The development of methods that pre- 
cisely determine the mass in well defined definitions is an 
important problem in current top quark physics. 

Many interaction properties have been challenged by 
the Tevatron experiment in the recent years. Production 
properties such as the quark and gluon induced production 
rates, the electrical charge and the various decay proper- 
ties did not show significant deviation from the expec- 
tation and thus confirm that we observe the top quark 
expected in the SM. The V — A structure of the weak 
top quark decay has been tested in the W boson helicity 
measurements to the 5% level and weak decays through 
flavour changing neutral currents are constrained to be 
below the 4% level. The largest deviation (3.4cr) from the 
SM was found in the forward backward charge asymmetry 
at high rritf. While the sheer number of tests documents 
a great progress the precision of these measurements is 
generally limited by statistics. Searches for new particles 
have been performed in events with top quark like sig- 
natures. Neither specific searches for supersymmetric top 
quark partners, for charged Higgs bosons nor searches for 
W bosons, t' quarks or generic resonances found any sig- 
nificant deviations from the SM expectations. Also these 
studies are limited by the statistics available so far. 

The Tevatron is planned to be shutdown in autumn 
2011. Until then it is expected to increase the total in- 
tegrated luminosity for the CDF and D0 experiments to 



more than 10 fb~ . The analyses of the Tevatron take ad- 
vantage of several years of optimising the reconstruction 
of the recorded events. Thus further updates of existing 
analyses, but also investigation of yet untested proper- 
ties can be expected. In addition the LHC has started 
to operate at y/s = 7TeV well above that of the Teva- 
tron centre-of-mass energy. It is expected to deliver an 
integrated luminosity of at least lfb~^ throughout 2011. 
Due to the increased collision energy at the LHC the top 
production cross-sections are significantly higher than at 
the Tevatron. With 1 fb~^ the LHC experiments will have 
collected about twice as many top quark pairs and three 
times as many single top quarks than the Tevatron exper- 
iments at 10 fb^^. 

Given the impressive performance of the LHC experi- 
ments shown at recent conferences it is natural to assume 
that the systematic uncertainties can be controlled com- 
parably well as at the Tevatron experiments, if not better. 
We will thus see a plethora of competitive results from the 
LHC experiments. Clearly some searches profit from the 
increased phase space and will be dominated by the LHC. 
There are also results where the Tevatron has an advan- 
tage, like the charge asymmetry and the single top quark 
W search, however, most results will be comparable. Also 
at the LHC many will be statistically limited which makes 
it useful to foresee a combination of the results. 

After an update to the design LHC with proton-proton 
collisions at -^/i = 14TeV the field will be taken over by 
the LHC experiments. With a top quark pair cross-section 
which is 100 times higher than at the Tevatron and domi- 
nant backgrounds like M^-f jets and Z-|-jets increasing only 
by factors of about 10 the LHC will become a real top 
quark factory. Again many searches for new particles will 
profit from the enlarged phase space. Measurements of 
interaction properties and searches for new particles in 
top quark decays will become limited by systematics after 
only a short running at design luminosity. At the same 
time with the huge expected statistics new types of mea- 
surement methods on rare subsamples will become feasi- 
ble. For the top quark mass methods using leptonic J/ip 
decays or using events with very high top quark momen- 
tum will have very different experimental and theoretical 
uncertainties. Such alternative methods with different sys- 
tematic uncertainties will contribute to the combination 
of results. And may even help to solve the puzzle of the 
meaning of the top quark mass. 
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